Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,951,334 C→A 94.4% 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,3340CA94.4% 627.2 / 25.3 215T352T (ACG→ACTtorYTMAO reductase III (TorYZ), cytochrome c‑type subunit
Reads supporting (aligned to +/- strand):  ref base C (7/5);  new base A (106/97);  total (113/102)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 215 total aligned reads displayed.

GTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGC  >  CP009273/1951189‑1951477
                                                                                                                                                 |                                                                                                                                               
gTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGt                                                                                                                                               <  1:2281768/148‑1 (MQ=255)
gTGCTGCCGATGTCACCACCAGCGCTACAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGt                                                                                                                                               <  1:248612/148‑1 (MQ=255)
    tgcCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTaaa                                                                                                                                           <  1:380121/148‑1 (MQ=255)
    tgcCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTaaa                                                                                                                                           <  1:1845013/148‑1 (MQ=255)
     gcCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTaaaa                                                                                                                                          <  1:1152044/148‑1 (MQ=255)
     gcCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTaaaa                                                                                                                                          <  1:3106955/148‑1 (MQ=255)
     gcCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTaaaa                                                                                                                                          <  1:1649294/148‑1 (MQ=255)
       cGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATg                                                                                                                                        <  2:1703714/148‑1 (MQ=255)
        gATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGt                                                                                                                                       <  1:2196855/148‑1 (MQ=255)
           gTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCa                                                                                                                                    <  1:501779/148‑1 (MQ=255)
                 accaGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGt                                                                                                                              <  1:1710392/148‑1 (MQ=255)
                      cgcTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCg                                                                                                                         <  1:2939943/148‑1 (MQ=255)
                      cgcTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCg                                                                                                                         <  1:1396030/148‑1 (MQ=255)
                        cTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCt                                                                                                                       >  2:2266361/1‑148 (MQ=255)
                               ggCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTg                                                                                                                <  2:3471037/148‑1 (MQ=255)
                                 cAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGt                                                                                                              >  2:3388699/1‑148 (MQ=255)
                                 cAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGt                                                                                                              >  1:3300943/1‑148 (MQ=255)
                                 cAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGt                                                                                                              >  1:1423056/1‑148 (MQ=255)
                                   aTACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTc                                                                                                            >  2:2780032/1‑148 (MQ=255)
                                    tACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAACTACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCg                                                                                                           >  2:1541368/1‑148 (MQ=255)
                                    tACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCg                                                                                                           <  1:584141/148‑1 (MQ=255)
                                     aCCGCTGTGTTTAATAAATTCACGTCTTTTTAATGTCATATCTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCATTGGCGTTGTACTTGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGt                                                                                                          <  2:601923/148‑1 (MQ=255)
                                        gCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGcc                                                                                                       <  2:1774762/148‑1 (MQ=255)
                                         cTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCTCTCATGCTGGTTCGTGccc                                                                                                      >  2:3341159/1‑148 (MQ=255)
                                         cTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGccc                                                                                                      <  2:2684960/148‑1 (MQ=255)
                                           gtgtTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGccccc                                                                                                    >  2:2882850/1‑148 (MQ=255)
                                           gtgtTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGccccc                                                                                                    >  2:551271/1‑148 (MQ=255)
                                           gtgtTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGccccc                                                                                                    >  1:2121206/1‑148 (MQ=255)
                                               ttAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTc                                                                                                >  2:1498195/1‑148 (MQ=255)
                                                taataaATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTcc                                                                                               >  1:1336840/1‑148 (MQ=255)
                                                 aataaATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCt                                                                                              >  1:3241443/1‑148 (MQ=255)
                                                 aataaATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCt                                                                                              <  1:984996/148‑1 (MQ=255)
                                                   taaATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCttt                                                                                            <  2:2430055/148‑1 (MQ=255)
                                                   taaATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCttt                                                                                            <  2:2429883/148‑1 (MQ=255)
                                                      aTTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGc                                                                                         >  1:706959/1‑148 (MQ=255)
                                                       ttCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCa                                                                                        >  2:3007900/1‑148 (MQ=255)
                                                        tCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCaa                                                                                       <  2:832950/148‑1 (MQ=255)
                                                        tCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCaa                                                                                       >  1:1956905/1‑148 (MQ=255)
                                                          aCGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATg                                                                                     >  1:2388089/1‑148 (MQ=255)
                                                          aCGTCTTGTTAATGTAATAACTTCCTCCCAGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGTCGTTGTACTGGAAATACATAGTTAAAATGTCAAGTACGATTTCGCTCATGCTGGTTCTTGACACCATTCCATTGTCAATg                                                                                     >  2:913044/1‑148 (MQ=255)
                                                           cGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATgg                                                                                    >  1:106039/1‑148 (MQ=255)
                                                           cGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATg                                                                                     >  1:1218293/1‑147 (MQ=255)
                                                           cGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATg                                                                                     <  2:1218293/147‑1 (MQ=255)
                                                              cTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACg                                                                                 >  2:2313878/1‑148 (MQ=255)
                                                              cTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACg                                                                                 <  2:926175/148‑1 (MQ=255)
                                                                tGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGc                                                                               <  2:489506/148‑1 (MQ=255)
                                                                tGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGc                                                                               >  2:1914712/1‑148 (MQ=255)
                                                                tGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGc                                                                               <  2:2608447/148‑1 (MQ=255)
                                                                tGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGc                                                                               >  2:1586916/1‑148 (MQ=255)
                                                                      tGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCa                                                                         >  1:1796371/1‑148 (MQ=255)
                                                                        tCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCAtt                                                                       <  1:617046/148‑1 (MQ=255)
                                                                         cATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTg                                                                      >  2:1992619/1‑148 (MQ=255)
                                                                          aTAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGa                                                                     <  1:2744937/148‑1 (MQ=255)
                                                                          aTAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGa                                                                     >  1:2556594/1‑148 (MQ=255)
                                                                              cTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCCTTGACGgt                                                                 >  2:2371678/1‑148 (MQ=255)
                                                                               ttCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCAGTCCTTTGGCAATGGACGGCCACGCATTGACGgtg                                                                <  1:1861665/148‑1 (MQ=255)
                                                                                 cctccCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTa                                                              <  1:1417038/148‑1 (MQ=255)
                                                                                   tccCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAAt                                                            >  2:3286452/1‑148 (MQ=255)
                                                                                   tccCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAAt                                                            <  1:1192387/148‑1 (MQ=255)
                                                                                      cTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGt                                                         >  1:3039158/1‑148 (MQ=255)
                                                                                          tCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTcggc                                                     >  1:2395095/1‑148 (MQ=255)
                                                                                            aaCGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTcggcgg                                                   >  2:3266790/1‑148 (MQ=255)
                                                                                               gAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCaa                                                <  2:1423689/148‑1 (MQ=255)
                                                                                               gAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCaa                                                <  2:3300943/148‑1 (MQ=255)
                                                                                                aGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCact                                               >  1:1334848/1‑146 (MQ=255)
                                                                                                  gATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAg                                             <  2:2365584/148‑1 (MQ=255)
                                                                                                  gATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAg                                             <  1:729708/148‑1 (MQ=255)
                                                                                                      aCTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGc                                         >  1:1815119/1‑148 (MQ=255)
                                                                                                           ttCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGAc                                    >  2:674210/1‑148 (MQ=255)
                                                                                                             ctcGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTAGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACaa                                  <  1:1138694/148‑1 (MQ=255)
                                                                                                               cGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAAcc                                <  1:2665488/148‑1 (MQ=255)
                                                                                                               cGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAAcc                                <  2:1151018/148‑1 (MQ=255)
                                                                                                                 gTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAg                              >  1:728790/1‑148 (MQ=255)
                                                                                                                 gTAATATCTTTGGCGTTGTACTGGAAAGACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAg                              <  2:473792/148‑1 (MQ=255)
                                                                                                                  tAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAgc                             >  1:521474/1‑148 (MQ=255)
                                                                                                                  tAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAgc                             >  2:2225499/1‑148 (MQ=255)
                                                                                                                  tAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAgc                             >  2:431169/1‑148 (MQ=255)
                                                                                                                  tAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAgc                             >  2:3252843/1‑148 (MQ=255)
                                                                                                                   aaTATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAgcg                            >  1:2414289/1‑148 (MQ=255)
                                                                                                                    atatCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAgcgc                           >  2:2100736/1‑148 (MQ=255)
                                                                                                                     tatCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCa                          >  1:2262482/1‑148 (MQ=255)
                                                                                                                     tatCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCa                          <  1:90408/148‑1 (MQ=255)
                                                                                                                        cTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATa                       >  1:2886282/1‑148 (MQ=255)
                                                                                                                        cTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATa                       <  2:3427126/148‑1 (MQ=255)
                                                                                                                         tttGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAg                      >  2:3566156/1‑148 (MQ=255)
                                                                                                                             gCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTAt                  <  2:1608069/148‑1 (MQ=255)
                                                                                                                              cGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATc                 >  2:2637767/1‑148 (MQ=255)
                                                                                                                                ttGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGt               >  1:2067795/1‑148 (MQ=255)
                                                                                                                                 tGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTc              >  1:1470850/1‑148 (MQ=255)
                                                                                                                                     cTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGt          <  1:3348682/148‑1 (MQ=255)
                                                                                                                                     cTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGt          <  2:3203041/148‑1 (MQ=255)
                                                                                                                                       ggAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGttt        <  2:2666660/148‑1 (MQ=255)
                                                                                                                                         aaaTACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGttttt      >  2:1851893/1‑148 (MQ=255)
                                                                                                                                          aaTACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTc     >  2:2510721/1‑148 (MQ=255)
                                                                                                                                           aTACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCa    <  1:2925836/148‑1 (MQ=255)
                                                                                                                                           aTACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCa    <  1:2867221/148‑1 (MQ=255)
                                                                                                                                            tACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAg   <  2:1832604/148‑1 (MQ=255)
                                                                                                                                            tACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAg   <  2:1233278/148‑1 (MQ=255)
                                                                                                                                            tACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAg   <  2:1233227/148‑1 (MQ=255)
                                                                                                                                             aCCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGc  >  1:826130/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGC  >  CP009273/1951189‑1951477

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: