Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,438,869 A→C 5.6% Q307P (CAA→CCA)  flk → putative flagella assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,438,8690AC 5.6% 325.8 / 8.0 108Q307P (CAA→CCA) flkputative flagella assembly protein
Reads supporting (aligned to +/- strand):  ref base A (52/49);  new base C (4/2);  total (56/52)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.68e-01

Only 100 of 115 total aligned reads displayed.

AATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAAG  >  NC_000913/2438722‑2439018
                                                                                                                                                   |                                                                                                                                                     
      gTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATgg                                                                                                                                               >  1:798417/1‑150 (MQ=255)
       tCTTTCTGCGCAATCATCGCTTTCTTCTGCCCCGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTGCACCGATGATGGACACGTTGCAAAATATGGc                                                                                                                                              <  2:571088/150‑1 (MQ=255)
       tCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGc                                                                                                                                              >  2:164076/1‑150 (MQ=255)
       tCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGc                                                                                                                                              <  1:875504/150‑1 (MQ=255)
       gcTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGc                                                                                                                                              <  2:1653910/149‑1 (MQ=255)
        cTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCa                                                                                                                                             <  1:1151381/150‑1 (MQ=255)
         tttCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCaa                                                                                                                                            >  1:79329/1‑150 (MQ=255)
           tCTGCGCAATCATCGCTTGCTTCTGCCCAGGCCGAGTATATTCTAAATCAGCGTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGCTGATGGACACGTTGCAAAATATGGCAAcg                                                                                                                                          <  2:302605/150‑1 (MQ=255)
             tGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAAcgcg                                                                                                                                        >  1:1132442/1‑150 (MQ=255)
             tGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAAcgcg                                                                                                                                        <  2:1513017/150‑1 (MQ=255)
             tGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAAcgcg                                                                                                                                        >  1:1645902/1‑150 (MQ=255)
             ggcgcACTCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATGAGCTTTATCAACGGCGGGTTAAAGCGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAAcgcg                                                                                                                                        <  1:300250/149‑1 (MQ=255)
              gcgcAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCCATATTTTGCAACgtgt                                                                                                                                       <  1:513567/150‑1 (MQ=255)
              gcgcAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGt                                                                                                                                       >  2:1440086/1‑150 (MQ=255)
                gcAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTTTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTcc                                                                                                                                     <  1:1381374/150‑1 (MQ=255)
                gcAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACGCGTTGCCATATTTTGCAACGTGTcc                                                                                                                                     <  2:781996/150‑1 (MQ=255)
                 cAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTccc                                                                                                                                    <  1:1404730/150‑1 (MQ=255)
                   atcatcGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCcgc                                                                                                                                  <  1:967177/150‑1 (MQ=255)
                       tcGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCtg                                                                                                                              >  2:211404/1‑150 (MQ=255)
                        cGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCtgt                                                                                                                             >  2:1679776/1‑150 (MQ=255)
                         gCTTTCTTCTGCCCAGTCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCtgtg                                                                                                                            <  2:1293121/150‑1 (MQ=255)
                         gCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCtgtg                                                                                                                            >  1:1019298/1‑150 (MQ=255)
                         gCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCtgtg                                                                                                                            >  2:1397345/1‑150 (MQ=255)
                          cTTTCTTCTGACAAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGGTATGCACCCGCTGCGTAATCCTTTGCCCCAGATGATGGACACGTTGCAAAATATATCAACGCGTCCCgcgcgttgg                                                                                                                           >  2:957422/1‑145 (MQ=255)
                           tttCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGa                                                                                                                          <  2:1557978/150‑1 (MQ=255)
                            ttcttcTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGAt                                                                                                                         >  2:985655/1‑150 (MQ=255)
                            ttcttcTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGAt                                                                                                                         >  1:1258751/1‑150 (MQ=255)
                              cttctGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATAc                                                                                                                       >  2:1570541/1‑150 (MQ=255)
                                tctGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTg                                                                                                                     >  1:1295850/1‑150 (MQ=255)
                                tctGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTg                                                                                                                     <  2:436719/150‑1 (MQ=255)
                                         gCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATTCAACCGCTTCTTAATCCTTTTCCACCGATGATGTACACGTTGCAAATTATGGCAACTCGTCAAGCTCTGTGGATACTGTTAGTctct                                                                                                            >  1:1841990/1‑147 (MQ=255)
                                         gCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGCATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCACCGCTTCACGCGCTGTGGATACTGTTAGTcgcg                                                                                                            >  2:159692/1‑150 (MQ=255)
                                          ccGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGa                                                                                                           <  2:1464055/150‑1 (MQ=255)
                                          ccGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGa                                                                                                           <  1:930482/150‑1 (MQ=255)
                                          ccGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGAGTGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACAGTTAGTCGCGa                                                                                                           <  1:523688/150‑1 (MQ=255)
                                                 tatTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCt                                                                                                    <  2:642429/150‑1 (MQ=255)
                                                   ttCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCCACCGCTGCTTAATCCTTATCCCCCGATGATGGAAACGTTGCAAAATATGGCAACGCGTCACGCGCTGTGGATACTGTTAGTCGCGATTATCCTgt                                                                                                  >  2:1276295/1‑149 (MQ=255)
                                                   ttCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGa                                                                                                  >  1:2006663/1‑150 (MQ=255)
                                                     cgtAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATg                                                                                                >  2:1591045/3‑150 (MQ=255)
                                                     cgtAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATg                                                                                                >  2:1616819/3‑150 (MQ=255)
                                                        aaTCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTg                                                                                             >  1:1815813/1‑150 (MQ=255)
                                                         aTCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTgg                                                                                            >  2:934431/1‑150 (MQ=255)
                                                            aGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTct                                                                                         >  2:1500211/1‑150 (MQ=255)
                                                                ttATCAACGGCGGGTTAAAGGGATTGATCCGCGATTTATGCAACCGCTGCTTAAACCTTTTCCACCGATGATGGACACGTTGCAAAATAATGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCTATTATCCTGATGCTGGTctggct                                                                                     >  2:1407846/1‑150 (MQ=255)
                                                                    cAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGtt                                                                                 >  1:813715/1‑150 (MQ=255)
                                                                      aCGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCg                                                                               <  1:1916893/150‑1 (MQ=255)
                                                                      aCGGCGGGTTAAAGCGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCCCCGAGGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCg                                                                               <  2:299343/150‑1 (MQ=255)
                                                                       cggcggGTTACAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGt                                                                              <  2:1060031/150‑1 (MQ=255)
                                                                        ggcggGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGtt                                                                             <  1:665349/150‑1 (MQ=255)
                                                                           gggTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAATATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATGCTGATGCTGGTCTGGCTGGTTCGGTAAc                                                                          >  1:1859577/1‑150 (MQ=255)
                                                                            ggTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAAcc                                                                         <  1:1025068/150‑1 (MQ=255)
                                                                             gTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTTCAAAATATGGCAACGCGTCCCTCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAAccc                                                                        <  2:907770/150‑1 (MQ=255)
                                                                               tAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCcg                                                                      >  2:103024/1‑150 (MQ=255)
                                                                                 aaGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCcgac                                                                    <  1:649772/150‑1 (MQ=255)
                                                                                          aTCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACa                                                           <  1:532180/150‑1 (MQ=255)
                                                                                            ccGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGt                                                         >  2:598942/1‑150 (MQ=255)
                                                                                              gcgTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTGGCTAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTAt                                                       <  2:1873988/150‑1 (MQ=255)
                                                                                                        cAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACaata                                             >  2:197144/1‑150 (MQ=255)
                                                                                                          aCCTCTGCTTAATCCATTTCCACCGAGGATTGACACTTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACaataat                                           <  1:159692/150‑1 (MQ=255)
                                                                                                           ccGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAataata                                          >  2:470283/1‑150 (MQ=255)
                                                                                                              ctgctTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATAccc                                       <  1:1035568/150‑1 (MQ=255)
                                                                                                                 ctTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGc                                    <  1:992555/150‑1 (MQ=255)
                                                                                                                  ttAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCCGCa                                   >  2:355020/1‑150 (MQ=255)
                                                                                                                        cTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATTCTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCag                             >  2:596912/1‑150 (MQ=255)
                                                                                                                           ttCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAgaga                          <  2:1132442/150‑1 (MQ=255)
                                                                                                                            tCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGAt                         >  1:1404365/1‑150 (MQ=255)
                                                                                                                            tCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGAt                         >  1:647637/1‑150 (MQ=255)
                                                                                                                             ccACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATc                        <  2:1241959/150‑1 (MQ=255)
                                                                                                                              cACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAATTGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCg                       >  1:1262663/1‑150 (MQ=255)
                                                                                                                              cACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCg                       >  1:1937359/1‑150 (MQ=255)
                                                                                                                              cACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCg                       >  2:34089/1‑150 (MQ=255)
                                                                                                                              cACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCg                       >  1:2017802/1‑150 (MQ=255)
                                                                                                                              cACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCACGCGCTGTGGCTACTGGTATTAGCGATTAACCTGATGCAGGTCTCGCTGGTACGTTAACCCCGACGAAAGGAAATTATCGTGAAAATAATACCCAGAACCGCAGAGATCg                       >  1:942616/1‑150 (MQ=255)
                                                                                                                                ccGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGcc                     <  2:676328/150‑1 (MQ=255)
                                                                                                                                 cGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGccc                    >  1:1993684/1‑150 (MQ=255)
                                                                                                                                 cGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGccc                    <  2:939844/150‑1 (MQ=255)
                                                                                                                                  gatgatGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGcccc                   >  2:2016785/1‑150 (MQ=255)
                                                                                                                                   ttgatgGAAACCTTGCAAAATATGGCAACGCGTCCCGCTCTGTGGATACTGTTAGTCGCGATTATCCTGATGATGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCg                  <  2:791753/149‑1 (MQ=255)
                                                                                                                                   gcgcgggACGCGTTGCCATATTTTGCAACGTGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCg                  >  2:398257/6‑150 (MQ=255)
                                                                                                                                     gcgggACGCGTTGCCATATTTTGCAACGTGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGc                >  2:1902899/4‑150 (MQ=255)
                                                                                                                                       gggACGCGTTGCCATATTTTGCAACGTGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCaa              >  1:1924731/2‑150 (MQ=255)
                                                                                                                                        ggACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAACGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGctat             >  2:342781/1‑147 (MQ=255)
                                                                                                                                         gACGCGTTGCCATATTTTGCAACGTGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGa            >  1:1166430/1‑150 (MQ=255)
                                                                                                                                         gACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGa            >  2:1538199/1‑150 (MQ=255)
                                                                                                                                         gACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGa            <  1:861730/150‑1 (MQ=255)
                                                                                                                                          acacGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGaa           <  1:1532734/150‑1 (MQ=255)
                                                                                                                                             cGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATa        >  2:974484/1‑150 (MQ=255)
                                                                                                                                              gTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATAc       >  1:1464811/1‑150 (MQ=255)
                                                                                                                                               ttGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATAcc      <  1:197144/150‑1 (MQ=255)
                                                                                                                                                tGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCg     <  2:1443232/150‑1 (MQ=255)
                                                                                                                                                 gCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGa    >  1:978436/1‑150 (MQ=255)
                                                                                                                                                 gCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCAGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGAAAGAAATACCGa    >  1:1649278/1‑150 (MQ=255)
                                                                                                                                                  cAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGaa   <  2:628384/150‑1 (MQ=255)
                                                                                                                                                  cAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGaa   >  2:1790092/1‑150 (MQ=255)
                                                                                                                                                   aaatTATGGCAACGCGTTCCGCGCTGTGTATACTTTTATTCGCGATTTTCCTGATGCTGGTCTGGCTGGTTCGTTAACCACGACGAAATTACAGTATCGTGACAATAATACACAGCACCGCAGAGAGCGCCCCGGCTAGAAATACCGAag  <  1:1407846/150‑1 (MQ=255)
                                                                                                                                                   aaaaTATGTCAACGCGTCACGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAag  <  2:1019298/150‑1 (MQ=255)
                                                                                                                                                   aaaaTATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAag  <  2:1703014/150‑1 (MQ=255)
                                                                                                                                                   aaaaTATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAag  <  1:1679776/150‑1 (MQ=255)
                                                                                                                                                   aaaaTATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAag  <  1:103024/150‑1 (MQ=255)
                                                                                                                                                   aaaaTATGGCAACGCGGCCCGCGCTGTGTATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGATATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAag  <  2:1711102/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAAG  >  NC_000913/2438722‑2439018

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