Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,186,367 G→T 15.5% Q124K (CAA→AAA)  rsd ← stationary phase protein, binds sigma 70 RNA polymerase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,186,3670GT15.5% 1050.8 / 111.9 320Q124K (CAA→AAA) rsdstationary phase protein, binds sigma 70 RNA polymerase subunit
Reads supporting (aligned to +/- strand):  ref base G (124/146);  new base T (28/21);  total (152/168)
Fisher's exact test for biased strand distribution p-value = 1.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 320 total aligned reads displayed.

ATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAATCTTCGCGGCGCGCGCGAGTTGCCCGTTGCCTTCCAGCTTATGAAGAATGCG  >  CP009273/4186219‑4186514
                                                                                                                                                    |                                                                                                                                                   
  ggggCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGaaa                                                                                                                                                  >  2:3604844/1‑150 (MQ=255)
  ggggCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGaaa                                                                                                                                                  >  2:1638818/1‑150 (MQ=255)
    ggCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTTAAACt                                                                                                                                                >  2:778864/1‑150 (MQ=255)
       aTTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCaa                                                                                                                                             >  1:3641895/1‑150 (MQ=255)
        ttGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATTTCAGATAAAACCTGTTGAAACTCAAg                                                                                                                                            <  1:2308794/150‑1 (MQ=255)
          gAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGt                                                                                                                                          <  1:397996/150‑1 (MQ=255)
           aaTGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTa                                                                                                                                         >  1:75227/1‑150 (MQ=255)
            aTGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTTAAACTCAAGGTaa                                                                                                                                        >  2:4308237/1‑150 (MQ=255)
             tGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAAt                                                                                                                                       <  2:2341406/150‑1 (MQ=255)
               tAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTa                                                                                                                                     <  1:3448912/150‑1 (MQ=255)
               tAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTa                                                                                                                                     <  2:2840489/150‑1 (MQ=255)
               tAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTa                                                                                                                                     <  1:2532387/150‑1 (MQ=255)
                aaaTTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTAt                                                                                                                                    <  1:3409589/150‑1 (MQ=255)
                 aaTTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATc                                                                                                                                   <  2:540094/150‑1 (MQ=255)
                   ttACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCAt                                                                                                                                 <  2:3404614/150‑1 (MQ=255)
                     aCGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGa                                                                                                                               >  2:1061593/1‑150 (MQ=255)
                      cgcgTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTTAAACTCAAGGTAATTATCATGAt                                                                                                                              >  2:1322504/1‑150 (MQ=255)
                        cgTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCg                                                                                                                            >  1:950225/1‑150 (MQ=255)
                          ttAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGAt                                                                                                                          <  2:3493711/150‑1 (MQ=255)
                           tAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTTAAACTCAAGGTAATTATCATGATCGATg                                                                                                                         <  1:3270886/150‑1 (MQ=255)
                           tAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATg                                                                                                                         >  1:1312829/1‑150 (MQ=255)
                            aaCAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATgg                                                                                                                        <  1:119433/150‑1 (MQ=255)
                             aCAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGc                                                                                                                       >  1:1003807/1‑150 (MQ=255)
                              cAGCGAAGAACTGAATCCGGATGTTTGTCGCGGGCGGCGTAAAGCACAAGCAGAATCAGCTTATCTTCCAGCACACAGCGCGCTTCCAGCGCCTCGCCGTTGTCAGATAATTCCTGTTGATACTCAAGGTAATTATCATGATCGAGGGCg                                                                                                                      <  2:2124162/150‑1 (MQ=255)
                               aGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCgg                                                                                                                     >  1:613263/1‑150 (MQ=255)
                                 cgcAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGtt                                                                                                                   <  2:2016069/150‑1 (MQ=255)
                                 cgcAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGtt                                                                                                                   >  2:402732/1‑150 (MQ=255)
                                  gcAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGttt                                                                                                                  >  2:1586019/1‑150 (MQ=255)
                                   cAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTc                                                                                                                 >  1:95886/1‑150 (MQ=255)
                                    aGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTcc                                                                                                                >  1:4260068/1‑150 (MQ=255)
                                    aGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTcc                                                                                                                >  2:4128380/1‑150 (MQ=255)
                                     gAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCa                                                                                                               >  1:3826715/1‑150 (MQ=255)
                                     gAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCa                                                                                                               <  1:3715055/150‑1 (MQ=255)
                                     gAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCa                                                                                                               <  2:1318311/150‑1 (MQ=255)
                                           aaGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTgg                                                                                                         <  2:3214706/150‑1 (MQ=255)
                                           aaGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTgg                                                                                                         >  1:1310122/1‑150 (MQ=255)
                                              cAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAAt                                                                                                      >  1:3883151/1‑150 (MQ=255)
                                               aGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATc                                                                                                     <  2:2995532/150‑1 (MQ=255)
                                               aGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATc                                                                                                     <  2:2493616/150‑1 (MQ=255)
                                               aGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATc                                                                                                     >  1:1289337/1‑150 (MQ=255)
                                               aGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATc                                                                                                     <  1:3072540/150‑1 (MQ=255)
                                               aGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATc                                                                                                     >  1:88120/1‑150 (MQ=255)
                                                ggATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCg                                                                                                    <  1:3340589/150‑1 (MQ=255)
                                                   tGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCgtag                                                                                                 >  1:813898/1‑150 (MQ=255)
                                                        gACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATc                                                                                            <  1:528842/150‑1 (MQ=255)
                                                         aCGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATcc                                                                                           <  1:1124113/150‑1 (MQ=255)
                                                         aCGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATcc                                                                                           <  2:1562139/150‑1 (MQ=255)
                                                            cgggcggCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATa                                                                                        >  1:683978/1‑150 (MQ=255)
                                                             gggcggCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATaa                                                                                       <  2:2663737/150‑1 (MQ=255)
                                                             gggcggCGTCAAGCACCAGCAGAATCAGCTTACCTTCCAGCACACAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTTAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATac                                                                                       <  1:4589445/150‑2 (MQ=255)
                                                                    gTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCANACTGGAATCGTAGTAATCCATAATCTGTTg                                                                                <  1:4285190/150‑1 (MQ=255)
                                                                    gTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTg                                                                                >  2:930024/1‑150 (MQ=255)
                                                                     tCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGa                                                                               <  1:367462/150‑1 (MQ=255)
                                                                       aaGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGt                                                                             >  1:195164/1‑150 (MQ=255)
                                                                        aGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTa                                                                            <  1:3967019/150‑1 (MQ=255)
                                                                        aGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTa                                                                            <  1:1581952/150‑1 (MQ=255)
                                                                           aCCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTg                                                                         <  2:1877942/150‑1 (MQ=255)
                                                                             cagcagAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGc                                                                       <  2:3688477/150‑1 (MQ=255)
                                                                                 agAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATTATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCg                                                                   >  1:2326172/1‑150 (MQ=255)
                                                                                   aaTCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAg                                                                 <  1:4214747/150‑1 (MQ=255)
                                                                                      cAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTg                                                              <  1:437793/150‑1 (MQ=255)
                                                                                       aGCTTATCTTCCAGCACAAAGCGCGCGTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGc                                                             <  1:4344254/150‑1 (MQ=255)
                                                                                         cTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTTAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCgg                                                           >  2:4250417/1‑150 (MQ=255)
                                                                                          ttATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGc                                                          >  2:2189183/1‑150 (MQ=255)
                                                                                           tATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGcc                                                         <  2:3708792/150‑1 (MQ=255)
                                                                                            aTCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCa                                                        <  2:1224674/150‑1 (MQ=255)
                                                                                            aTCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCa                                                        <  1:2831681/150‑1 (MQ=255)
                                                                                               ttCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAAt                                                     <  2:287289/150‑1 (MQ=255)
                                                                                               ttCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAAt                                                     >  1:3489218/1‑150 (MQ=255)
                                                                                                     cacaAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTTAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAATCTTcgc                                               <  2:3930833/150‑1 (MQ=255)
                                                                                                       caAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTTAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAATCTTCgcgg                                             <  2:2240451/150‑1 (MQ=255)
                                                                                                       caAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAATCTTCgcgg                                             >  2:3048167/1‑150 (MQ=255)
                                                                                                       caAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAATCTTCgcgg                                             >  2:2396336/1‑150 (MQ=255)
                                                                                                       caAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAATCTTCgcgg                                             >  1:2858967/1‑150 (MQ=255)
                                                                                                            cgcgcTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAATCTTCGCGgcgcgc                                        >  1:300892/1‑150 (MQ=255)
                                                                                                             gcgcTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAATCTTCGCGgcgcgcg                                       <  1:1721115/150‑1 (MQ=255)
                                                                                                             gcgcTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTAGTGGCTTCGAGTTGCGGCCAAATCTTCGCGgcgcgcg                                       >  2:2040360/1‑150 (MQ=255)
                                                                                                                cTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAATCTTCGCGGCGCGCGCGa                                    <  2:527591/150‑1 (MQ=255)
                                                                                                                     aGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGGTGCGGCCAAATCTTCGCGGCGCGCGCGAGTTGc                               >  2:1796695/1‑150 (MQ=255)
                                                                                                                       cgcCTCGCCGATGTCAGATAAAACCTGTTTAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAATCTTCGCGGCGCGCGCGAGTTGccc                             <  2:2894998/150‑1 (MQ=255)
                                                                                                                       cgcCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAATCTTCGCGGCGCGCGCGAGTTGccc                             <  2:2881952/150‑1 (MQ=255)
                                                                                                                       cgactctCCGATGTCAGATAAAAACTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATAGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGTGTTGCGGCCAAATCTTCGCGGCGCGCGCGAGTTGccc                             <  1:2228138/147‑1 (MQ=255)
                                                                                                                               cGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAATCTTCGCGGCGCGCGCGAGTTGCCCGTTGCCtt                     >  2:4541187/1‑150 (MQ=255)
                                                                                                                                gATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAATCTTCGCGGCGCGCGCGAGTTGCCCGTTGCCTTc                    >  2:629753/1‑150 (MQ=255)
                                                                                                                                gATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAATCTTCGCGGCGCGCGCGAGTTGCCCGTTGCCTTc                    >  1:558157/1‑150 (MQ=255)
                                                                                                                                gATGTCAGATAAAACCTGGTGAAAATCAAGTTAATTATCATGATCGATGGCGGATTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGATTCGAGTTGCGGCCAAATCTTCGCGGCGCGCGCGAGTTGCCCGTTGCCTTc                    >  1:1202074/1‑150 (MQ=255)
                                                                                                                                 aTGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAATCTTCGCGGCGCGCGCGAGTTGCCCGTTGCCTTcc                   >  1:1872366/1‑150 (MQ=255)
                                                                                                                                  tGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAATCTTCGCGGCGCGCGCGAGTTGCCCGTTGCCTTCCa                  <  1:2789915/150‑1 (MQ=255)
                                                                                                                                   gTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAATCTTCGCGGCGCGCGCGAGTTGCCCGTTGCCTTCCAg                 <  1:2758416/150‑1 (MQ=255)
                                                                                                                                   gTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAAGCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAATCTTCGCGGCGCGCGCGAGTTGCCCGTTGCCTTCCAg                 >  2:1432333/1‑150 (MQ=255)
                                                                                                                                     cAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAATCTTCGCGGCGCGCGCGAGTTGCCCGTTGCCTTCCAGCt               <  1:2875848/150‑1 (MQ=255)
                                                                                                                                      aGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTTATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAATCTTCGCGGCGCGCGCGAGTTGCCCGTTGCCTTCCAGCtt              <  2:3838274/150‑1 (MQ=255)
                                                                                                                                        aTAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAATCTTCGCGGCGCGCGCGAGTTGCCCGTTGCCTTCCAGCTTAt            >  1:4133439/1‑150 (MQ=255)
                                                                                                                                        aTAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAATCTTCGCGGCGCGCGCGAGTTGCCCGTTGCCTTCCAGCTTAt            <  1:1496112/150‑1 (MQ=255)
                                                                                                                                         tAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAATCTTCGCGGCGCGCGCGAGTTGCCCGTTGCCTTCCAGCTTATg           <  1:2407937/150‑1 (MQ=255)
                                                                                                                                            aaCCTGTTTAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAATCTTCGCGGCGCGCGCGAGTTGCCCGTTGCCTTCCAGCTTATgaag        >  1:1221936/1‑150 (MQ=255)
                                                                                                                                            aaCCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAATCTTCGCGGCGCGCGCGAGTTGCCCGTTGCCTTCCAGCTTATgaag        >  1:2616306/1‑150 (MQ=255)
                                                                                                                                                tGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAATCTTCGCGGCGCGCGCGAGTTGCCCGTTGCCTTCCAGCTTATGAAGAATg    >  1:1792999/1‑150 (MQ=255)
                                                                                                                                                 gTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAATCTTCGCGGCGCGCGCGAGTTGCCCGTTGCCTTCCAGCTTATGAAGAATgc   >  2:4395594/1‑150 (MQ=255)
                                                                                                                                                 gTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAATCTTCGCGGCGCGCGCGAGTTGCCCGTTGCCTTCCAGATTATGAAGAATgc   >  2:1696087/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATCTGTTGAGTATTGGCTTCGAGTTGCGGCCAAATCTTCGCGGCGCGCGCGAGTTGCCCGTTGCCTTCCAGCTTATGAAGAATGCG  >  CP009273/4186219‑4186514

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: