Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 2,122,874 T→G 9.6% Q794H (CAA→CAC)  fimD_1 ← Outer membrane usher protein FimD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1432,122,8740TG 9.6% 314.1 / 2.6 134Q794H (CAA→CAC) fimD_1Outer membrane usher protein FimD
Reads supporting (aligned to +/- strand):  ref base T (26/91);  new base G (14/3);  total (40/94)
Fisher's exact test for biased strand distribution p-value = 2.33e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 135 total aligned reads displayed.

TGCACCGCGTGAAGGTTGCTAGATACAGGTCAGGGTTTGCTGGCGATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTGGGTATCGTGCAGGGTCAGTACAAGG  >  ACB143/2122742‑2123009
                                                                                                                                    |                                                                                                                                       
tGCACCGCGTGAAGGTTGCTAGATACAGGTCAGGGTTTGCTGGCGATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGGGGGTGGCCCCCCCAAACCCGCCGCGGTTgg                                                                                                                                      >  5:223535/1‑135 (MQ=255)
tGCACCGCGTGAAGGTTGCTAGATACAGGTCAGGGTTTGCTGGCGATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCCACACACCTGCAGCGGGTGc                                                                                                                                      >  1:277581/1‑136 (MQ=255)
tGCACCGCGTGAAGGTTGCTAGATACAGGTCAGGGTTTGCTGGCGATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCGGTTCGGCCCCACACCCCTGGCGCGGGTgg                                                                                                                                      >  3:308891/1‑135 (MQ=255)
tGCACCGCGTGAAGGTTGCTAGATACAGGTCAGGGTTTGCTGGCGATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCGGTTAGCCCCCACGCACCAGCCGCGGGTgg                                                                                                                                      >  4:381227/1‑135 (MQ=255)
     cgcgTGAAGGTTGCTAGATACAGGTCAGGGTTTGCTGGCGATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGGTTTTAGCCCCAACACACCCGGGGCGGGTGcggctg                                                                                                                                 >  1:47717/1‑131 (MQ=255)
     cgcgTGAAGGTTGCTAGATACAGGTCAGGGTTTGCTGGCGATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGGGGGTGGGCCCCACACACCAGCAGCCGGTGcggccg                                                                                                                                 >  5:342208/1‑131 (MQ=255)
     cgcgTGAAGGCTGCTAGATACAGGTCAGGGTTTGCTGGCGATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGGTTCTCctg                                                                                                                                 >  7:308523/1‑136 (MQ=255)
          gAAGGTTGCTAGATACAGGTCAGGGTTTGCTGGCGATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGGTGTTCGCCCCCAAACACCAGCAGCGGTTGCTCcggggggg                                                                                                                            >  1:349118/1‑132 (MQ=255)
          gAAGGTTGCTAGATACAGGTCAGGGTTTGCTGGCGATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTAGCCCCAACTCACCTGCAGCGGTTTCTCctggctgg                                                                                                                            >  8:247165/1‑136 (MQ=255)
            aGGTTGCTAGATACAGGTCAGTGTTTGCTGGCGATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTgggg                                                                                                                          <  7:183860/136‑1 (MQ=255)
            aGGTTGCTAGATACAGGTCAGGGTTTGCTGGCGATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTgggg                                                                                                                          <  2:113900/136‑1 (MQ=255)
              gTTGCTAGATACAGGTCAGGGTTTGCTGGCTATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGt                                                                                                                         <  5:95911/135‑1 (MQ=255)
              gTTGCTAGATACAGGTCAGGGTTTGCTGGCGATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGtt                                                                                                                        <  3:127635/136‑1 (MQ=255)
              gTTGCTAGATACAGGTCAGGGTTTGCTGGCGATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGtt                                                                                                                        <  6:289401/136‑1 (MQ=255)
                    aGATACAGGTCAGGGTTTGCTGGCGATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTTGGGTGTTCGCCCCAAAAAACCAGCAGCGGGTGAGACTGGCGGGGGGTGGTGAg                                                                                                                  >  1:65313/1‑136 (MQ=255)
                    aGATACAGGTCAGGGTTTGCTGGCGATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGGTGTTGGCCCCAACAAACCAGCAGCGGGTTCTACAGGGTGGGGTTGGTGAg                                                                                                                  >  1:118455/1‑136 (MQ=255)
                    aGATACAGGTCAGGGTTTGCTGGCGATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTAACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAg                                                                                                                  >  4:315276/1‑136 (MQ=255)
                     gATACAGGTCAGGGTTTGCTGGCGATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGa                                                                                                                   <  5:391213/134‑1 (MQ=255)
                     gATACAGGTCAGGGTTTGCTGGCGATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCGGTTGGCCCCCACACCCCTGCAGCGGGTGGGACAGGGGGGGGTTGGGGAGc                                                                                                                 >  4:22535/1‑136 (MQ=255)
                          aGGTCAGGGTTTGCTGGCGATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTAc                                                                                                            <  2:266421/136‑1 (MQ=255)
                          aGGTCAGGGTTTGCTGGCGATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTAc                                                                                                            <  6:329446/136‑1 (MQ=255)
                               aGGGTTTGCTGGCGATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGGTGTTGGCCCCAACTCACCTGCAGCGGTTGCTCCTGGGTGGGGGTGGTGAGCATTAACACGc                                                                                                       >  2:15672/1‑136 (MQ=255)
                               aGGGTTTGCTGGCGATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGc                                                                                                       >  8:318890/1‑136 (MQ=255)
                               aGGGTTTGCTGGCGATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCGGTTAGCCCCAAAAAACCAGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCAATAACAGGc                                                                                                       >  5:75458/1‑136 (MQ=255)
                                  gTTTGCTGGCGATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGGGGGTGGCCCCCACACACCAGCACCCGGTGGTGCGGGCGGGGGGTGGTGTGGATTACCattgtag                                                                                                    >  2:183108/1‑130 (MQ=255)
                                     tGCTGGCGATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCt                                                                                                 <  7:279557/136‑1 (MQ=255)
                                      gCTGGCGATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCt                                                                                                 <  6:87477/135‑1 (MQ=255)
                                         ggCGATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGGGGTTTGCCCCCCCACACCAGCAGCGGGTGCTCCAGGGGGGGGGTGGGTTGCATTAACCGGCCCCCCCGccc                                                                                             >  5:394198/1‑136 (MQ=255)
                                             aTAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGAcc                                                                                         <  4:51758/136‑1 (MQ=255)
                                              tAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCGGGCTGGGGTTGGTGCGCATTACCAGGCCCGCCTgcccgccc                                                                                         >  2:124486/1‑135 (MQ=255)
                                               aGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTTGTGAGCATTACCAGGCCCGCCTGCCCGACCAt                                                                                       <  1:370911/136‑1 (MQ=255)
                                               aGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACCCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCAt                                                                                       <  8:85257/136‑1 (MQ=255)
                                                gcggTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCAt                                                                                       <  3:216717/132‑1 (MQ=255)
                                                gCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCAt                                                                                       <  7:104520/135‑1 (MQ=255)
                                                gCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCAt                                                                                       <  4:221849/135‑1 (MQ=255)
                                                gCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCAt                                                                                       <  5:1950/135‑1 (MQ=255)
                                                gCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCAt                                                                                       <  6:67571/135‑1 (MQ=255)
                                                gCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCAt                                                                                       <  8:206347/135‑1 (MQ=255)
                                                gCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCAt                                                                                       <  7:311992/135‑1 (MQ=255)
                                                    tCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCa                                                                                  <  7:280449/136‑1 (MQ=255)
                                                       acgacgGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGtt                                                                               >  2:258780/1‑136 (MQ=255)
                                                                    ccAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTtctc                                                                  <  1:394096/136‑1 (MQ=255)
                                                                           aTCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAg                                                            <  1:118610/135‑1 (MQ=255)
                                                                              gagaTGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCAc                                                         <  8:248205/135‑1 (MQ=255)
                                                                                    cACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGcc                                                  <  1:218668/136‑1 (MQ=255)
                                                                                    cACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGcc                                                  <  3:378642/136‑1 (MQ=255)
                                                                                     aCCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGGTGCTCCTGGGGGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGGTGGGTTACGTTCTCGGCCCGCCCCAGGGccc                                                 >  3:100992/1‑136 (MQ=255)
                                                                                          gCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCaaa                                             <  6:355873/135‑1 (MQ=255)
                                                                                          gCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAg                                            <  7:188109/136‑1 (MQ=255)
                                                                                                ggCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGc                                      <  7:78647/136‑1 (MQ=255)
                                                                                                ggCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGc                                      <  1:4661/136‑1 (MQ=255)
                                                                                                ggCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGc                                      <  2:312229/136‑1 (MQ=255)
                                                                                                ggCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGc                                      <  3:182815/136‑1 (MQ=255)
                                                                                                  cGTTGGGGTGTTCGCCCCAACTCCCCTGCAGCGGTTGGGGCGGGGGGGGGGTGGTGGGCATTACCAGGCCCGCCTGCCCCCCCATGCCCCGTTGGGTAACGTTCTCGGCCCGCACCTGGGCCCCAAAAGGGAGccg                                    >  8:307749/1‑134 (MQ=255)
                                                                                                  cGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCgg                                    <  3:178312/136‑1 (MQ=255)
                                                                                                  cGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGGTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCgg                                    <  1:304161/136‑1 (MQ=255)
                                                                                                    ttGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGtt                                  <  4:60164/136‑1 (MQ=255)
                                                                                                    ttGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGtt                                  <  8:18107/136‑1 (MQ=255)
                                                                                                       ggCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGAc                               <  4:202843/136‑1 (MQ=255)
                                                                                                       ggCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGAc                               <  4:205519/136‑1 (MQ=255)
                                                                                                         cTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGAccc                             <  4:301520/136‑1 (MQ=255)
                                                                                                            ttCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTgg                          <  7:318890/136‑1 (MQ=255)
                                                                                                            ttCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTgg                          <  7:91846/136‑1 (MQ=255)
                                                                                                              cGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTggg                         <  3:381227/135‑1 (MQ=255)
                                                                                                              cGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTggg                         <  4:179915/135‑1 (MQ=255)
                                                                                                              cGCCCCAACGCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTggg                         <  2:265545/135‑1 (MQ=255)
                                                                                                               gCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGc                                                                                                       >  5:361619/1‑56 (MQ=255)
                                                                                                               gCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGc                                                                                                       <  6:361619/56‑1 (MQ=255)
                                                                                                               gCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTGGGTa                       >  1:27121/1‑136 (MQ=255)
                                                                                                                 cccAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTGGGTATc                     >  5:202655/1‑136 (MQ=255)
                                                                                                                  ccAACTCAGCTGCAGCGGTTGCTCCTGGCTGGGGGTGGTGAGCATTACCAGGCGCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTGTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGTTTGACCCTGGGTATCg                    <  1:300851/136‑1 (MQ=255)
                                                                                                                  ccAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCTACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTGGGTATCg                    <  5:335059/136‑1 (MQ=255)
                                                                                                                  ccAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTGGGTATCg                    <  1:231379/136‑1 (MQ=255)
                                                                                                                  ccAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTGGGTATCg                    <  6:39238/136‑1 (MQ=255)
                                                                                                                  ccAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTGGGTATCg                    <  3:68449/136‑1 (MQ=255)
                                                                                                                  ccAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTGGGTATCg                    <  4:113455/136‑1 (MQ=255)
                                                                                                                  ccAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTGGGTATCg                    <  7:125245/136‑1 (MQ=255)
                                                                                                                  ccAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTGGGTATCg                    <  4:275375/136‑1 (MQ=255)
                                                                                                                  ccAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTGGGTATCg                    <  7:220043/136‑1 (MQ=255)
                                                                                                                  ccAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTGGGTATCg                    <  2:126765/136‑1 (MQ=255)
                                                                                                                  ccAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTGGGTATCg                    <  6:113435/136‑1 (MQ=255)
                                                                                                                  ccAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTGGGTATCg                    <  7:307749/136‑1 (MQ=255)
                                                                                                                  ccAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGGTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTGGGTATCg                    <  5:264193/136‑1 (MQ=255)
                                                                                                                  ccAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGATGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTGGGTATCg                    <  5:209719/136‑1 (MQ=255)
                                                                                                                      cTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTGGGTATCGTGCa                <  5:102591/136‑1 (MQ=255)
                                                                                                                       tCACCTGCAGCGGTTGCTGCTGGCTGGTGTTTGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTGTGTATCGTGCAg               <  1:140225/136‑1 (MQ=255)
                                                                                                                         aCCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCt                                                                                                 <  5:204646/52‑1 (MQ=255)
                                                                                                                         aCCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCt                                                                                                 >  6:204646/1‑52 (MQ=255)
                                                                                                                         aCCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTGGGTATCGTGCAggg             >  5:276974/1‑136 (MQ=255)
                                                                                                                         aCCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTGGGTATCGTGCAg               <  2:318994/134‑1 (MQ=255)
                                                                                                                         aCCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTGGGTATCGTGCAg               <  5:148909/134‑1 (MQ=255)
                                                                                                                         aCCTGCAGCGGTTGCTACTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTGGGTATCGTGCAg               <  6:343038/134‑1 (MQ=255)
                                                                                                                              cAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCTGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTGGGTATCGTGCAGGGTCAg         <  8:318189/135‑1 (MQ=255)
                                                                                                                              cAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTGGGTATCGTGCAGGGTCAg         <  2:71566/135‑1 (MQ=255)
                                                                                                                              cAGCGGTTGCTCCTGGCTGGGGGTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTGGGTATCGTGCAGGGTCAg         <  2:309353/135‑1 (MQ=255)
                                                                                                                                gCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCg                                                                                            <  7:330908/50‑1 (MQ=255)
                                                                                                                                gCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCg                                                                                            >  8:330908/1‑50 (MQ=255)
                                                                                                                                 cGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTGGGTATCGTGCAGGGTCAGTAc      >  3:33928/1‑135 (MQ=255)
                                                                                                                                 cGGTTGCTCCTGGATGTGGTTGGAGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTGGGTATCGTGCAGGGTCAGTAc      >  7:353638/1‑135 (MQ=255)
                                                                                                                                    ttGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTGGGTATCGTGCAGGGTCAGTACAAgg  <  6:223535/136‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                       
TGCACCGCGTGAAGGTTGCTAGATACAGGTCAGGGTTTGCTGGCGATAGCCGTCGACGACGGGAGCGTCCAGTTCATCGAGATGCACCTGGCAGCCGGCGTTGGGCTGTTCGCCCCAACTCACCTGCAGCGGTTGCTCCTGGCTGGGGTTGGTGAGCATTACCAGGCCCGCCTGCCCGACCATGCCCAGTTGGGTACCGTTCTCGGCCAGCACCTGGGCCCCAAAGGGCAGCGGTTGACCCTGGGTATCGTGCAGGGTCAGTACAAGG  >  ACB143/2122742‑2123009

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑