Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,951,334 C→A 100% T352T (ACG→ACT torY ← TMAO reductase III (TorYZ), cytochrome c‑type subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,951,3340CA98.6% 749.2 / 0.8 233T352T (ACG→ACTtorYTMAO reductase III (TorYZ), cytochrome c‑type subunit
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base A (126/103);  total (128/105)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 233 total aligned reads displayed.

GTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCG  >  CP009273/1951189‑1951478
                                                                                                                                                 |                                                                                                                                                
gTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGt                                                                                                                                                <  2:2826293/148‑1 (MQ=255)
  gctgcCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTa                                                                                                                                              <  2:1631012/148‑1 (MQ=255)
    tgcCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTaaa                                                                                                                                            <  1:1148046/148‑1 (MQ=255)
    tgcCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTaaa                                                                                                                                            <  2:1667882/148‑1 (MQ=255)
        gATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGt                                                                                                                                        >  1:242177/1‑148 (MQ=255)
         aTGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTc                                                                                                                                       <  2:3342213/148‑1 (MQ=255)
          tGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTACTAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAAAGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTcc                                                                                                                                      <  2:1012489/148‑1 (MQ=255)
           gTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCa                                                                                                                                     >  1:374309/1‑148 (MQ=255)
            tCACCCCCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAg                                                                                                                                    <  2:2011926/148‑1 (MQ=255)
             caccacCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGt                                                                                                                                   <  1:3157257/148‑1 (MQ=255)
                caccaGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCg                                                                                                                                >  1:2446065/1‑148 (MQ=255)
                 accaGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGt                                                                                                                               <  2:3430677/148‑1 (MQ=255)
                 accaGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGt                                                                                                                               <  2:2356322/148‑1 (MQ=255)
                  ccaGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGtt                                                                                                                              >  2:1231877/1‑148 (MQ=255)
                  ccaGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGtt                                                                                                                              >  2:54626/1‑148 (MQ=255)
                   caGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGttt                                                                                                                             <  1:738969/148‑1 (MQ=255)
                    aGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGtttt                                                                                                                            <  2:2550681/148‑1 (MQ=255)
                    aGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGtttt                                                                                                                            <  2:3378162/148‑1 (MQ=255)
                        cTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAAAGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCt                                                                                                                        <  1:867913/148‑1 (MQ=255)
                        cTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAAAGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCt                                                                                                                        <  1:2648828/148‑1 (MQ=255)
                           cAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCAt                                                                                                                     >  1:2282948/1‑148 (MQ=255)
                            aGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATg                                                                                                                    <  2:685911/148‑1 (MQ=255)
                             gCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGc                                                                                                                   >  1:1486173/1‑148 (MQ=255)
                             gCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGc                                                                                                                   >  2:1958455/1‑148 (MQ=255)
                               ggCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTg                                                                                                                 <  2:1393606/148‑1 (MQ=255)
                                gCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTgg                                                                                                                >  2:3220351/1‑148 (MQ=255)
                                gCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTgg                                                                                                                <  2:85362/148‑1 (MQ=255)
                                 cAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGt                                                                                                               <  1:3314838/148‑1 (MQ=255)
                                 cAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGt                                                                                                               >  2:3518658/1‑148 (MQ=255)
                                 cAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGt                                                                                                               >  2:590177/1‑148 (MQ=255)
                                     cccGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGt                                                                                                           <  2:1975681/147‑1 (MQ=255)
                                      ccGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTg                                                                                                          <  2:1542610/148‑1 (MQ=255)
                                      ccGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTg                                                                                                          >  2:2597483/1‑148 (MQ=255)
                                       cGCTGTGTTTAATAAAGTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGc                                                                                                         <  2:2235366/148‑1 (MQ=255)
                                          tgtgTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGcccc                                                                                                      >  1:495362/1‑148 (MQ=255)
                                           gtgtTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGccccc                                                                                                     <  1:1612660/148‑1 (MQ=255)
                                              tttAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCAtt                                                                                                  >  1:2486696/1‑148 (MQ=255)
                                              tttAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCAtt                                                                                                  >  1:396271/1‑148 (MQ=255)
                                                taataaATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTcc                                                                                                >  2:2394448/1‑148 (MQ=255)
                                                taataaATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTcc                                                                                                >  2:1271249/1‑148 (MQ=255)
                                                 aataaATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCt                                                                                               <  1:1855931/148‑1 (MQ=255)
                                                  ataaATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCtt                                                                                              >  1:994429/1‑148 (MQ=255)
                                                    aaaTTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTg                                                                                            >  1:351738/1‑148 (MQ=255)
                                                       ttCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCTCTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCa                                                                                         <  2:2661216/148‑1 (MQ=255)
                                                       ttCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCa                                                                                         >  2:2486808/1‑148 (MQ=255)
                                                       ttCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCa                                                                                         >  2:3326641/1‑148 (MQ=255)
                                                         cACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAAt                                                                                       <  2:859106/148‑1 (MQ=255)
                                                           cGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATgg                                                                                     >  1:1333778/1‑148 (MQ=255)
                                                            gTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGa                                                                                    <  1:2573490/148‑1 (MQ=255)
                                                             tCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTNGAAATTCCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGAc                                                                                   <  1:1801082/148‑1 (MQ=255)
                                                               ttGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACgg                                                                                 <  1:2746185/148‑1 (MQ=255)
                                                                  ttAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCa                                                                              >  2:3280763/1‑148 (MQ=255)
                                                                     aTGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGc                                                                           >  1:3050275/1‑148 (MQ=255)
                                                                        tCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCAtt                                                                        >  1:2736608/1‑148 (MQ=255)
                                                                         cATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTg                                                                       <  1:939884/148‑1 (MQ=255)
                                                                         cATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTg                                                                       >  2:2810718/1‑148 (MQ=255)
                                                                         cATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTg                                                                       >  2:1873187/1‑148 (MQ=255)
                                                                            aaCTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACg                                                                    >  1:384324/1‑148 (MQ=255)
                                                                            aaCTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACg                                                                    >  2:2761173/1‑148 (MQ=255)
                                                                            aaCTTCCTCCCTGATCAACGAGGATCACTGTGTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACg                                                                    <  1:3098207/148‑1 (MQ=255)
                                                                              cTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGgt                                                                  >  1:2255757/1‑148 (MQ=255)
                                                                                    cccTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTa                                                               >  2:2860852/1‑145 (MQ=255)
                                                                                    cccTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTa                                                               <  1:2860852/145‑1 (MQ=255)
                                                                                      cTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGt                                                          >  2:1563533/1‑148 (MQ=255)
                                                                                      cTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGt                                                          <  1:1894903/148‑1 (MQ=255)
                                                                                      cTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGt                                                          <  2:1298347/148‑1 (MQ=255)
                                                                                        gATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTcg                                                        <  2:1206910/148‑1 (MQ=255)
                                                                                              cGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCa                                                  >  1:1439495/1‑148 (MQ=255)
                                                                                                aGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAAt                                                <  2:926228/148‑1 (MQ=255)
                                                                                                aGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAAt                                                >  2:1459174/1‑148 (MQ=255)
                                                                                                aGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAAt                                                <  1:2941063/148‑1 (MQ=255)
                                                                                                  gATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAg                                              >  1:1280526/1‑148 (MQ=255)
                                                                                                   aTCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAgg                                             >  2:2966556/1‑148 (MQ=255)
                                                                                                    tCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAggg                                            >  2:2724016/1‑148 (MQ=255)
                                                                                                      aCTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGc                                          <  2:2828076/148‑1 (MQ=255)
                                                                                                      aCTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCAGTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGc                                          <  2:1328720/148‑1 (MQ=255)
                                                                                                       cTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCa                                         >  1:454006/1‑148 (MQ=255)
                                                                                                        tGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCAt                                        >  1:1418680/1‑148 (MQ=255)
                                                                                                            tctcGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACa                                    <  2:374309/148‑1 (MQ=255)
                                                                                                              tcGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAAc                                  <  1:312088/148‑1 (MQ=255)
                                                                                                                   aaTATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAgcg                             >  1:2009391/1‑148 (MQ=255)
                                                                                                                   aaTATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAgcg                             >  1:2009002/1‑148 (MQ=255)
                                                                                                                   aaTATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAgcg                             >  1:2082638/1‑148 (MQ=255)
                                                                                                                   aaTATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAgcg                             >  1:2083437/1‑148 (MQ=255)
                                                                                                                    atatCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAgcgc                            >  1:822032/1‑148 (MQ=255)
                                                                                                                    atatCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAgcgc                            >  1:1052652/1‑148 (MQ=255)
                                                                                                                    atatCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAgcgc                            >  1:2831320/1‑148 (MQ=255)
                                                                                                                       tCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAAt                         >  1:261893/1‑148 (MQ=255)
                                                                                                                       tCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAAt                         <  2:1450450/148‑1 (MQ=255)
                                                                                                                        cTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATa                        <  1:870317/148‑1 (MQ=255)
                                                                                                                        cTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATa                        >  2:2708275/1‑148 (MQ=255)
                                                                                                                         tttGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAg                       <  2:659504/148‑1 (MQ=255)
                                                                                                                            ggCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTa                    <  2:1813314/148‑1 (MQ=255)
                                                                                                                               gTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCg                 <  1:948520/148‑1 (MQ=255)
                                                                                                                                  gTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCa              <  2:2097708/148‑1 (MQ=255)
                                                                                                                                    aCTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAg            <  2:3042282/148‑1 (MQ=255)
                                                                                                                                     cTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGt           <  2:2952424/148‑1 (MQ=255)
                                                                                                                                      tGGAAATACCGAGTTAAAATTTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGtt          >  1:2156883/1‑148 (MQ=255)
                                                                                                                                         aaaTACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGttttt       >  2:1638585/1‑148 (MQ=255)
                                                                                                                                            tACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAg    >  2:337661/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCG  >  CP009273/1951189‑1951478

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: