Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,951,127 C→A 36.6% 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,1270CA36.6% 369.7 / 234.2 251A46A (GCG→GCTtorZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base C (77/82);  new base A (48/44);  total (125/126)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 251 total aligned reads displayed.

AGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTT  >  CP009273/1950981‑1951269
                                                                                                                                                  |                                                                                                                                              
agaTAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCg                                                                                                                                               <  2:549980/148‑1 (MQ=255)
 gaTAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGc                                                                                                                                              <  1:1692870/148‑1 (MQ=255)
  aTAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCt                                                                                                                                             <  2:4386/148‑1 (MQ=255)
   tAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTc                                                                                                                                            <  1:2388941/148‑1 (MQ=255)
   tAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTc                                                                                                                                            >  2:2355125/1‑148 (MQ=255)
   tAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTc                                                                                                                                            >  1:950818/1‑148 (MQ=255)
    aGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTcc                                                                                                                                           >  1:3483975/1‑148 (MQ=255)
    aGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTcc                                                                                                                                           >  2:634451/1‑148 (MQ=255)
     gCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTccc                                                                                                                                          <  1:3394616/148‑1 (MQ=255)
     gCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTccc                                                                                                                                          <  1:1576026/148‑1 (MQ=255)
        ttttCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCaa                                                                                                                                       >  1:2786023/1‑148 (MQ=255)
        ttttCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCaa                                                                                                                                       >  1:1860837/1‑148 (MQ=255)
        ttttCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCaa                                                                                                                                       <  1:1683668/148‑1 (MQ=255)
            ctcACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTc                                                                                                                                   >  2:203586/1‑148 (MQ=255)
             tcACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTcc                                                                                                                                  <  1:1010658/148‑1 (MQ=255)
              cACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTccc                                                                                                                                 >  1:3050748/1‑148 (MQ=255)
               aCCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCcg                                                                                                                                <  2:816629/148‑1 (MQ=255)
                ccATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCcgc                                                                                                                               >  1:1206929/1‑148 (MQ=255)
                   tCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCATTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGt                                                                                                                            >  2:110228/1‑148 (MQ=255)
                    cGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTg                                                                                                                           >  1:2550681/1‑148 (MQ=255)
                     ggATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGa                                                                                                                          <  2:2612863/148‑1 (MQ=255)
                      gATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAg                                                                                                                         >  1:528170/1‑148 (MQ=255)
                      gATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAg                                                                                                                         >  1:2963249/1‑148 (MQ=255)
                       aTGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAgg                                                                                                                        >  1:2940611/1‑148 (MQ=255)
                         gCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGAt                                                                                                                      >  2:867913/1‑148 (MQ=255)
                           tGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGAttt                                                                                                                    <  1:392855/148‑1 (MQ=255)
                            gAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGAtttt                                                                                                                   >  1:2169905/1‑148 (MQ=255)
                             aaTACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTa                                                                                                                  <  1:2502464/148‑1 (MQ=255)
                              aTACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTAc                                                                                                                 >  1:2356322/1‑148 (MQ=255)
                                  gcgcCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGcc                                                                                                             <  1:675661/148‑1 (MQ=255)
                                     ccGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCctt                                                                                                          >  2:2941063/1‑148 (MQ=255)
                                      cGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCcttc                                                                                                         <  2:3306161/148‑1 (MQ=255)
                                      cGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCcttc                                                                                                         <  2:1800764/148‑1 (MQ=255)
                                        tggtgTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCcttctc                                                                                                       >  2:1445219/1‑148 (MQ=255)
                                         ggtgTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTGACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCcttctct                                                                                                      <  2:1706464/148‑1 (MQ=255)
                                         ggtgTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCcttctct                                                                                                      >  1:429090/1‑148 (MQ=255)
                                         ggtgTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCcttctct                                                                                                      >  1:2241211/1‑148 (MQ=255)
                                          gtgtgtACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCcttctctt                                                                                                     >  2:2041285/1‑148 (MQ=255)
                                           tgtgtACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCcttctcttc                                                                                                    <  1:1411229/148‑1 (MQ=255)
                                           tgtgtACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCcttctcttc                                                                                                    <  1:3178878/148‑1 (MQ=255)
                                             tgtACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTg                                                                                                  >  1:1220577/1‑148 (MQ=255)
                                                 ccTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCa                                                                                              <  2:2440214/148‑1 (MQ=255)
                                                     aTCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCg                                                                                          >  1:1975681/1‑148 (MQ=255)
                                                       ccgccgccgTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCggg                                                                                        <  1:2309613/148‑1 (MQ=255)
                                                           cgccgTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAAGCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCACCGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGc                                                                                    <  2:3071736/148‑1 (MQ=255)
                                                                tAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGtgc                                                                               >  2:2524050/1‑148 (MQ=255)
                                                                  gACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGtgctg                                                                             >  1:859106/1‑148 (MQ=255)
                                                                  gACTGTAAGGAATTCGGTATGGTTTTCGCCAGAGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGAGTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGtgctg                                                                             >  2:1862825/1‑148 (MQ=255)
                                                                    cTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTTTTCTGCCCATGCGGGGAGCGGTGCTGcc                                                                           >  2:738969/1‑148 (MQ=255)
                                                                      gTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGa                                                                         <  1:2632472/148‑1 (MQ=255)
                                                                        aaGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATg                                                                       <  2:1050726/148‑1 (MQ=255)
                                                                        aaGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATg                                                                       >  2:1971849/1‑148 (MQ=255)
                                                                         aGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGt                                                                      <  2:274454/148‑1 (MQ=255)
                                                                           gAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTca                                                                    <  1:2138702/148‑1 (MQ=255)
                                                                             aTTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTcacc                                                                  >  1:2747284/1‑148 (MQ=255)
                                                                             aTTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTcacc                                                                  >  1:1298347/1‑148 (MQ=255)
                                                                                    aTGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAgcg                                                           <  1:3189206/148‑1 (MQ=255)
                                                                                     tGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAgcgc                                                          <  2:1926358/148‑1 (MQ=255)
                                                                                       gTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTc                                                        <  1:534608/148‑1 (MQ=255)
                                                                                          ttCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAg                                                     >  1:1450450/1‑148 (MQ=255)
                                                                                             gCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCgg                                                  <  1:799981/148‑1 (MQ=255)
                                                                                               cAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCa                                                <  2:2141644/148‑1 (MQ=255)
                                                                                                  cgcgcCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATa                                             <  1:3441282/148‑1 (MQ=255)
                                                                                                  cgcgcCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATa                                             >  2:1604138/1‑148 (MQ=255)
                                                                                                  cgcgcCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATa                                             >  2:416364/1‑148 (MQ=255)
                                                                                                   gcgcCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATAc                                            <  2:2848790/148‑1 (MQ=255)
                                                                                                     gcCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCg                                          <  2:2169905/148‑1 (MQ=255)
                                                                                                        tGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCtg                                       >  1:2823768/1‑148 (MQ=255)
                                                                                                        tGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCtg                                       <  1:1968285/148‑1 (MQ=255)
                                                                                                        tGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCtg                                       <  2:247901/148‑1 (MQ=255)
                                                                                                        tGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCtg                                       <  2:1274163/148‑1 (MQ=255)
                                                                                                         gTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCtgt                                      >  2:411131/1‑148 (MQ=255)
                                                                                                         gTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCtgt                                      >  2:3401861/1‑148 (MQ=255)
                                                                                                            gaagaaACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGtt                                   >  1:2070586/1‑148 (MQ=255)
                                                                                                              agaaACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTta                                 <  2:3482223/148‑1 (MQ=255)
                                                                                                               gaaACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTtaa                                >  1:1493521/1‑148 (MQ=255)
                                                                                                               gaaACAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTtaa                                >  1:606520/1‑148 (MQ=255)
                                                                                                                aaaCAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTtaat                               <  1:165455/148‑1 (MQ=255)
                                                                                                                   cAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATaaa                            >  1:1328651/1‑148 (MQ=255)
                                                                                                                   cAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATaaa                            >  1:2828076/1‑148 (MQ=255)
                                                                                                                   cAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATaaa                            <  2:3204784/148‑1 (MQ=255)
                                                                                                                   cAATCTTGCCGTCCTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATaaa                            <  2:3529479/148‑1 (MQ=255)
                                                                                                                      tCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTc                         <  2:265616/148‑1 (MQ=255)
                                                                                                                       cTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCa                        <  2:185199/148‑1 (MQ=255)
                                                                                                                        ttGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCAc                       >  1:1044826/1‑148 (MQ=255)
                                                                                                                        ttGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCAc                       <  1:634960/148‑1 (MQ=255)
                                                                                                                           ccGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTc                    <  1:285280/148‑1 (MQ=255)
                                                                                                                           ccGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTc                    <  1:2079579/148‑1 (MQ=255)
                                                                                                                             gTCCTTCACTTCGACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCtt                  <  1:1828846/148‑1 (MQ=255)
                                                                                                                              tCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTg                 <  2:3456322/148‑1 (MQ=255)
                                                                                                                               ccTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGt                <  1:3226761/148‑1 (MQ=255)
                                                                                                                               ccTTCACTTCTACGTTCATAGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGt                >  1:976609/1‑148 (MQ=255)
                                                                                                                                cTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGtt               <  2:2872876/148‑1 (MQ=255)
                                                                                                                                  tCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTaa             >  1:3378162/1‑148 (MQ=255)
                                                                                                                                  tCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTaa             <  1:1343117/148‑1 (MQ=255)
                                                                                                                                    aCTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATg           <  1:2207426/148‑1 (MQ=255)
                                                                                                                                     cTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGt          <  2:3351340/148‑1 (MQ=255)
                                                                                                                                      ttCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTc         <  2:2786023/148‑1 (MQ=255)
                                                                                                                                         tACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATa      >  2:2573490/1‑148 (MQ=255)
                                                                                                                                             ttCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACtt  >  1:3342213/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
AGATAGCTTTTTCTCACCATCGGATGCTGAATACGCGCCGTGGTGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTT  >  CP009273/1950981‑1951269

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: