Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,589,490 C→T 5.7% pseudogene (2755/5323 nt) yneO ← pseudogene, AidA homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,589,4900CT 5.7% 1000.5 / 29.5 279pseudogene (2755/5323 nt)yneOpseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base C (106/157);  new base T (7/9);  total (113/166)
Fisher's exact test for biased strand distribution p-value = 7.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 283 total aligned reads displayed.

TAACGTCAGCTTTTGCCCACCATGGATGTTCAGCGCTTTCAGCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCAT  >  CP009273/1589341‑1589639
                                                                                                                                                     |                                                                                                                                                     
tAACGTCAGCTTTTGCCCACCATGGATGTTCAGCGCTTTCAGCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGc                                                                                                                                                       >  2:1935314/1‑150 (MQ=255)
 aaCGTCAGCTTTTGCCCACCATGGATGTTCAGCGCTTTCAGCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCt                                                                                                                                                      >  1:2868764/1‑150 (MQ=255)
    gTCAGCTTTTGCCCACCATGGATGTTCAGCGCTTTCAGCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTcc                                                                                                                                                   <  2:2594658/150‑1 (MQ=255)
       aGCTTTTGCCCACCATGGATGTTCAGCGCTTTCAGCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCAcc                                                                                                                                                <  1:1885890/150‑1 (MQ=255)
         cTTTTGCCCACCATGGATGTTCAGCGCTTTCAGCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCtc                                                                                                                                              <  1:10698/150‑1 (MQ=255)
         cTTTTGCCCACCATGGATGTTCAGCGCTTTCAGCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCtc                                                                                                                                              <  2:1731510/150‑1 (MQ=255)
         cTTTTGCCCACCATGGATGTTCAGCGCTTTCAGCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCtc                                                                                                                                              <  1:238141/150‑1 (MQ=255)
            ttGCCCACCATGGATGTTCAGCGCTTTCAGCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCtcgtc                                                                                                                                           <  2:1629199/150‑1 (MQ=255)
            ttGCCCACCATGGATGTTCAGCGCTTTCAGCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCtcgtc                                                                                                                                           <  1:2477549/150‑1 (MQ=255)
              gCCCACCATGGATGTTCAGCGCTTTCAGCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGc                                                                                                                                         <  2:4689003/150‑1 (MQ=255)
               cccACCATGGATGTTCAGCGCTTTCAGCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGTTTCCACCTCGTCGcc                                                                                                                                        >  2:441023/1‑150 (MQ=255)
                ccaccaTGGATGTTCAGCGCTTTCAGCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGccc                                                                                                                                       <  2:635508/150‑1 (MQ=255)
                ccaccaTGGATGTTCAGCGCTTTCAGCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGccc                                                                                                                                       <  2:3006340/150‑1 (MQ=255)
                 caccaTGGATGTTCAGCGCTTTCAGCCCATAGTTAACGTACAAACCCTTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCt                                                                                                                                      <  2:860285/150‑1 (MQ=255)
                 caccaTGGATGTTCAGCGCTTTCAGCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGTTTCCACCTCGTCGCCCt                                                                                                                                      <  1:2513214/150‑1 (MQ=255)
                 caccaTGGATGTTCAGCGCTTTCAGCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCt                                                                                                                                      <  2:3040599/150‑1 (MQ=255)
                  accaTGGATGTTCAGCGCTTTCAGCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTg                                                                                                                                     <  2:820572/150‑1 (MQ=255)
                   ccaTGGATGTTCAGCGCTTTCAGCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGt                                                                                                                                    <  1:3216060/150‑1 (MQ=255)
                       ggATGTTCAGCGCTTTCAGCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAAt                                                                                                                                <  2:3222121/150‑1 (MQ=255)
                          tGTTCAGCGCTTTCAGCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCAc                                                                                                                             <  2:72042/150‑1 (MQ=255)
                                gcgcTTTCAGCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACg                                                                                                                       <  1:4610179/150‑1 (MQ=255)
                                gcgcTTTCAGCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACg                                                                                                                       <  2:3736072/150‑1 (MQ=255)
                                  gcTTTCAGCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGtt                                                                                                                     <  1:2182585/150‑1 (MQ=255)
                                   cTTTCAGCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGTTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTa                                                                                                                    <  1:355862/150‑1 (MQ=255)
                                       cAGCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGa                                                                                                                >  2:450235/1‑150 (MQ=255)
                                        aGCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAg                                                                                                               >  1:732769/1‑150 (MQ=255)
                                        aGCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCAGAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAg                                                                                                               >  2:4088409/1‑150 (MQ=255)
                                         gCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGc                                                                                                              >  2:3554052/1‑150 (MQ=255)
                                         gCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGc                                                                                                              >  1:2979508/1‑150 (MQ=255)
                                          cccATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGTTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGcc                                                                                                             >  2:4134329/1‑150 (MQ=255)
                                          cccATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGcc                                                                                                             >  2:1231615/1‑150 (MQ=255)
                                            cATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCAc                                                                                                           <  1:723268/150‑1 (MQ=255)
                                                gTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGa                                                                                                       >  1:169714/1‑150 (MQ=255)
                                                   aaCGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGAccc                                                                                                    >  1:780486/1‑150 (MQ=255)
                                                   aaCGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGAccc                                                                                                    >  1:1842222/1‑150 (MQ=255)
                                                      gTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAAt                                                                                                 <  2:2818431/150‑1 (MQ=255)
                                                       tACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATa                                                                                                <  1:2724138/150‑1 (MQ=255)
                                                         cAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAAc                                                                                              >  2:2002644/1‑150 (MQ=255)
                                                          aaaCCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGCAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAAcc                                                                                             >  2:1440352/1‑150 (MQ=255)
                                                              cATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGttt                                                                                         <  1:2052789/150‑1 (MQ=255)
                                                               aTTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTg                                                                                        <  2:2098735/150‑1 (MQ=255)
                                                                ttACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTgc                                                                                       <  2:1965213/150‑1 (MQ=255)
                                                                    cTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCt                                                                                   <  1:2545557/150‑1 (MQ=255)
                                                                     tGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCtt                                                                                  <  2:482313/150‑1 (MQ=255)
                                                                       gTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCa                                                                                <  2:328875/150‑1 (MQ=255)
                                                                       gTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCa                                                                                <  2:3349349/150‑1 (MQ=255)
                                                                        tGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGTTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCAc                                                                               >  2:15464/1‑150 (MQ=255)
                                                                           ggTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAg                                                                            >  2:1836162/1‑150 (MQ=255)
                                                                              cGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTg                                                                         <  2:765039/150‑1 (MQ=255)
                                                                                 gAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAAc                                                                      >  1:4193098/1‑150 (MQ=255)
                                                                                 gAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAAc                                                                      >  1:4029911/1‑150 (MQ=255)
                                                                                 gAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGGTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAAc                                                                      >  1:2878039/1‑150 (MQ=255)
                                                                                  aGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACg                                                                     <  2:714917/150‑1 (MQ=255)
                                                                                    gCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGcc                                                                   <  2:3634477/150‑1 (MQ=255)
                                                                                     cGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCg                                                                  <  1:1089813/150‑1 (MQ=255)
                                                                                      ggAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGa                                                                 >  2:2330510/1‑150 (MQ=255)
                                                                                        aaTCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACa                                                               <  1:4357665/150‑1 (MQ=255)
                                                                                        aaTCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACa                                                               <  1:1995955/150‑1 (MQ=255)
                                                                                         aTCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACAt                                                              <  2:2646736/150‑1 (MQ=255)
                                                                                          tCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATg                                                             <  2:150648/150‑1 (MQ=255)
                                                                                          tCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATg                                                             <  2:1341512/150‑1 (MQ=255)
                                                                                          tCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAAACGTTTGCGCCTTCACCAGTTGAACGCCGACATg                                                             >  1:3447838/1‑150 (MQ=255)
                                                                                            cATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGc                                                           >  1:2225377/1‑150 (MQ=255)
                                                                                            cATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGc                                                           <  2:2395115/150‑1 (MQ=255)
                                                                                            cATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGc                                                           >  2:4244139/1‑150 (MQ=255)
                                                                                              tAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTa                                                         >  1:3098498/1‑150 (MQ=255)
                                                                                              tAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTa                                                         <  2:3894715/150‑1 (MQ=255)
                                                                                              tAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAAGGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTa                                                         <  2:1404290/150‑1 (MQ=255)
                                                                                                   aCCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTGTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTAtcgtc                                                    >  1:4146660/1‑150 (MQ=255)
                                                                                                    cccTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTAtcgtcg                                                   <  1:4042429/150‑1 (MQ=255)
                                                                                                     ccTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATcgtcgt                                                  <  1:1330261/150‑1 (MQ=255)
                                                                                                           cGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGTTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCa                                            <  1:4436643/150‑1 (MQ=255)
                                                                                                           cGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCa                                            <  2:609320/150‑1 (MQ=255)
                                                                                                           cGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCa                                            >  2:1681548/1‑150 (MQ=255)
                                                                                                           cGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCa                                            <  2:1589381/150‑1 (MQ=255)
                                                                                                            gACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAg                                           >  1:4169023/1‑150 (MQ=255)
                                                                                                            gACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAg                                           <  2:2382485/150‑1 (MQ=255)
                                                                                                                 ccgttccgttTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAAt                                      >  1:3547098/1‑150 (MQ=255)
                                                                                                                  cgttccgttTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAAtt                                     >  1:3693597/1‑150 (MQ=255)
                                                                                                                   gttccgttTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTg                                    >  2:2572824/1‑150 (MQ=255)
                                                                                                                    ttccgttTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGa                                   >  2:4340704/1‑150 (MQ=255)
                                                                                                                     tccgttTTGCGCAATGTGCAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAg                                  <  1:1440352/150‑1 (MQ=255)
                                                                                                                       cgttTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGcg                                <  1:957399/150‑1 (MQ=255)
                                                                                                                       cgttTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGcg                                >  1:3359975/1‑150 (MQ=255)
                                                                                                                         ttttGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGcgtc                              <  1:604050/150‑1 (MQ=255)
                                                                                                                                 aaTGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAgg                      <  1:3906605/150‑1 (MQ=255)
                                                                                                                                  aTGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAggg                     <  2:3958283/150‑1 (MQ=255)
                                                                                                                                   tgtgTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGt                    <  2:935614/150‑1 (MQ=255)
                                                                                                                                    gtgtAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGtt                   <  2:271292/150‑1 (MQ=255)
                                                                                                                                    gtgtAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGtt                   <  2:272021/150‑1 (MQ=255)
                                                                                                                                     tgtAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTa                  <  1:3473425/150‑1 (MQ=255)
                                                                                                                                     tgtAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTa                  >  1:2340697/1‑150 (MQ=255)
                                                                                                                                       tAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTAtt                >  1:471994/1‑150 (MQ=255)
                                                                                                                                        aaCGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTc               <  2:4464231/150‑1 (MQ=255)
                                                                                                                                            ttttGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCa           >  2:1696534/1‑150 (MQ=255)
                                                                                                                                                  ccGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGAt     <  1:379377/150‑1 (MQ=255)
                                                                                                                                                  ccGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGAt     <  2:4071742/150‑1 (MQ=255)
                                                                                                                                                    gCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCa   <  2:2900526/150‑1 (MQ=255)
                                                                                                                                                    gCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCa   <  2:2900576/150‑1 (MQ=255)
                                                                                                                                                     cTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCAt  >  1:4222026/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TAACGTCAGCTTTTGCCCACCATGGATGTTCAGCGCTTTCAGCCCATAGTTAACGTACAAACCATTACCTGGTGCGGTCGTGAGGCGGAATCCATAATCACCCTCGGCGACCACCGTTCCGTTTTGCGCAATGTGTAACGTTTTGTCCGCTTCCACCTCGTCGCCCTGTAAATCACGTAACGTTAATGAGCCACCCGACCCAATAACCGTTTGCGCCTTCACCAGTTGAACGCCGACATGGCTATCGTCGTGTTCCAGCAAATTGAGCGTCGTCAGAGGGTTATTCGCCATGGGATCAT  >  CP009273/1589341‑1589639

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: