Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,618,536 C→A 100% R267R (CGG→AGG)  yhhI → putative transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,618,5360CA97.8% 626.9 / 5.8 232R267R (CGG→AGG) yhhIputative transposase
Reads supporting (aligned to +/- strand):  ref base C (4/1);  new base A (114/113);  total (118/114)
Fisher's exact test for biased strand distribution p-value = 3.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 232 total aligned reads displayed.

GAGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACG  >  CP009273/3618392‑3618680
                                                                                                                                                |                                                                                                                                                
 agCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTc                                                                                                                                              >  2:1110701/1‑148 (MQ=17)
     tGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCata                                                                                                                                          >  1:3265205/1‑148 (MQ=17)
     tGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCata                                                                                                                                          >  1:580683/1‑148 (MQ=17)
      gACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCataa                                                                                                                                         <  2:1811205/148‑1 (MQ=17)
        cAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCataata                                                                                                                                       >  1:1006385/1‑148 (MQ=17)
        cAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCataata                                                                                                                                       <  2:760301/148‑1 (MQ=17)
           ttACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCa                                                                                                                                    <  2:515135/148‑1 (MQ=17)
              cGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGaa                                                                                                                                 <  1:771728/148‑1 (MQ=17)
                  aTCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACaaa                                                                                                                             >  1:639579/1‑148 (MQ=17)
                  aTCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACaaa                                                                                                                             <  2:1586678/148‑1 (MQ=17)
                   tCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACaaaa                                                                                                                            >  2:109395/1‑148 (MQ=14)
                    cAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAaaag                                                                                                                           >  1:1737950/1‑148 (MQ=11)
                            agagTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGcc                                                                                                                   <  2:2395971/148‑1 (MQ=11)
                                tCACGGCAGAGAAGAAATCCGTCTTCATATTGTATGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGaa                                                                                                               >  1:259828/1‑148 (MQ=18)
                                  aCGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAAt                                                                                                             <  2:2148613/148‑1 (MQ=21)
                                  aCGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAAt                                                                                                             <  1:2057690/148‑1 (MQ=21)
                                   cGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATg                                                                                                            <  1:878379/148‑1 (MQ=21)
                                   cGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATg                                                                                                            <  1:2596808/148‑1 (MQ=21)
                                    ggCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGa                                                                                                           <  1:225875/148‑1 (MQ=21)
                                      cAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACg                                                                                                         >  2:2019337/1‑148 (MQ=21)
                                      cAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACg                                                                                                         >  2:260092/1‑148 (MQ=21)
                                        gagaAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGt                                                                                                       <  1:38915/148‑1 (MQ=21)
                                         agaagaAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTc                                                                                                      >  2:2410452/1‑148 (MQ=21)
                                           aagaaATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAg                                                                                                    <  1:962723/148‑1 (MQ=21)
                                           aagaaATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAg                                                                                                    <  2:1284895/148‑1 (MQ=21)
                                           aagaaATCCGTCTTCATATGGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAg                                                                                                    <  2:3090836/148‑1 (MQ=21)
                                            agaaATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGa                                                                                                   <  1:236996/148‑1 (MQ=21)
                                             gaaATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGat                                                                                                  <  2:644375/148‑1 (MQ=21)
                                             gaaATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGat                                                                                                  >  2:3344460/1‑148 (MQ=21)
                                              aaaTCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGata                                                                                                 <  2:885967/148‑1 (MQ=21)
                                              aaaTCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGata                                                                                                 <  1:2737865/148‑1 (MQ=21)
                                              aaaTCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGata                                                                                                 <  2:2790040/148‑1 (MQ=21)
                                               aaTCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGatat                                                                                                <  1:3511825/148‑1 (MQ=21)
                                                aTCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGAtatt                                                                                               >  1:1543115/1‑148 (MQ=21)
                                                aTCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGAtatt                                                                                               <  1:28879/148‑1 (MQ=21)
                                                 tCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGAtatta                                                                                              >  2:1777308/1‑148 (MQ=21)
                                                  ccGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTTTGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGAtattat                                                                                             >  2:1892450/1‑148 (MQ=21)
                                                  ccGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGAtattat                                                                                             <  1:3000143/148‑1 (MQ=21)
                                                    gTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTatat                                                                                           <  2:2922507/148‑1 (MQ=21)
                                                    gTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTatat                                                                                           <  2:2547367/148‑1 (MQ=21)
                                                     tCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATc                                                                                          >  1:799743/1‑148 (MQ=21)
                                                     tCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATc                                                                                          <  2:501215/148‑1 (MQ=21)
                                                     tCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATc                                                                                          >  1:716884/1‑148 (MQ=21)
                                                      cTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAACAGAAAGAGCCAGAAATGACGGTCAGATATTATATCa                                                                                         >  2:2141013/1‑148 (MQ=21)
                                                      cTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCa                                                                                         <  2:629378/148‑1 (MQ=21)
                                                         cATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGtt                                                                                      <  1:670715/148‑1 (MQ=21)
                                                         cATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGtt                                                                                      <  1:1434204/148‑1 (MQ=21)
                                                          atatTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATTACGGTCAGATATTATATCAGTTc                                                                                     <  2:2715730/148‑1 (MQ=21)
                                                          atatTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTc                                                                                     <  1:3394344/148‑1 (MQ=21)
                                                            atTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAACTTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTAGATCAGTTctg                                                                                   >  2:1843041/1‑148 (MQ=255)
                                                            atTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTctg                                                                                   >  1:2172103/1‑148 (MQ=21)
                                                            atTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTctg                                                                                   <  1:1509299/148‑1 (MQ=21)
                                                              tGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTctgct                                                                                 >  1:43913/1‑148 (MQ=21)
                                                               gTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTctgctg                                                                                <  2:1227364/148‑1 (MQ=21)
                                                                tttGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGa                                                                               >  1:3030132/1‑148 (MQ=21)
                                                                tttGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGa                                                                               <  2:3259372/148‑1 (MQ=21)
                                                                    cGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTa                                                                           <  1:1000388/148‑1 (MQ=21)
                                                                      aTGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAAc                                                                         >  2:1347169/1‑148 (MQ=21)
                                                                      aTGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAAc                                                                         >  1:3503856/1‑148 (MQ=21)
                                                                       tGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAAcc                                                                        >  1:3486241/1‑148 (MQ=21)
                                                                         tCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGc                                                                      >  1:1127023/1‑148 (MQ=21)
                                                                         tCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGc                                                                      >  2:1749732/1‑148 (MQ=21)
                                                                           ccTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTg                                                                    >  1:1007508/1‑148 (MQ=25)
                                                                              gatgaACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTgaga                                                                 >  1:21522/1‑148 (MQ=37)
                                                                                 gaACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAAt                                                              >  2:321371/1‑148 (MQ=38)
                                                                                 gaACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAAt                                                              <  2:2030839/148‑1 (MQ=38)
                                                                                  aaCTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAAtt                                                             >  2:1364886/1‑148 (MQ=255)
                                                                                  aaCTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAAtt                                                             <  2:1560521/148‑1 (MQ=255)
                                                                                  aaCTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAAtt                                                             <  2:1560530/148‑1 (MQ=255)
                                                                                    cTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCg                                                           >  1:2758199/1‑148 (MQ=255)
                                                                                       aTTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCca                                                        >  2:1824999/1‑148 (MQ=255)
                                                                                       aTTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCca                                                        >  2:1257111/1‑148 (MQ=255)
                                                                                             ttCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGa                                                  >  1:1106447/1‑148 (MQ=255)
                                                                                             ttCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGa                                                  >  1:1653348/1‑148 (MQ=255)
                                                                                               cACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATc                                                <  2:1661458/148‑1 (MQ=255)
                                                                                               cACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATc                                                <  1:1832998/148‑1 (MQ=255)
                                                                                               cACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATc                                                <  2:2758736/148‑1 (MQ=255)
                                                                                                    ttGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGaaa                                           <  2:3455427/148‑1 (MQ=255)
                                                                                                    ttGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGaaa                                           >  1:2015642/1‑148 (MQ=255)
                                                                                                    ttGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGaaa                                           >  1:1798472/1‑148 (MQ=255)
                                                                                                       aaTGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCa                                        <  1:2019337/148‑1 (MQ=255)
                                                                                                        aTGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCAc                                       <  1:2651917/148‑1 (MQ=255)
                                                                                                         tGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACt                                      <  2:3042761/148‑1 (MQ=255)
                                                                                                          ggAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTg                                     >  2:1396460/1‑148 (MQ=255)
                                                                                                                ggCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGt                               >  2:456349/1‑148 (MQ=255)
                                                                                                                 gCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTg                              <  2:610378/148‑1 (MQ=255)
                                                                                                                    gaagaaATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGgag                           >  1:2448368/1‑148 (MQ=255)
                                                                                                                      agaaATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGaa                         >  1:2394259/1‑148 (MQ=255)
                                                                                                                       gaaATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAAt                        <  1:1850049/148‑1 (MQ=255)
                                                                                                                          aTTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAg                     <  2:3283401/148‑1 (MQ=255)
                                                                                                                           ttATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAATAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGc                    >  1:3072154/1‑148 (MQ=255)
                                                                                                                           ttATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATTTTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGc                    <  1:1299144/148‑1 (MQ=255)
                                                                                                                           ttATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGc                    <  1:2602148/148‑1 (MQ=255)
                                                                                                                                cGTGGCAGTCTCCTTTAGGTCCATAATAGCACAACAAAAGAAAGAGCAAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCAc               <  1:11563/148‑1 (MQ=255)
                                                                                                                                  tGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTg             <  2:3361583/148‑1 (MQ=255)
                                                                                                                                  tGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTg             <  1:3344460/148‑1 (MQ=255)
                                                                                                                                     cAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCg          >  2:227313/1‑148 (MQ=255)
                                                                                                                                        tctcCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCt       <  2:1403673/148‑1 (MQ=255)
                                                                                                                                        tctcCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCt       >  2:839338/1‑148 (MQ=255)
                                                                                                                                           ccTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGa    >  1:1542562/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
GAGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACG  >  CP009273/3618392‑3618680

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: ‑