Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,438,874 A→T 5.3% M309L (ATG→TTG) ‡ flk → putative flagella assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,438,8740AT 5.3% 341.6 / 6.5 114M309L (ATG→TTG) ‡flkputative flagella assembly protein
Reads supporting (aligned to +/- strand):  ref base A (52/56);  new base T (4/2);  total (56/58)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

Only 100 of 115 total aligned reads displayed.

TGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAAGAG  >  NC_000913/2438726‑2439020
                                                                                                                                                    |                                                                                                                                                  
  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                                                                                                                                          <  2:1513017/150‑1 (MQ=255)
         tGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAAcgcg                                                                                                                                          >  1:1645902/1‑150 (MQ=255)
         tGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAAcgcg                                                                                                                                          >  1:1132442/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)
                     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                                                                                                                       <  2:436719/150‑1 (MQ=255)
                            tctGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTg                                                                                                                       >  1:1295850/1‑150 (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)
                                                tCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGAt                                                                                                   <  2:1504267/150‑1 (MQ=255)
                                                 cgtAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATg                                                                                                  >  2:1591045/3‑150 (MQ=255)
                                                    aaTCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTg                                                                                               >  1:1815813/1‑150 (MQ=255)
                                                     aTCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTgg                                                                                              >  2:934431/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)
                                                                     gcggGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGTTGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTa                                                                              <  2:1352923/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:647637/1‑150 (MQ=255)
                                                                                                                        tCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGAt                           >  1:1404365/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                         >  1:2017802/1‑150 (MQ=255)
                                                                                                                          cACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCg                         >  2:34089/1‑150 (MQ=255)
                                                                                                                          cACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCACGCGCTGTGGCTACTGGTATTAGCGATTAACCTGATGCAGGTCTCGCTGGTACGTTAACCCCGACGAAAGGAAATTATCGTGAAAATAATACCCAGAACCGCAGAGATCg                         >  1:942616/1‑150 (MQ=255)
                                                                                                                            ccGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGcc                       <  2:676328/150‑1 (MQ=255)
                                                                                                                             cGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGccc                      <  2:939844/150‑1 (MQ=255)
                                                                                                                             cGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGccc                      >  1:1993684/1‑150 (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)
                                                                                                                                                 aaTATGGCAACGCGGCCCGCGCTGTGGATACTGTTAGGCGCGTTTAGCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAagag  <  2:1251459/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAAGAG  >  NC_000913/2438726‑2439020

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: