Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,062,424 G→T 84.6% I71I (ATC→ATA yihQ ← alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,062,4240GT84.6% 802.8 / 151.9 323I71I (ATC→ATAyihQalpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base G (17/33);  new base T (117/155);  total (135/188)
Fisher's exact test for biased strand distribution p-value = 2.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 324 total aligned reads displayed.

GCGCAGCCAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACGTTGTTGAAAATGTAG  >  CP009273/4062277‑4062573
                                                                                                                                                   |                                                                                                                                                     
  gcAGCCAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATTTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATgg                                                                                                                                                   >  1:3748786/1‑150 (MQ=255)
  gcAGCCAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATgg                                                                                                                                                   <  1:499793/150‑1 (MQ=255)
       cAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCg                                                                                                                                              <  2:2914635/150‑1 (MQ=255)
         gATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACTGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGt                                                                                                                                            >  2:5442178/1‑150 (MQ=255)
         gATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGt                                                                                                                                            >  2:4665872/1‑150 (MQ=255)
         gATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGt                                                                                                                                            >  1:5148753/1‑150 (MQ=255)
           tACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTaa                                                                                                                                          >  2:928008/1‑150 (MQ=255)
            aCGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAg                                                                                                                                         >  2:1010076/1‑150 (MQ=255)
              ggTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAgcg                                                                                                                                       <  1:488838/150‑1 (MQ=255)
               gTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAgcgc                                                                                                                                      >  1:614974/1‑150 (MQ=255)
                    ggTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAAttt                                                                                                                                 <  1:5462437/150‑1 (MQ=255)
                     gTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAAtttt                                                                                                                                <  1:4412905/150‑1 (MQ=255)
                       tAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTtct                                                                                                                              >  1:413983/1‑150 (MQ=255)
                         aGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTcc                                                                                                                            <  2:2109647/150‑1 (MQ=255)
                         aGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTcc                                                                                                                            >  1:3244453/1‑150 (MQ=255)
                          ggTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCt                                                                                                                           <  2:508259/150‑1 (MQ=255)
                           gTTGTCGTTTTGTTGTTCCTGCAATAAACGCCCCTGATCGTCGGCAGAGATTTTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTg                                                                                                                          <  2:4021933/150‑1 (MQ=255)
                           gTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTg                                                                                                                          <  1:1095309/150‑1 (MQ=255)
                           gTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGGAAGCGCAATTTTCTCCTg                                                                                                                          >  1:5404989/1‑150 (MQ=255)
                            ttGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGt                                                                                                                         >  1:2217936/1‑150 (MQ=255)
                             tGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCACAATTTTCTCCTGTa                                                                                                                        <  1:1211228/150‑1 (MQ=255)
                             tGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTa                                                                                                                        <  2:3200469/150‑1 (MQ=255)
                               tCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGt                                                                                                                      >  1:4799985/1‑150 (MQ=255)
                                cGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGtt                                                                                                                     <  1:3428892/150‑1 (MQ=255)
                                  ttttGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTa                                                                                                                   <  2:2192157/150‑1 (MQ=255)
                                  ttttGGAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTa                                                                                                                   <  2:373594/150‑1 (MQ=255)
                                   tttGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTAt                                                                                                                  >  1:5304107/1‑150 (MQ=255)
                                    ttGTAGTTCCAGCAATATACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATc                                                                                                                 <  2:4740388/150‑1 (MQ=255)
                                            ccAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGc                                                                                                         <  2:5544671/150‑1 (MQ=255)
                                             cAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGTCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCt                                                                                                        <  2:4874966/150‑1 (MQ=255)
                                             cAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCt                                                                                                        <  1:2729319/150‑1 (MQ=255)
                                                 aaTAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGaaa                                                                                                    <  2:222227/150‑1 (MQ=255)
                                                 aaTAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGaaa                                                                                                    >  1:1772321/1‑150 (MQ=255)
                                                   tAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAAtt                                                                                                  <  2:2238345/150‑1 (MQ=255)
                                                    aaaCGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTAAAATTa                                                                                                 <  2:1734560/150‑1 (MQ=255)
                                                        gCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACcgc                                                                                             >  1:4689934/1‑150 (MQ=255)
                                                        gCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACcgc                                                                                             >  2:3786087/1‑150 (MQ=255)
                                                           ccTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGa                                                                                          <  2:1546981/150‑1 (MQ=255)
                                                            cTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGaa                                                                                         <  2:832665/150‑1 (MQ=255)
                                                               aTCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACat                                                                                      >  1:2072135/1‑150 (MQ=255)
                                                                tcgtcgGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACata                                                                                     >  2:4513261/1‑150 (MQ=255)
                                                                 cgtcgGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACatat                                                                                    >  2:1822557/1‑150 (MQ=255)
                                                                  gtcgGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATc                                                                                   >  2:1069849/1‑150 (MQ=255)
                                                                  gtcgGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATc                                                                                   <  1:5196728/150‑1 (MQ=255)
                                                                   tcgGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCg                                                                                  <  1:45245/150‑1 (MQ=255)
                                                                       cAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGatat                                                                              <  1:1550902/150‑1 (MQ=255)
                                                                         gagaTATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATcc                                                                            <  2:3502886/150‑1 (MQ=255)
                                                                         gagaTATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATcc                                                                            >  1:4080012/1‑150 (MQ=255)
                                                                           gaTATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAATCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGc                                                                          >  1:1910288/1‑150 (MQ=255)
                                                                           gaTATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGc                                                                          >  2:680858/1‑150 (MQ=255)
                                                                           gaTATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGc                                                                          >  1:2949490/1‑150 (MQ=255)
                                                                            atatTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATACGCt                                                                         >  2:4844024/1‑150 (MQ=255)
                                                                              atTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCtat                                                                       <  2:3001905/150‑1 (MQ=255)
                                                                              atTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCtat                                                                       >  2:5005685/1‑150 (MQ=255)
                                                                               ttCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCtata                                                                      <  1:2474637/150‑1 (MQ=255)
                                                                               ttCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCtata                                                                      <  1:4238373/150‑1 (MQ=255)
                                                                                  aGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCt                                                                   <  2:3538820/150‑1 (MQ=255)
                                                                                   gCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTg                                                                  <  1:770582/150‑1 (MQ=255)
                                                                                   gCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTg                                                                  <  1:2038210/150‑1 (MQ=255)
                                                                                      tGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGc                                                               >  1:1048922/1‑150 (MQ=255)
                                                                                      tGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGc                                                               <  2:2021892/150‑1 (MQ=255)
                                                                                         cgcTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCaa                                                            <  1:738082/150‑1 (MQ=255)
                                                                                          gcTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAAt                                                           <  2:611749/150‑1 (MQ=255)
                                                                                             aaTGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCa                                                        >  2:113074/1‑150 (MQ=255)
                                                                                              aTGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCAt                                                       >  2:3838707/1‑150 (MQ=255)
                                                                                              aTGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCAt                                                       <  2:1995060/150‑1 (MQ=255)
                                                                                               tGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATa                                                      <  2:5657676/150‑1 (MQ=255)
                                                                                               tGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATa                                                      <  1:3926028/150‑1 (MQ=255)
                                                                                               tGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATa                                                      <  1:4039960/150‑1 (MQ=255)
                                                                                                   aGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACa                                                  >  2:1933853/1‑150 (MQ=255)
                                                                                                   aGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACa                                                  >  2:1625835/1‑150 (MQ=255)
                                                                                                      aCCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAgg                                               <  2:4235629/150‑1 (MQ=255)
                                                                                                      aCCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAgg                                               <  2:1302928/150‑1 (MQ=255)
                                                                                                      aCCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAgg                                               >  2:3112714/1‑150 (MQ=255)
                                                                                                             cTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCt                                        <  2:3640737/150‑1 (MQ=255)
                                                                                                                aaaTGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTg                                     <  2:3055082/150‑1 (MQ=255)
                                                                                                                    gAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTAt                                 <  2:1336934/150‑1 (MQ=255)
                                                                                                                    gAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTAt                                 <  2:2638965/150‑1 (MQ=255)
                                                                                                                    gAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTAt                                 <  2:5450487/150‑1 (MQ=255)
                                                                                                                     aaTTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATg                                >  1:2914559/1‑150 (MQ=255)
                                                                                                                     aaTTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATg                                >  1:5611850/1‑150 (MQ=255)
                                                                                                                     aaTTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATg                                <  2:2476269/150‑1 (MQ=255)
                                                                                                                         aaccaaccATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGtt                            <  1:949897/150‑1 (MQ=255)
                                                                                                                          accaaccaTCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTa                           <  2:3286700/150‑1 (MQ=255)
                                                                                                                           ccaaccaTCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTaa                          >  2:5311949/1‑150 (MQ=255)
                                                                                                                            caaccaTCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTaaa                         <  2:4287537/150‑1 (MQ=255)
                                                                                                                             aaccaTCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTaaaa                        <  1:2305665/150‑1 (MQ=255)
                                                                                                                               ccaTCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATa                      >  1:3840214/1‑150 (MQ=255)
                                                                                                                                caTCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATaa                     >  1:1931496/1‑150 (MQ=255)
                                                                                                                                   ccGGTGACTGGCTGACTATGGCGTCGGTAAGCTCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGATGAACATATCGATATCCGCTATACCTGATCCATTCCATAAACATGGATTATCTTTGCTATGGGTTAAAATAAGAc                  <  1:1140556/150‑1 (MQ=255)
                                                                                                                                   ccGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGAc                  <  1:3981378/150‑1 (MQ=255)
                                                                                                                                   ccGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGAc                  >  2:510299/1‑150 (MQ=255)
                                                                                                                                    cGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTAGACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACg                 <  1:3158504/150‑1 (MQ=255)
                                                                                                                                       tGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACgttg              <  1:3670339/150‑1 (MQ=255)
                                                                                                                                         aCTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACgttgtt            >  1:1326554/1‑150 (MQ=255)
                                                                                                                                         aCTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACgttgtt            >  1:3038127/1‑150 (MQ=255)
                                                                                                                                          cTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACGttgttg           <  2:1374308/150‑1 (MQ=255)
                                                                                                                                             gCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACGTTGTTGaaa        >  2:4355123/1‑150 (MQ=255)
                                                                                                                                             gCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACGTTGTTGaaa        <  2:25655/150‑1 (MQ=255)
                                                                                                                                              cTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACGTTGTTGaaaa       >  1:679945/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCGCAGCCAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACGTTGTTGAAAATGTAG  >  CP009273/4062277‑4062573

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: