Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,438,876 G→T 6.5% M309I (ATG→ATT) ‡ flk → putative flagella assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,438,8760GT 6.5% 317.4 / 10.8 115M309I (ATG→ATT) ‡flkputative flagella assembly protein
Reads supporting (aligned to +/- strand):  ref base G (50/57);  new base T (5/2);  total (56/59)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-02

Only 100 of 116 total aligned reads displayed.

GTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAAGAGTAAC  >  NC_000913/2438728‑2439024
                                                                                                                                                    |                                                                                                                                                    
gTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATgg                                                                                                                                                     >  1:798417/1‑150 (MQ=255)
 tCTTTCTGCGCAATCATCGCTTTCTTCTGCCCCGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTGCACCGATGATGGACACGTTGCAAAATATGGc                                                                                                                                                    <  2:571088/150‑1 (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                                                                                                                                              >  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)
             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)
                   gCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGTTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCtgtg                                                                                                                                  <  1:142342/150‑1 (MQ=255)
                   gCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCtgtg                                                                                                                                  >  1:1019298/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)
                             gCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGCTAAAGGGAATCATCCGCGTGGTATGCAACCGCGGCTGAATCCTTTTCCACCGATGATGGACACGATGCAAAATATGGCAACGCGGCCCGCGCTGTGGATACTGTTa                                                                                                                        <  1:997002/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)
                                    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:1616819/3‑150 (MQ=255)
                                               cgtAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATg                                                                                                      >  2:1591045/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)
                                                           tATCAACGGCGGGTTAAAGGGAATGATCCGCGTGTTATGCAACCGCTGCTTAAACCTTTTCCACCGATGATGGACACGTTCCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTctggctg                                                                                          <  2:321049/150‑1 (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)
                                                                   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)
                                                                        ttAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAAcccc                                                                             <  1:178254/150‑1 (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)
                                                                                               aTGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCGGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACa                                                      <  1:816309/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)
                                                                                                           ctTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTATGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGc                                          <  1:588171/150‑1 (MQ=255)
                                                                                                            ttAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCCGCa                                         >  2:355020/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:2017802/1‑150 (MQ=255)
                                                                                                                        cACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCg                             <  1:1977720/150‑1 (MQ=255)
                                                                                                                        cACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCg                             >  1:1937359/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)
                                                                                                                                  ggACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAACGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGctat                   >  2:342781/1‑147 (MQ=255)
                                                                                                                                   gACGCGTTGCCATATTTTGCAACGTGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGa                  >  1:1166430/1‑150 (MQ=255)
                                                                                                                                   gACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGa                  <  1:861730/150‑1 (MQ=255)
                                                                                                                                   gACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGa                  >  2:1538199/1‑150 (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)
                                                                                                                                            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)
                                                                                                                                                   tGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTGCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAAGAGTAAc  >  1:1246065/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAAGAGTAAC  >  NC_000913/2438728‑2439024

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: