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.9% 674.6 / 4.6 242R267R (CGG→AGG) yhhIputative transposase
Reads supporting (aligned to +/- strand):  ref base C (1/4);  new base A (129/108);  total (130/112)
Fisher's exact test for biased strand distribution p-value = 1.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-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 242 total aligned reads displayed.

CCAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGT  >  CP009273/3618389‑3618681
                                                                                                                                                   |                                                                                                                                                 
 cAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGt                                                                                                                                                <  1:2486780/150‑1 (MQ=21)
  agagCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTc                                                                                                                                               >  2:723796/1‑150 (MQ=21)
   gagCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTcc                                                                                                                                              >  2:1509355/1‑150 (MQ=18)
   gagCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTcc                                                                                                                                              <  2:388730/150‑1 (MQ=18)
    agCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCa                                                                                                                                             >  1:764847/1‑150 (MQ=17)
    agCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCa                                                                                                                                             <  1:3787880/150‑1 (MQ=17)
      cATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCata                                                                                                                                           <  2:1118494/150‑1 (MQ=32)
      cATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCata                                                                                                                                           <  2:4740498/150‑1 (MQ=32)
           cAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGc                                                                                                                                      >  2:68136/1‑150 (MQ=17)
            aGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCa                                                                                                                                     <  1:2155671/150‑1 (MQ=32)
                   cAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACaaa                                                                                                                              <  2:4450811/150‑1 (MQ=17)
                   cAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACaaa                                                                                                                              >  1:1908186/1‑150 (MQ=17)
                   cAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACaaa                                                                                                                              >  1:1821156/1‑150 (MQ=17)
                      tCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAaaaga                                                                                                                           <  2:2952054/150‑1 (MQ=14)
                       cAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTTAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAaaagaa                                                                                                                          <  1:965205/150‑1 (MQ=11)
                       cAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAaaagaa                                                                                                                          <  2:1788172/150‑1 (MQ=11)
                             aaaGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGcc                                                                                                                    >  1:3547793/1‑150 (MQ=11)
                              aaGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCa                                                                                                                   >  1:3864836/1‑150 (MQ=14)
                                 agTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGaa                                                                                                                >  2:2343769/1‑150 (MQ=18)
                                   tCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAAt                                                                                                              >  1:2927911/1‑150 (MQ=21)
                                      cGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGAc                                                                                                           <  1:1250440/150‑1 (MQ=21)
                                      cGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGAc                                                                                                           <  1:4175624/150‑1 (MQ=21)
                                       ggCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACg                                                                                                          <  1:315623/150‑1 (MQ=21)
                                        gCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACgg                                                                                                         <  1:3713586/150‑1 (MQ=21)
                                        gCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACgg                                                                                                         >  2:3746132/1‑150 (MQ=21)
                                         cAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGt                                                                                                        <  1:3612695/150‑1 (MQ=21)
                                             gaagaaATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGa                                                                                                    >  2:4013084/1‑150 (MQ=21)
                                             gaagaaATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGa                                                                                                    >  1:1015123/1‑150 (MQ=21)
                                             gaagaaATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGGTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACCAAAGAAAGAGCCAGAAATGACGGTCAGa                                                                                                    >  1:2605808/1‑150 (MQ=255)
                                                 aaaTCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGAtatt                                                                                                >  1:4364005/1‑150 (MQ=21)
                                                 aaaTCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGAtatt                                                                                                >  2:1458752/1‑150 (MQ=21)
                                                  aaTCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGAtattc                                                                                               >  2:4629544/1‑149 (MQ=21)
                                                   aTCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGAtattat                                                                                              <  1:84105/150‑1 (MQ=21)
                                                     ccGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGGTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTatat                                                                                            >  2:5149061/1‑150 (MQ=21)
                                                         cTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGt                                                                                        <  1:2754292/150‑1 (MQ=21)
                                                           tCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTc                                                                                      >  1:2992140/1‑150 (MQ=21)
                                                            cATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTct                                                                                     >  1:1311813/1‑150 (MQ=21)
                                                                 tGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGa                                                                                <  2:2396419/150‑1 (MQ=21)
                                                                 tGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGa                                                                                >  2:59004/1‑150 (MQ=21)
                                                                  gTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGAt                                                                               <  1:2790692/150‑1 (MQ=21)
                                                                   tttGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGAtt                                                                              <  1:1018219/150‑1 (MQ=21)
                                                                   tttGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGAtt                                                                              <  1:3362464/150‑1 (MQ=21)
                                                                    ttGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGAttt                                                                             >  2:3071520/1‑150 (MQ=21)
                                                                     tGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGATAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTa                                                                            <  1:2067174/150‑1 (MQ=21)
                                                                     tGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTa                                                                            <  1:1731212/150‑1 (MQ=21)
                                                                     tGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTa                                                                            <  2:3203299/150‑1 (MQ=21)
                                                                     tGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTa                                                                            <  1:2466445/150‑1 (MQ=21)
                                                                      gCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTaa                                                                           >  1:1413705/1‑150 (MQ=21)
                                                                       cGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAAc                                                                          >  2:1999220/1‑150 (MQ=21)
                                                                            tCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTg                                                                     <  2:848140/150‑1 (MQ=25)
                                                                            tCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTg                                                                     >  1:441924/1‑150 (MQ=25)
                                                                             cccTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTga                                                                    <  1:2028054/150‑1 (MQ=37)
                                                                             cccTGATGAACTTATTGATTTCACGTTTGAAGGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAATAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTga                                                                    >  1:3298718/1‑150 (MQ=255)
                                                                              ccTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTgag                                                                   >  1:4259484/1‑150 (MQ=37)
                                                                                tgatgaACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGaa                                                                 <  2:2158989/150‑1 (MQ=37)
                                                                                tgatgaACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGaa                                                                 >  2:1935786/1‑150 (MQ=37)
                                                                                 gatgaACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGaaa                                                                >  1:487092/1‑150 (MQ=38)
                                                                                 gatgaACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGaaa                                                                >  2:2237334/1‑150 (MQ=38)
                                                                                 gatgaACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGaaa                                                                >  1:1112156/1‑150 (MQ=38)
                                                                                    gaACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTc                                                             >  1:3379390/1‑150 (MQ=255)
                                                                                     aaCTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCg                                                            <  2:3500593/150‑1 (MQ=255)
                                                                                      aCTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGc                                                           >  2:807210/1‑150 (MQ=255)
                                                                                       cTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGcc                                                          >  2:421311/1‑150 (MQ=255)
                                                                                        ttATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCca                                                         <  1:1277709/150‑1 (MQ=255)
                                                                                         tATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCTCcac                                                        >  2:4379620/1‑150 (MQ=255)
                                                                                           ttGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAg                                                      >  1:1043627/1‑150 (MQ=255)
                                                                                           ttGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAg                                                      >  1:4699838/1‑150 (MQ=255)
                                                                                               tttCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGAt                                                  <  1:5148399/150‑1 (MQ=255)
                                                                                                 tCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATcc                                                >  2:3793294/1‑150 (MQ=255)
                                                                                                  cACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCg                                               <  2:1808901/150‑1 (MQ=255)
                                                                                                    cGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGaa                                             >  2:455528/1‑150 (MQ=255)
                                                                                                    cGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGaa                                             <  1:394547/150‑1 (MQ=255)
                                                                                                     gTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGaaa                                            <  2:3452209/150‑1 (MQ=255)
                                                                                                     gTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGaaa                                            >  2:113871/1‑150 (MQ=255)
                                                                                                     gTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGaaa                                            <  1:3431441/150‑1 (MQ=255)
                                                                                                          aaTGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACt                                       >  2:4323288/1‑150 (MQ=255)
                                                                                                           aTGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTg                                      >  1:2340365/1‑150 (MQ=255)
                                                                                                            tGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTgg                                     <  2:1248356/150‑1 (MQ=255)
                                                                                                             ggATAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGc                                    <  2:2662784/150‑1 (MQ=255)
                                                                                                               aaaGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCAc                                  >  2:2891851/1‑150 (MQ=255)
                                                                                                                   ggCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTgg                              <  1:333486/150‑1 (MQ=255)
                                                                                                                      tGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGgaga                           >  1:3486571/1‑150 (MQ=255)
                                                                                                                       gaagaaATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGaa                          >  1:2517291/1‑150 (MQ=255)
                                                                                                                       gaagaaATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGaa                          >  2:3938439/1‑150 (MQ=255)
                                                                                                                        aagaaATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAAt                         >  1:4654767/1‑150 (MQ=255)
                                                                                                                           aaaTTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAg                      >  1:2686657/1‑150 (MQ=255)
                                                                                                                             aTTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCt                    >  1:4196333/1‑150 (MQ=255)
                                                                                                                               tATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGc                  >  1:2888801/1‑150 (MQ=255)
                                                                                                                               tATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGc                  >  1:1482770/1‑150 (MQ=255)
                                                                                                                                aTGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCa                 >  1:3508495/1‑150 (MQ=255)
                                                                                                                                 tGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCAc                <  2:2298025/150‑1 (MQ=255)
                                                                                                                                  gCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACt               <  1:1509355/150‑1 (MQ=255)
                                                                                                                                     tGGCAGTATCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGc            <  2:3842328/150‑1 (MQ=255)
                                                                                                                                      ggCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCg           <  2:603253/150‑1 (MQ=255)
                                                                                                                                          gTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTg       <  1:4122413/150‑1 (MQ=255)
                                                                                                                                          gTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTg       <  2:3332221/150‑1 (MQ=255)
                                                                                                                                            ctcCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGa     >  2:4974834/1‑150 (MQ=255)
                                                                                                                                              ccTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACg   <  2:1900331/150‑1 (MQ=255)
                                                                                                                                               cTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGt  <  2:352710/150‑1 (MQ=255)
                                                                                                                                               cTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGt  >  1:4462734/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGT  >  CP009273/3618389‑3618681

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: