Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,034,169 T→A 16.9% Q5L (CAA→CTA) 
A87A (GCA→GCT
cptA ←
sdhE ←
toxin of CptAB toxin‑antitoxin pair
flavinator of succinate dehydrogenase; antitoxin of CptAB toxin‑antitoxin pair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,034,1690TA16.9% 727.9 / 126.8 278Q5L (CAA→CTA) 
A87A (GCA→GCT
cptA
sdhE
toxin of CptAB toxin‑antitoxin pair
flavinator of succinate dehydrogenase; antitoxin of CptAB toxin‑antitoxin pair
Reads supporting (aligned to +/- strand):  ref base T (117/114);  new base A (26/21);  total (143/135)
Fisher's exact test for biased strand distribution p-value = 6.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 280 total aligned reads displayed.

ACGCAATCAAACACCACCAGCGAAAGTAACACCATCCATAACGGGGTGTAACTGAGTGGCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTC  >  CP009273/3034022‑3034315
                                                                                                                                                   |                                                                                                                                                  
 cgaaaacAAACACCACCAGCGAAATTAACACCATCCATACCGGGGTGTAACTGAGTGGCCCGGGCATGATTAAAATAACCGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTg                                                                                                                                                   <  1:3226335/142‑1 (MQ=255)
 cgaaaTCAAACACCACCAGCGAAAGTAACACCATCCATAACGGGGTGTAACTGAGGGGCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTg                                                                                                                                                   <  1:100092/145‑1 (MQ=255)
        aaaCACCACCAGCGAAAGTAACACCATCCATAACGGGGTGTAACTGAGTGGCCAGGGCATGAGTAAAATAACCGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAgg                                                                                                                                            <  2:2782678/148‑1 (MQ=255)
          acacCACCAGCGAAAGTAACACCATCCATAACGGGGTGTAACTGAGTGGCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGAc                                                                                                                                          >  1:313806/1‑148 (MQ=255)
          acacCACCAGCGAAAGTAACACCATCCATAACGGGGTGTAACTGAGTGGCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGAc                                                                                                                                          >  1:3207406/1‑148 (MQ=255)
               accaGCGAAAGTAACACCATCCATAACGGGGTGTAACTGAGTGGCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGt                                                                                                                                     <  2:3712782/148‑1 (MQ=255)
                  aGCGAAAGTAACACCATCCATAACGGGGTGTAACTGAGTGGCCAGGGCATGAGTAAAATCACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTcc                                                                                                                                  <  2:3253805/148‑1 (MQ=255)
                  aGCGAAAGTAACACCATCCATAACGGGGTGTAACTGAGTGGCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATAGCCACAGGACCACGTTcc                                                                                                                                  <  2:3776337/148‑1 (MQ=255)
                    cGAAAGTAACACCATCCATAACGGGGTGTAACTGAGGGGCCAGGGCATGAGGTAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGCCGCGCCAGGCAGATTGCCACAGGACCACGTTCCCg                                                                                                                                >  2:2358951/1‑148 (MQ=255)
                       aaGTAACACCATCCATAACGGGGTGTAACTGAGTGGCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGtt                                                                                                                             >  1:2341624/1‑148 (MQ=255)
                       aaGTAACACCATCCATAACGGGGTGTAACTGAGTGGCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGtt                                                                                                                             <  1:1893842/148‑1 (MQ=255)
                        aGTAACACCATCCATAACGGGGTGTAACTGAGTGGCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGtt                                                                                                                                    >  2:1286901/1‑140 (MQ=255)
                        aGTAACACCATCCATAACGGGGTGTAACTGAGTGGCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTc                                                                                                                            >  1:3742597/1‑148 (MQ=255)
                        aGTAACACCATCCATAACGGGGTGTAACTGAGTGGCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTc                                                                                                                            >  2:1684883/1‑148 (MQ=255)
                          tAACACCATCCATAACGGGGTGTAACTGAGTGGCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCg                                                                                                                          >  2:2231994/1‑148 (MQ=255)
                             cacCATCCATAACGGGGTGTAACTGAGTGGCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCgtgt                                                                                                                       >  1:2692535/1‑148 (MQ=255)
                               ccatccatAACGGGGTGTAACTGAGTGGCCAGGGCATGAGTAAAAGAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCt                                                                                                                     <  2:2000663/148‑1 (MQ=255)
                                 atccatAACGGGGTGTAACTGAGTGGCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTgg                                                                                                                   >  1:716783/1‑148 (MQ=255)
                                  tccatAACGGGGTGTAACTGAGTGGCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGa                                                                                                                  >  1:211616/1‑148 (MQ=255)
                                    catAACGGGGTGTAACTGAGTGGCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATg                                                                                                                <  1:2033935/148‑1 (MQ=255)
                                     atAACGGGGTGTAACTGAGTGGCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCTGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGa                                                                                                               >  1:2512350/1‑148 (MQ=255)
                                     atAACGGGGTGTAACTGAGTGGCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGa                                                                                                               >  2:1390096/1‑148 (MQ=255)
                                      tAACGGGGTGTAACTGAGTGGCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAg                                                                                                              <  2:2028379/148‑1 (MQ=255)
                                      tAACGGGGTGTAACTGAGTGGCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAg                                                                                                              >  2:3733216/1‑148 (MQ=255)
                                       aaCGGGGTGTAACTGAGTGGCCAGGGCATGAGTAAAATAACCGCGGCAACCCGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGt                                                                                                             <  1:1135952/148‑1 (MQ=255)
                                           gggTGTAACTGAGTGGCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATAGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGa                                                                                                         >  1:1770609/1‑148 (MQ=255)
                                             gtgtAACTGAGTGGCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATAGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGAcc                                                                                                       <  2:1344596/148‑1 (MQ=255)
                                              tgtAACTGAGTGGCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACca                                                                                                      >  2:154785/1‑148 (MQ=255)
                                              tgtAACTGAGTGGCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACca                                                                                                      >  2:852827/1‑148 (MQ=255)
                                               gtAACTGAGTGGCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATAGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACcat                                                                                                     <  2:2669681/148‑1 (MQ=255)
                                               gtAACTGAGTGGCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATAGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACcat                                                                                                     <  1:1549388/148‑1 (MQ=255)
                                                  aCTGAGTGGCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACcatcat                                                                                                  >  1:1229256/1‑148 (MQ=255)
                                                      aGTGGCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATAGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTcc                                                                                              >  1:835156/1‑148 (MQ=255)
                                                          gCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGtt                                                                                          <  2:1855963/148‑1 (MQ=255)
                                                              gggCATGAGTAAAATAACCGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGc                                                                                      <  2:1910964/148‑1 (MQ=255)
                                                              gggCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGc                                                                                      >  1:952679/1‑148 (MQ=255)
                                                               ggCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCa                                                                                     >  1:1737732/1‑148 (MQ=255)
                                                                gCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCAt                                                                                    <  2:2882317/148‑1 (MQ=255)
                                                                  aTGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCg                                                                                  >  2:641074/1‑148 (MQ=255)
                                                                   tGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATAGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCgg                                                                                 >  2:3741536/1‑148 (MQ=255)
                                                                    gAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGc                                                                                >  1:1489423/1‑148 (MQ=255)
                                                                    gAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGc                                                                                >  2:527619/1‑148 (MQ=255)
                                                                     aGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGcc                                                                               >  1:1710252/1‑148 (MQ=255)
                                                                       tAAAATAACCGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCgg                                                                             <  1:2184211/148‑1 (MQ=255)
                                                                        aaaaTAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGt                                                                            <  1:2010922/148‑1 (MQ=255)
                                                                          aaTAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGttt                                                                          >  1:283950/1‑148 (MQ=255)
                                                                               cAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGt                                                                     >  1:3484839/1‑148 (MQ=255)
                                                                                       aCCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATtcatc                                                             <  1:2425857/148‑1 (MQ=255)
                                                                                        ccAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATAGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATtcatca                                                            >  1:222014/1‑148 (MQ=255)
                                                                                            cccATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCa                                                        >  2:3342700/1‑148 (MQ=255)
                                                                                            cccATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCa                                                        >  1:75052/1‑148 (MQ=255)
                                                                                            cccATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCa                                                        >  1:232030/1‑148 (MQ=255)
                                                                                            cccATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGACACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCa                                                        >  2:1760961/1‑148 (MQ=255)
                                                                                            cccATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGACACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCa                                                        >  2:3496687/1‑148 (MQ=255)
                                                                                                tGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTGTACCGTGATTCATCAGCCAGTTa                                                    >  2:3071008/1‑148 (MQ=255)
                                                                                                 gAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTaa                                                   <  2:2411017/148‑1 (MQ=255)
                                                                                                 gAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTaa                                                   <  2:488647/148‑1 (MQ=255)
                                                                                                   aTCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAAc                                                 <  2:1516911/148‑1 (MQ=255)
                                                                                                   aTCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAAc                                                 <  1:2176851/148‑1 (MQ=255)
                                                                                                   aTCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAAc                                                 <  1:1804916/148‑1 (MQ=255)
                                                                                                   aTCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAAc                                                 <  2:2211538/148‑1 (MQ=255)
                                                                                                    tCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACa                                                <  2:47270/148‑1 (MQ=255)
                                                                                                    tCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACa                                                <  1:760000/148‑1 (MQ=255)
                                                                                                     cAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAg                                               >  1:1312061/1‑148 (MQ=255)
                                                                                                     cAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAg                                               <  1:2758320/148‑1 (MQ=255)
                                                                                                     cAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAg                                               >  1:2620386/1‑148 (MQ=255)
                                                                                                     cAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATAGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAg                                               <  2:835156/148‑1 (MQ=255)
                                                                                                       gCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGt                                             <  2:211616/148‑1 (MQ=255)
                                                                                                       gCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGt                                             >  2:3335341/1‑148 (MQ=255)
                                                                                                           ggAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCgg                                         >  2:3403616/1‑148 (MQ=255)
                                                                                                           ggAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCgg                                         <  1:3170297/148‑1 (MQ=255)
                                                                                                           ggAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCgg                                         >  2:3813469/1‑148 (MQ=255)
                                                                                                           ggAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATAGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCgg                                         <  1:1525365/148‑1 (MQ=255)
                                                                                                            gAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGa                                        >  1:2436253/1‑148 (MQ=255)
                                                                                                             aaaGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATAGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGAt                                       >  1:1050756/1‑148 (MQ=255)
                                                                                                                 ccACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTc                                   <  2:2519209/148‑1 (MQ=255)
                                                                                                                 ccACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGGTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTc                                   <  2:600491/148‑1 (MQ=255)
                                                                                                                       tgCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTc                             <  2:3524779/148‑1 (MQ=255)
                                                                                                                       tgCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTc                             <  1:273287/148‑1 (MQ=255)
                                                                                                                        gcgcgcCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTcc                            <  1:1390096/148‑1 (MQ=255)
                                                                                                                          gcgcCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCcag                          <  2:861477/148‑1 (MQ=255)
                                                                                                                                ggAGACGCGCAAATCAGATAGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCNAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACg                    >  1:1936465/1‑148 (MQ=255)
                                                                                                                                 gagaCGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGa                   >  1:1365527/1‑148 (MQ=255)
                                                                                                                                 gagaCGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGa                   >  2:2555679/1‑148 (MQ=255)
                                                                                                                                 gagaCGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGa                   >  2:2273624/1‑148 (MQ=255)
                                                                                                                                 gagaCGCGCAAATCAGATAGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGa                   >  1:35428/1‑148 (MQ=255)
                                                                                                                                  agaCGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGaa                  >  1:392806/1‑148 (MQ=255)
                                                                                                                                    aCGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATa                >  1:505362/1‑148 (MQ=255)
                                                                                                                                     cgcgCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATaa               <  2:3207406/148‑1 (MQ=255)
                                                                                                                                     cgcgCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATaa               <  2:313806/148‑1 (MQ=255)
                                                                                                                                       cgcAAATCAGATAGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAg             >  2:2635111/1‑148 (MQ=255)
                                                                                                                                        gcAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGa            <  1:3341906/148‑1 (MQ=255)
                                                                                                                                         cAAATCAGATAGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGAt           >  1:2745025/1‑148 (MQ=255)
                                                                                                                                            aTCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCg        <  1:1625416/148‑1 (MQ=255)
                                                                                                                                              cAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGtt      >  2:3527795/1‑148 (MQ=255)
                                                                                                                                              cAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGtt      <  2:665500/148‑1 (MQ=255)
                                                                                                                                              cAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGtt      <  2:2647519/148‑1 (MQ=255)
                                                                                                                                                gATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGtttt    <  2:2737812/148‑1 (MQ=255)
                                                                                                                                                gATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGtttt    >  2:1196739/1‑148 (MQ=255)
                                                                                                                                                gATAGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGtttt    <  2:453348/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ACGCAATCAAACACCACCAGCGAAAGTAACACCATCCATAACGGGGTGTAACTGAGTGGCCAGGGCATGAGTAAAATAACAGCGGCAACCAGCCCATGAATCAGCAAGGAAAGCCACTGTGCGCGCCAGGAGACGCGCAAATCAGATTGCCACAGGACCACGTTCCCGGTTCCGTGTCTGGATGAGTCGGACCATCATTTCCAGTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTC  >  CP009273/3034022‑3034315

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: