Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,860,327 C→A 58.4% intergenic (‑70/+42) ibpB ← / ← ibpA heat shock chaperone/heat shock chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,860,3270CA58.4% 73.3 / inf 286intergenic (‑70/+42)ibpB/ibpAheat shock chaperone/heat shock chaperone
Reads supporting (aligned to +/- strand):  ref base C (57/62);  new base A (86/81);  total (143/143)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 291 total aligned reads displayed.

CCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAACCATTTACCAGGTTAGCACCACGAACATGAATGTT  >  CP009273/3860178‑3860473
                                                                                                                                                     |                                                                                                                                                  
  ggCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGCGa                                                                                                                                                  >  1:1950152/1‑150 (MQ=255)
     gTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCa                                                                                                                                               >  2:232306/1‑150 (MQ=255)
         tGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGga                                                                                                                                           >  1:3688177/1‑150 (MQ=255)
             gTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATg                                                                                                                                       <  1:535190/150‑1 (MQ=255)
                cGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCa                                                                                                                                    >  1:3254480/1‑150 (MQ=255)
                 gTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAg                                                                                                                                   >  2:1602658/1‑150 (MQ=255)
                           ttGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAg                                                                                                                         <  2:3076560/150‑1 (MQ=255)
                                aaaaCCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAggcggc                                                                                                                    >  2:794298/1‑150 (MQ=255)
                                  aaCCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCt                                                                                                                  >  2:254760/1‑150 (MQ=255)
                                  aaCCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCt                                                                                                                  <  2:1831995/150‑1 (MQ=255)
                                   aCCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGACGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCtt                                                                                                                 <  1:1228063/150‑1 (MQ=255)
                                       aTCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAggg                                                                                                             <  2:1629628/150‑1 (MQ=255)
                                       aTCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAggg                                                                                                             <  2:3419572/150‑1 (MQ=255)
                                          cATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAAt                                                                                                          <  1:3016721/150‑1 (MQ=255)
                                             tGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAg                                                                                                       >  2:2626174/1‑150 (MQ=255)
                                             tGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAg                                                                                                       >  2:2077744/1‑150 (MQ=255)
                                              gACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGt                                                                                                      <  1:2777384/150‑1 (MQ=255)
                                                  cATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGAt                                                                                                  >  1:3958577/1‑150 (MQ=255)
                                                   aTCAGTGGGGATAAATCTAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTCACGGCGAGCCTGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGAtt                                                                                                 >  2:2933819/1‑150 (MQ=255)
                                                   aTCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGAtt                                                                                                 <  1:3094880/150‑1 (MQ=255)
                                                   aTCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGAtt                                                                                                 >  1:3893220/1‑150 (MQ=255)
                                                    tCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGAttt                                                                                                <  1:2594413/150‑1 (MQ=255)
                                                    tCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGAttt                                                                                                <  1:1025416/150‑1 (MQ=255)
                                                    tCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGAttt                                                                                                <  2:1360925/150‑1 (MQ=255)
                                                      aGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCg                                                                                              <  2:1161473/150‑1 (MQ=255)
                                                       gTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGACCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGa                                                                                             >  1:379803/1‑150 (MQ=255)
                                                       gTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGAAGTCAGGCCGCGCCAGGCGGCCTTAGGGCATTAGTTGATTTCGa                                                                                             >  2:2739445/1‑150 (MQ=255)
                                                         ggggATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATa                                                                                           <  1:466260/150‑1 (MQ=255)
                                                          gggATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATAc                                                                                          >  1:1440055/1‑150 (MQ=255)
                                                              tAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGgcg                                                                                      >  1:3632112/1‑150 (MQ=255)
                                                                 aTCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCgg                                                                                   >  2:4151078/1‑150 (MQ=255)
                                                                   cGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGtt                                                                                 <  2:565264/150‑1 (MQ=255)
                                                                      aGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGACTTAGTTGATTTCGATACGGCGCGGttttt                                                                              >  1:155467/1‑150 (MQ=255)
                                                                      aGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGttttt                                                                              <  1:1619340/150‑1 (MQ=255)
                                                                       gTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGtttttt                                                                             <  2:3280142/150‑1 (MQ=255)
                                                                        ttACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTc                                                                            <  2:2399796/150‑1 (MQ=255)
                                                                           cGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCCGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCt                                                                         >  1:3526392/1‑150 (MQ=255)
                                                                             cATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTc                                                                       >  2:3609432/1‑150 (MQ=255)
                                                                              aTAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTcc                                                                      <  2:4065303/150‑1 (MQ=255)
                                                                               tAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCATGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGTCGCGGTTTTTTCGCTTCCg                                                                     <  2:1907402/150‑1 (MQ=255)
                                                                                  tCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGaa                                                                  >  1:784499/1‑150 (MQ=255)
                                                                                   cATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAAt                                                                 >  2:119696/1‑150 (MQ=255)
                                                                                   cATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAAt                                                                 >  1:594637/1‑150 (MQ=255)
                                                                                    aTTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATc                                                                <  1:4090349/150‑1 (MQ=255)
                                                                                     tttCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCa                                                               >  1:522148/1‑150 (MQ=255)
                                                                                      ttCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCATCCGGAATCAc                                                              >  1:1793725/1‑150 (MQ=255)
                                                                                        ctcCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCAcgc                                                            >  2:3750201/1‑150 (MQ=255)
                                                                                         tcCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCAcgcg                                                           >  2:935153/1‑150 (MQ=255)
                                                                                          ccTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGt                                                          <  2:1554607/150‑1 (MQ=255)
                                                                                          ccTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGt                                                          >  1:3379544/1‑150 (MQ=255)
                                                                                            ttCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTc                                                        <  2:61572/150‑1 (MQ=255)
                                                                                             tCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCg                                                       <  1:1872532/150‑1 (MQ=255)
                                                                                             tCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCg                                                       <  1:2669622/150‑1 (MQ=255)
                                                                                               tAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGTCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCgag                                                     >  2:3572294/1‑150 (MQ=255)
                                                                                                aagaagCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCgaga                                                    <  2:3826087/150‑1 (MQ=255)
                                                                                                  gaagCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATc                                                  <  1:4059690/150‑1 (MQ=255)
                                                                                                   aagCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCg                                                 >  2:3536016/1‑150 (MQ=255)
                                                                                                   aagCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCAGGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCg                                                 >  2:3376706/1‑150 (MQ=255)
                                                                                                    agCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGa                                                >  1:1966807/1‑150 (MQ=255)
                                                                                                    agCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGa                                                >  2:452304/1‑150 (MQ=255)
                                                                                                     gCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGat                                               <  1:2161396/150‑1 (MQ=255)
                                                                                                     gCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGat                                               <  1:2950827/150‑1 (MQ=255)
                                                                                                      cGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGata                                              <  2:237094/150‑1 (MQ=255)
                                                                                                       gAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGatat                                             <  1:455237/150‑1 (MQ=255)
                                                                                                        agtaagtaCCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGAtata                                            <  2:52136/150‑1 (MQ=255)
                                                                                                         gtaagtaCCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGACTTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATAc                                           >  1:958503/1‑150 (MQ=255)
                                                                                                           aagtaCCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAg                                         >  2:3319547/1‑150 (MQ=255)
                                                                                                             gtaCCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCa                                       >  2:2189345/1‑150 (MQ=255)
                                                                                                              taCCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCaa                                      <  2:2250235/150‑1 (MQ=255)
                                                                                                               aCCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCaaa                                     <  2:2568231/150‑1 (MQ=255)
                                                                                                               aCCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCaaa                                     >  2:4109008/1‑150 (MQ=255)
                                                                                                                ccTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAAc                                    <  1:3244316/150‑1 (MQ=255)
                                                                                                                 cTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAAcc                                   >  1:410905/1‑150 (MQ=255)
                                                                                                                  tGCAAATCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAACCa                                  >  2:3094327/1‑150 (MQ=255)
                                                                                                                     aaaTCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAACCAttt                               >  1:3561978/1‑150 (MQ=255)
                                                                                                                     aaaTCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAACCAttt                               <  1:889289/150‑1 (MQ=255)
                                                                                                                     aaaTCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAACCAttt                               <  1:1646278/150‑1 (MQ=255)
                                                                                                                     aaaTCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAACCAttt                               >  1:3188016/1‑150 (MQ=255)
                                                                                                                     aaaTCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAACCAttt                               <  2:2725332/150‑1 (MQ=255)
                                                                                                                     aaaTCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAACCAttt                               <  1:3745031/150‑1 (MQ=255)
                                                                                                                      aaTCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAACCATTTa                              >  2:118441/1‑150 (MQ=255)
                                                                                                                       aTCCGAAGATTCGCATCTGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAACCATTTAc                             <  1:3845259/150‑1 (MQ=255)
                                                                                                                       aTCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAACCATTTAc                             <  1:1022911/150‑1 (MQ=255)
                                                                                                                        tCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAACCATTTAcc                            >  2:1214499/1‑150 (MQ=255)
                                                                                                                        tCCGAAGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAACCATTTAcc                            <  2:3887193/150‑1 (MQ=255)
                                                                                                                            aaGATTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAACCATTTACCAGGt                        >  2:3631897/1‑150 (MQ=255)
                                                                                                                               aTTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAACCATTTACCAGGTTAg                     <  2:2023600/150‑1 (MQ=255)
                                                                                                                               aTTCGCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAACCATTTACCAGGTTAg                     <  2:1696714/150‑1 (MQ=255)
                                                                                                                                 tCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAACCATTTACCAGGTTAGca                   <  2:3514963/150‑1 (MQ=255)
                                                                                                                                 tCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAACCATTTACCAGGTTAGca                   <  2:2767697/150‑1 (MQ=255)
                                                                                                                                   gCATATGCTCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAACCATTTACCAGGTTAGcacc                 <  2:2204963/150‑1 (MQ=255)
                                                                                                                                    cATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAACCATTTACCAGGTTAGcacca                >  2:1089129/1‑150 (MQ=255)
                                                                                                                                     atatGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAACCATTTACCAGGTTAGcaccac               >  1:4010982/1‑150 (MQ=255)
                                                                                                                                         gCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAACCATTTACCAGGTTAGCACCACGAAc           >  2:1327504/1‑150 (MQ=255)
                                                                                                                                           tCCCTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTAGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAACCATTTACCAGGTTAGCACCACGAACAt         <  1:2362549/150‑1 (MQ=255)
                                                                                                                                            cccTGACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAACCATTTACCAGGTTAGCACCACGAACATg        <  1:2494565/150‑1 (MQ=255)
                                                                                                                                             ccTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAACCATTTACCAGGTTAGCACCACGAACATGa       >  2:3222320/1‑150 (MQ=255)
                                                                                                                                               tGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAACCATTTACCAGGTTAGCACCACGAACATGAAt     <  2:4012659/150‑1 (MQ=255)
                                                                                                                                                gACGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAACCATTTACCAGGTTAGCACCACGAACATGAATg    <  1:2082354/150‑1 (MQ=255)
                                                                                                                                                  cGGAGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAACCATTTACCAGGTTAGCACCACGAACATGAATGtt  >  1:3938259/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CCGGCGTTTTGCAGTGCGTTGGCCAGTTTGTCAAAACCGATCCATTGACGCATCAGTGGGGATAAATCGAAGTTACGCATAGTCATTTCTCCTTCTAAGAAGCGAGTAAGTACCTGCAAATCCGAAGATTCGCATATGCTCCCTGACGGCGAGCATGGAGATGTCAGGCCGCGCCAGGCGGCCTTAGGGAATTAGTTGATTTCGATACGGCGCGGTTTTTTCGCTTCCGGAATCACGCGTTCGAGATCGATATACAGCAAACCATTTACCAGGTTAGCACCACGAACATGAATGTT  >  CP009273/3860178‑3860473

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: