Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,687,678 Δ1 bp 5.4% coding (1344/2061 nt) pilJ ← twitching motility protein PilJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,687,6780C. 5.4% 356.3 / 13.9 112coding (1344/2061 nt)pilJtwitching motility protein PilJ
Reads supporting (aligned to +/- strand):  ref base C (55/51);  new base . (2/4);  total (57/55)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.53e-01

Only 100 of 119 total aligned reads displayed.

CTTGGCGGTGTCCTGGATCTGCTCGCGAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTGCACTGCACTGTGGTTGATGGTGGCGACCAGGTC  >  NC_002947/5687531‑5687827
                                                                                                                                                    |                                                                                                                                                     
cTTGGCGGTGTCCTGGATCTGCTCGCGAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGAtcgactc                                                                                                                                                     <  1:1888814/150‑1 (MQ=255)
      ggTGTCCTGGATCTGCTCGCGAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCTGTCGAGTCGATCGAcctg                                                                                                                                               <  2:398792/150‑3 (MQ=255)
      ggTGTCCTGGATCTGCTCGCGAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCAt                                                                                                                                               <  1:1235837/150‑1 (MQ=255)
       gtgtCCTGGATCTGCTCGCGAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTTGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGAGTCGATCGAccggg                                                                                                                                              <  2:2179212/150‑4 (MQ=255)
              ggATCTGCTCGCGAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTcaccca                                                                                                                                       >  1:1847316/1‑150 (MQ=255)
                   tGCTCGCGAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGcc                                                                                                                                  >  2:834208/1‑150 (MQ=255)
                   tGCTCGCGAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGcc                                                                                                                                  <  2:1527140/150‑1 (MQ=255)
                     cTCGCGAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTc                                                                                                                                >  2:1194302/1‑150 (MQ=255)
                        gcgAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGa                                                                                                                             <  2:1920060/150‑1 (MQ=255)
                        gcgAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGa                                                                                                                             <  1:2019109/150‑1 (MQ=255)
                        gcgAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGa                                                                                                                             >  1:382055/1‑150 (MQ=255)
                          gAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAgg                                                                                                                           <  2:917203/150‑1 (MQ=255)
                           aaTGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGc                                                                                                                          >  2:2342927/1‑150 (MQ=255)
                                  tCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTg                                                                                                                   >  2:1079130/1‑150 (MQ=255)
                                     aTGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATc                                                                                                                <  2:684342/150‑1 (MQ=255)
                                         cGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCtgtg                                                                                                            <  1:103960/150‑1 (MQ=255)
                                             gATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGc                                                                                                        >  1:915660/1‑150 (MQ=255)
                                             gATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGc                                                                                                        <  1:1079130/150‑1 (MQ=255)
                                              aTGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGAGTCGATCGAccgg                                                                                                                                               >  1:1756378/1‑108 (MQ=255)
                                                    ttgtGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGc                                                                                                 >  2:385366/1‑150 (MQ=255)
                                                         ccgcaccCCGCTGCGCTTGTTGGCAATGGCCACCGTACGCTCGGCCACCTTGGCCGATTGGTAGGCAGGCGCTGAAACCCGGCTGATCGGCTCGACCATGTCCCCCACCGCCGCGGAGGCCTCGCTGGTCTGTGCAGCCTGGTGCTCggg                                                                                            <  2:707110/147‑2 (MQ=255)
                                                         caccacCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGa                                                                                            <  2:382055/150‑1 (MQ=255)
                                                         caccacCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACATTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGa                                                                                            <  2:1611161/150‑1 (MQ=255)
                                                          accaccTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCCACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAg                                                                                           <  1:1500205/150‑1 (MQ=255)
                                                           ccaccTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCGTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAgg                                                                                          >  2:1481072/1‑150 (MQ=255)
                                                            caccTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGc                                                                                         >  2:1734582/1‑150 (MQ=255)
                                                            caccTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCCTTCCAGTTTCGCGACCATGTCACCCACGGCCTAGGAGGCACAGCTGATCTGTGCAGCCTGGTGCTCGGAGGc                                                                                         >  1:1853962/1‑150 (MQ=255)
                                                              ccTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCt                                                                                       >  2:1202814/1‑150 (MQ=255)
                                                               cTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCtt                                                                                      >  1:684893/1‑150 (MQ=255)
                                                               cTCGTTGCCCTTGTTGGCAATGGCCACCGAACGATGGGGGAATCTGGACGACGGAGATGCCTTCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCtt                                                                                      >  1:46360/1‑150 (MQ=255)
                                                                tCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTg                                                                                     >  1:685971/1‑150 (MQ=255)
                                                                 cGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCGCCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGTTGCTCGGAGGCCTTgg                                                                                    <  2:215918/150‑1 (MQ=255)
                                                                  gTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGc                                                                                   <  2:1847316/150‑1 (MQ=255)
                                                                    tGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCa                                                                                 <  2:1215449/150‑1 (MQ=255)
                                                                     gCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCCGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTAATAGGCCTTGGCCAg                                                                                <  1:1181012/150‑1 (MQ=255)
                                                                       ccTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCa                                                                                                                                                >  1:1046127/1‑84 (MQ=255)
                                                                       ccTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCa                                                                                                                                                <  2:1046127/84‑1 (MQ=255)
                                                                       ccTTGTTAGAAATGGCCACCGAACGCTTGGCCACCGTGGCCGATTCGGAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCCCCCAACGCCGAAGAGGCCGCACGCATCGACTCAGCCGGGTGCTCGGATGCCTTGGCCAGtt                                                                              >  2:1728124/1‑150 (MQ=255)
                                                                        cTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGtt                                                                              >  1:1079683/1‑149 (MQ=255)
                                                                        cTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGtt                                                                              <  2:1079683/149‑1 (MQ=255)
                                                                        cTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTg                                                                             >  2:230591/1‑150 (MQ=255)
                                                                            ttgGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCcgc                                                                         >  2:1075759/1‑150 (MQ=255)
                                                                            ttctcAATGGCCACGGAACGCTCGGTCACCTTGGCCGGTTCGTAGACATGCGCCGAAACCCGGTTGATCGCCCCCACCATGTCACCAACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGGGCTCGGAGGCCTTGGCCAGTTGCcgc                                                                         >  2:970984/5‑150 (MQ=255)
                                                                                 aaTGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAgt                                                                    <  2:685971/150‑1 (MQ=255)
                                                                                 aaTGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAgt                                                                    >  1:1176367/1‑150 (MQ=255)
                                                                                      ccACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCt                                                               <  1:2337373/150‑1 (MQ=255)
                                                                                       cACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAANCCGGTCGGTCGACTCGACCATGTCACCCACCGCCTCGGGGGCCTCGCTGATCTGTGCAGCCTGGAGCTCGGGGGCCTTGGCCAGTTGCCGCGCAGTGTTCTg                                                              >  1:1037830/1‑150 (MQ=255)
                                                                                       cACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGAGTCGATCGAccgg                                                                                                                                               >  1:1095354/1‑67 (MQ=38)
                                                                                       cACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGAGTCGATCGAccgg                                                                                                                                               <  2:1095354/69‑3 (MQ=38)
                                                                                       cACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGGAGCCTGGTGCTCGGAGGCCTTGGCCAATTGTCGCGCAGTGTTCTg                                                              >  2:1640401/1‑150 (MQ=255)
                                                                                        aCCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCAt                                                                                                                                               >  1:408557/1‑68 (MQ=255)
                                                                                        aCCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGTCGTGGCCAGTTGCCGCGCAGTGTTCTGc                                                             <  1:1872927/150‑1 (MQ=255)
                                                                                        aCCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGc                                                             <  1:1194302/150‑1 (MQ=255)
                                                                                         ccGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGACGAGCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCg                                                            >  1:1377434/1‑150 (MQ=255)
                                                                                          cGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAATCCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCgt                                                           >  2:1748197/1‑150 (MQ=255)
                                                                                          cGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCCACACCCGGTCGAACTGCTCGACCATGTCACCACCAATCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCgt                                                           <  1:2342927/150‑1 (MQ=255)
                                                                                             aCGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTc                                                        <  1:204132/150‑1 (MQ=255)
                                                                                             aCGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTc                                                        <  1:1734582/150‑1 (MQ=255)
                                                                                               gCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGACCTGGCTGATCTGTGGAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCt                                                      >  1:151922/1‑150 (MQ=255)
                                                                                               gCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGACTCGGATCCTGCGTTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCTGCGCAGTGTTCTGCGTGTCCt                                                      >  1:760847/1‑150 (MQ=255)
                                                                                                 tCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGc                                                                                                                     <  1:1360264/85‑1 (MQ=255)
                                                                                                 tCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGc                                                                                                                     >  2:1360264/1‑85 (MQ=255)
                                                                                                   ggCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTGGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCAc                                                  >  1:1147834/1‑150 (MQ=255)
                                                                                                   ggCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGcc                                                                                                       >  2:254308/1‑97 (MQ=255)
                                                                                                   ggCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGcc                                                                                                       <  1:254308/97‑1 (MQ=255)
                                                                                                    gCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCAcg                                                 <  1:1202814/150‑1 (MQ=255)
                                                                                                      cACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGTAGCCTGGTGCTCgg                                                                                             >  1:291453/1‑104 (MQ=255)
                                                                                                      cACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGTAGCCTGGTGCTCgg                                                                                             <  2:291453/104‑1 (MQ=255)
                                                                                                      cACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCAcggc                                               >  1:2004126/1‑150 (MQ=255)
                                                                                                         cTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCtgtg                                                                                                            <  2:2194229/86‑1 (MQ=255)
                                                                                                         cTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCtgtg                                                                                                            >  1:2194229/1‑86 (MQ=255)
                                                                                                         cTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCa                                                                                                          <  1:604770/88‑1 (MQ=255)
                                                                                                         cTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCa                                                                                                          >  2:604770/1‑88 (MQ=255)
                                                                                                         cTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCAcggcggc                                            >  2:1266739/1‑150 (MQ=255)
                                                                                                             gCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGgcggcggcg                                        >  1:1736766/1‑150 (MQ=255)
                                                                                                              ccGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCTTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGa                                       >  2:1579763/1‑150 (MQ=255)
                                                                                                               cGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGAc                                      >  1:455912/1‑150 (MQ=255)
                                                                                                                      tAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACctgcact                               <  2:1777033/150‑1 (MQ=255)
                                                                                                                         gCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGa                                       >  1:1683994/1‑139 (MQ=255)
                                                                                                                         gCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGa                                       <  2:1683994/139‑1 (MQ=255)
                                                                                                                         gCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACctgcactgca                            <  2:1853962/150‑1 (MQ=255)
                                                                                                                              cgcCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTGCACTGCACtgtg                       <  2:1176367/150‑1 (MQ=255)
                                                                                                                              cgcCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTGCACTGCACtgtg                       <  1:2345303/150‑1 (MQ=255)
                                                                                                                               gcCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCt                                                                                                                    <  1:433118/56‑1 (MQ=255)
                                                                                                                               gcCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCt                                                                                                                    >  2:433118/1‑56 (MQ=255)
                                                                                                                               gcCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTGCACTGCACTGTgg                      >  1:2005777/1‑150 (MQ=255)
                                                                                                                               gcCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTGCACTGCACTGTgg                      >  1:688384/1‑150 (MQ=255)
                                                                                                                                ccGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTGCACTGCACTGTGGt                     <  2:1501238/150‑1 (MQ=255)
                                                                                                                                 cGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGTCCAGTTTGTGCGCATTGTTCGTCGTGTCCTTCACTTTGACGGCGACCTGCACTGCACTGTGGtt                    >  2:1109330/1‑150 (MQ=255)
                                                                                                                                 cGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACctgcact                               >  1:1428274/1‑139 (MQ=255)
                                                                                                                                         ggTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCCGTGTTCTGCGTGTCCTACAGCGCGTGGGCGACCGTCACTGCACGGTGGTTGAtggtgg            <  2:1835452/150‑1 (MQ=255)
                                                                                                                                            cgatcgaCTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTGCACTGCACTGTGGTTGATGGTGGCGa         <  1:385366/150‑1 (MQ=255)
                                                                                                                                             gatcgaCTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTGCACTGCACtgtg                       >  1:861818/1‑135 (MQ=255)
                                                                                                                                             gatcgaCTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTGCACTGCACtgtg                       <  2:861818/135‑1 (MQ=255)
                                                                                                                                              atcgaCTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTGCACTGCACTGTGGTTGATGGTGGCGAcc       >  1:1743869/1‑150 (MQ=255)
                                                                                                                                                 cgactcgacCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGCCCTGCACGGCGGCGGCGACCTGCACTGCACTGTGGTTGATGGTGGCGACCAg     <  1:1640401/150‑1 (MQ=255)
                                                                                                                                                   actcgacCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCCGTGTTCAGGATCGCCAGCACGGCGGAGGCCGCCGGCACTGCACTGTGGTTGATGGTCGCGACCAgg  >  2:1771168/1‑149 (MQ=255)
                                                                                                                                                   actcgacCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTGCACTGCACTGTGGTTGATGGTGGCGACCAg     >  1:87485/1‑148 (MQ=255)
                                                                                                                                                   actcgacCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTGCACTGCACTGTGGTTGATGGTGGCGACCAg     <  2:87485/148‑1 (MQ=255)
                                                                                                                                                    ctcgacCATGTCACCCACCGCCTCGGAGGCCTTGCTTATCTGTGCAGCCTGGTGCTCTGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTGCACTGCACTGTGGTTGATGGTGGCGACCAGGTc  <  2:1150182/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTTGGCGGTGTCCTGGATCTGCTCGCGAATGTTGTCCATGCCGTTGATGGTGTTGTGCACCACCTCGTTGCCCTTGTTGGCAATGGCCACCGAACGCTCGGCCACCTTGGCCGATTCGTAGGCATGCGCCGAAACCCGGTCGATCGACTCGACCATGTCACCCACCGCCTCGGAGGCCTCGCTGATCTGTGCAGCCTGGTGCTCGGAGGCCTTGGCCAGTTGCCGCGCAGTGTTCTGCGTGTCCTGCACGGCGGCGGCGACCTGCACTGCACTGTGGTTGATGGTGGCGACCAGGTC  >  NC_002947/5687531‑5687827

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: