Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,951,127 C→A 20.5% A46A (GCG→GCT torZ ← trimethylamine N‑oxide reductase system III, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,951,1270CA20.5% 634.1 / 107.8 234A46A (GCG→GCTtorZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base C (93/93);  new base A (28/20);  total (121/113)
Fisher's exact test for biased strand distribution p-value = 3.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 236 total aligned reads displayed.

AGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCC  >  CP009273/1950981‑1951274
                                                                                                                                                  |                                                                                                                                                   
agaTAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCg                                                                                                                                                    >  2:380936/1‑148 (MQ=255)
  aTAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCt                                                                                                                                                  >  1:258132/1‑148 (MQ=255)
  aTAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCt                                                                                                                                                  >  1:259230/1‑148 (MQ=255)
  aTAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTAACTGATCCGCCGCCGTAGACTGTAAGGAATTAGGTATGGTTTTCTCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCt                                                                                                                                                  >  1:3303383/1‑148 (MQ=255)
   tAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTc                                                                                                                                                 >  1:1159001/1‑148 (MQ=255)
   tAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTc                                                                                                                                                 >  1:1159247/1‑148 (MQ=255)
   tAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTc                                                                                                                                                 >  2:1464818/1‑148 (MQ=255)
    aGCTTTTTCTCACCATCGGATGCTGAATACGCTCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTcc                                                                                                                                                >  1:1430741/1‑148 (MQ=255)
     gCTTTTTCTCACCATCGGATTCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTccc                                                                                                                                               <  2:2757255/148‑1 (MQ=255)
      cTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTcccc                                                                                                                                              >  1:3706097/1‑148 (MQ=255)
      cTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTcccc                                                                                                                                              <  2:287181/148‑1 (MQ=255)
            ctcACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTc                                                                                                                                        <  2:302385/148‑1 (MQ=255)
              cACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTccc                                                                                                                                      >  2:1483224/1‑148 (MQ=255)
               aCCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCcg                                                                                                                                     >  2:86062/1‑148 (MQ=255)
               aCCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCcg                                                                                                                                     <  1:2135509/148‑1 (MQ=255)
               aCCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCcg                                                                                                                                     <  1:206502/148‑1 (MQ=255)
                  aTCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCgg                                                                                                                                  <  2:336515/148‑1 (MQ=255)
                  aTCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCgg                                                                                                                                  >  1:920213/1‑148 (MQ=255)
                    cGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTg                                                                                                                                <  2:2708482/148‑1 (MQ=255)
                    cGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTg                                                                                                                                >  1:718219/1‑148 (MQ=255)
                    cGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCAGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTg                                                                                                                                <  2:545310/148‑1 (MQ=255)
                       aTGCTGCAGACGCGCCGTGGTGTGGACCTGATCCGCCGCCGTAGCCTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACACTCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAgg                                                                                                                             <  2:2223725/148‑1 (MQ=255)
                       aTGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAgg                                                                                                                             <  1:643162/148‑1 (MQ=255)
                        tGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCCTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGa                                                                                                                            <  1:2732965/148‑1 (MQ=255)
                         gCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGAt                                                                                                                           <  1:2104667/148‑1 (MQ=255)
                         gCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGAt                                                                                                                           >  1:3491089/1‑148 (MQ=255)
                         gCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGAt                                                                                                                           >  1:2482382/1‑148 (MQ=255)
                         gCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGAt                                                                                                                           <  1:1740488/148‑1 (MQ=255)
                            gAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGAtttt                                                                                                                        >  1:436581/1‑148 (MQ=255)
                             aaTACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTa                                                                                                                       <  1:3614671/148‑1 (MQ=255)
                             aaTACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTa                                                                                                                       >  1:2664592/1‑148 (MQ=255)
                                   cgcCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGCCCCGCGGTGAGGATTTTACCGccc                                                                                                                 >  2:197193/1‑148 (MQ=255)
                                     ccGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCctt                                                                                                               <  1:1898961/148‑1 (MQ=255)
                                      cGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCcttc                                                                                                              <  2:2112570/148‑1 (MQ=255)
                                                    gATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGc                                                                                                >  1:80440/1‑148 (MQ=255)
                                                      tCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCgg                                                                                              >  1:2340148/1‑148 (MQ=255)
                                                        cgccgccgTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCgggg                                                                                            >  1:1378490/1‑148 (MQ=255)
                                                             ccgTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCgg                                                                                       <  2:673129/148‑1 (MQ=255)
                                                               gTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGtg                                                                                     <  2:1523659/148‑1 (MQ=255)
                                                                tAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGtgc                                                                                    >  1:3211602/1‑148 (MQ=255)
                                                                tAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGtgc                                                                                    >  1:1960520/1‑148 (MQ=255)
                                                                   aCTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGtgctgc                                                                                 >  2:1683888/1‑148 (MQ=255)
                                                                   aCTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGtgctgc                                                                                 >  2:3578004/1‑148 (MQ=255)
                                                                   aCTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGACCTTCTCTTCTGCCCATGCGGGGAGCGGtgctgc                                                                                 >  2:1080099/1‑148 (MQ=255)
                                                                     tGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCg                                                                               <  2:2231471/148‑1 (MQ=255)
                                                                     tGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCg                                                                               >  1:3790347/1‑148 (MQ=255)
                                                                      gTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTCCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGa                                                                              >  1:3760469/1‑148 (MQ=255)
                                                                        aaGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATg                                                                            >  1:236839/1‑148 (MQ=255)
                                                                           gAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTca                                                                         >  1:2008704/1‑148 (MQ=255)
                                                                           gAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTca                                                                         >  1:3773372/1‑148 (MQ=255)
                                                                           gAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTca                                                                         >  2:656169/1‑148 (MQ=255)
                                                                            aaTTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTcac                                                                        >  2:2950768/1‑148 (MQ=255)
                                                                                    aTGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAgcg                                                                >  2:25596/1‑148 (MQ=255)
                                                                                    aTGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAgcg                                                                >  1:1642863/1‑148 (MQ=255)
                                                                                     tGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAgcgc                                                               >  1:1427704/1‑148 (MQ=255)
                                                                                     tGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAgcgc                                                               >  2:3745114/1‑148 (MQ=255)
                                                                                      ggTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCt                                                              >  1:3341916/1‑148 (MQ=255)
                                                                                             gCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCgg                                                       <  2:3666792/148‑1 (MQ=255)
                                                                                             gCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCgg                                                       <  2:1775029/148‑1 (MQ=255)
                                                                                               cAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCa                                                     <  2:2405721/148‑1 (MQ=255)
                                                                                                  cgcgcCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATa                                                  <  1:1483224/148‑1 (MQ=255)
                                                                                                     gcCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCg                                               <  1:1783010/148‑1 (MQ=255)
                                                                                                     gcCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCg                                               <  2:1378490/148‑1 (MQ=255)
                                                                                                      ccTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGc                                              >  2:2877409/1‑148 (MQ=255)
                                                                                                      ccTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGc                                              >  1:2276669/1‑148 (MQ=255)
                                                                                                           tGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTgtgt                                         >  2:2324737/1‑148 (MQ=255)
                                                                                                            gaagaaACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGtt                                        <  2:1003404/148‑1 (MQ=255)
                                                                                                            gaagaaACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGtt                                        <  1:1166915/148‑1 (MQ=255)
                                                                                                            gaagaaACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGtt                                        <  2:2340092/148‑1 (MQ=255)
                                                                                                            gaagaaACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGtt                                        <  2:2441472/148‑1 (MQ=255)
                                                                                                             aagaaACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGttt                                       <  1:526287/148‑1 (MQ=255)
                                                                                                                aaaCAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTtaat                                    >  1:914379/1‑148 (MQ=255)
                                                                                                                  aCAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTtaataa                                  <  2:3129335/148‑1 (MQ=255)
                                                                                                                  aCAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTtaataa                                  <  2:3211682/148‑1 (MQ=255)
                                                                                                                  aCAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTtaataa                                  <  1:1083185/148‑1 (MQ=255)
                                                                                                                   cAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATaaa                                 <  1:527260/148‑1 (MQ=255)
                                                                                                                   cAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATaaa                                 >  2:3339133/1‑148 (MQ=255)
                                                                                                                     aTCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAAtt                               >  2:3404291/1‑148 (MQ=255)
                                                                                                                      tCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTc                              <  2:1159548/148‑1 (MQ=255)
                                                                                                                      tCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTc                              >  1:3745359/1‑148 (MQ=255)
                                                                                                                      tCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTc                              >  1:3678363/1‑148 (MQ=255)
                                                                                                                      tCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTc                              >  1:1657838/1‑148 (MQ=255)
                                                                                                                      tCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCGTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTc                              <  2:1159094/148‑1 (MQ=255)
                                                                                                                       cTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCa                             <  2:1290532/148‑1 (MQ=255)
                                                                                                                       cTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCa                             <  2:911631/148‑1 (MQ=255)
                                                                                                                         tGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACg                           >  2:2781900/1‑148 (MQ=255)
                                                                                                                           ccGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTc                         >  1:1081127/1‑148 (MQ=255)
                                                                                                                            cGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCt                        <  1:479233/148‑1 (MQ=255)
                                                                                                                            cGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCt                        <  2:343778/148‑1 (MQ=255)
                                                                                                                              tCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTg                      <  1:3367793/148‑1 (MQ=255)
                                                                                                                              tCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTg                      >  1:3273342/1‑148 (MQ=255)
                                                                                                                              tCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTg                      <  1:1510993/148‑1 (MQ=255)
                                                                                                                              tCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTg                      <  1:1510779/148‑1 (MQ=255)
                                                                                                                               ccTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGt                     >  1:2247052/1‑148 (MQ=255)
                                                                                                                                  tCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTaa                  <  1:3578004/148‑1 (MQ=255)
                                                                                                                                          aCGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATaa          <  1:2950768/148‑1 (MQ=255)
                                                                                                                                          aCGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGAGCTCCAGCGGCAAAACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATaa          >  1:1005796/1‑148 (MQ=255)
                                                                                                                                           cGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAAc         >  1:1704164/1‑148 (MQ=255)
                                                                                                                                              tCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTtc      <  1:3032312/148‑1 (MQ=255)
                                                                                                                                                  cGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTtcctcc  <  1:380936/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCC  >  CP009273/1950981‑1951274

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