Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 161,171 G→T 8.7% intergenic (+150/‑46) hrpB → / → mrcB putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273161,1710GT 8.7% 695.0 / 33.7 190intergenic (+150/‑46)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base G (92/81);  new base T (7/9);  total (100/90)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 208 total aligned reads displayed.

GTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACG  >  CP009273/161023‑161319
                                                                                                                                                    |                                                                                                                                                    
ttgggttTTCCCCCTCGCCGGGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaag                                                                                                                                                      <  1:5174599/146‑1 (MQ=255)
gttgtttTTCGCCCTCGCGTGGTGGGGGGGGGCGGGGTTGTGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaag                                                                                                                                                      <  1:2989920/149‑1 (MQ=255)
gttgtttTTCCCCCTCTCCCTGTGGTTTTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaag                                                                                                                                                      <  1:2131144/149‑1 (MQ=255)
          tccctcgccggGTGGGTGGGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAAt                                                                                                                                            <  1:4317723/149‑1 (MQ=255)
             ctctccctGTGGGAGAGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                                                                                                                         <  2:4267864/149‑1 (MQ=255)
                tccgtGGGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCt                                                                                                                                      <  2:3862170/149‑1 (MQ=255)
                 ccctGTGGGAGAGGGCCGGGGGGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtt                                                                                                                                     <  2:1762752/149‑1 (MQ=255)
                     gtgGTAGAGGGCCGGGGTTAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAATAACAGAAAATCGGGCTTTTgc                                                                                                                                 <  2:571469/149‑1 (MQ=255)
                        ggAGGGGCCCGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGcc                                                                                                                              <  2:4058820/149‑1 (MQ=255)
                            ggggCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGaa                                                                                                                          <  1:3335410/148‑1 (MQ=255)
                                   gggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                                                                                                   <  2:2918187/149‑1 (MQ=255)
                                          gccATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAATAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa                                                                                                            <  1:3386787/147‑1 (MQ=255)
                                              tccgcgcgcACGATCCCCCTCACCCGAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTCGTAACCGAAGAACAGAAAATCGGGCTTTTGCGCCGGAATATTGCGGAGAAAAAGCa                                                                                                        <  2:4171369/146‑1 (MQ=255)
                                              tCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCa                                                                                                        <  2:3625431/149‑1 (MQ=255)
                                               cAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCAt                                                                                                       <  2:4738227/149‑1 (MQ=255)
                                                  cgcgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGc                                                                                                    <  2:455660/149‑1 (MQ=255)
                                                    cgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCg                                                                                                  <  1:1379235/149‑1 (MQ=255)
                                                         gTTCACCCTCACCCTAAACCACACCCTCAAGGGAGAGGGGACCGATCGAGCACAAAATTTGAGAGATAAATTCATCTGTCTTGTAACAGAAGAAAAGAAAAACGGGCTTTTGCGCCAGAATATTGCGGAGAAAAAGCAAGGCCGGGaaa                                                                                             >  2:2286191/1‑148 (MQ=255)
                                                          ttCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATg                                                                                            <  2:682737/149‑1 (MQ=255)
                                                           tcaccctcacccTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGa                                                                                           <  2:4765445/149‑1 (MQ=255)
                                                             accctcaccctAACCCGCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGAcc                                                                                         <  2:3129358/149‑1 (MQ=255)
                                                              ccctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACcg                                                                                        >  2:2326896/1‑149 (MQ=255)
                                                              ccctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACcg                                                                                        >  1:1256874/1‑149 (MQ=255)
                                                               cctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACcgc                                                                                       >  1:4685199/1‑149 (MQ=255)
                                                                ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACcgcg                                                                                      >  1:4468898/1‑149 (MQ=255)
                                                                ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACcgcg                                                                                      >  1:1818540/1‑149 (MQ=255)
                                                                   accctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAATAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGc                                                                                   >  1:4578718/1‑149 (MQ=255)
                                                                   accctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGc                                                                                   >  2:4989811/1‑149 (MQ=255)
                                                                    ccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAATAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGcc                                                                                  >  2:3043198/1‑149 (MQ=255)
                                                                    ccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGcc                                                                                  <  2:718148/149‑1 (MQ=255)
                                                                     cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCa                                                                                 >  1:1553838/1‑149 (MQ=255)
                                                                      ctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa                                                                                <  1:1729428/149‑1 (MQ=255)
                                                                        aaCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt                                                                              >  2:334640/1‑149 (MQ=255)
                                                                          cccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTgg                                                                            >  1:1946284/1‑149 (MQ=255)
                                                                           ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa                                                                           <  1:2578198/149‑1 (MQ=255)
                                                                                cccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAATAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCaa                                                                      >  1:1866767/1‑149 (MQ=255)
                                                                                cccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCaa                                                                      >  2:3787146/1‑149 (MQ=255)
                                                                                  cTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAg                                                                    >  2:3395767/1‑149 (MQ=255)
                                                                                  cTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAg                                                                    <  1:1102069/149‑1 (MQ=255)
                                                                                      aGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGaa                                                                >  2:2508535/1‑149 (MQ=255)
                                                                                       gggAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGaaa                                                               <  2:1553838/149‑1 (MQ=255)
                                                                                        ggAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAATAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAAc                                                              <  2:4857473/149‑1 (MQ=255)
                                                                                        ggAGAGGGGACAGATCGAGCACAAATTTTGAGAGATCTCTTCGTCTGTCTTGGACCAGAAGAACAGAAAAGCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAAGGACCGCGAGCCAATTGGACGCAAAGGGAAAc                                                              <  2:3509445/149‑1 (MQ=255)
                                                                                             ggggACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGAcg                                                         >  2:2020935/1‑149 (MQ=255)
                                                                                              gggACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGACTTGTAACAGAAGAACATAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAAGGACCGCGAGCCAATTGGACGCAAAGGGAAACCGAcgc                                                        <  1:487464/149‑1 (MQ=255)
                                                                                               ggACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGAcgcg                                                       >  2:1055911/1‑149 (MQ=255)
                                                                                                gACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGt                                                      <  2:3698532/149‑1 (MQ=255)
                                                                                                 aCCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTc                                                     >  1:3949692/1‑149 (MQ=255)
                                                                                                 aCCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTc                                                     <  1:3308788/149‑1 (MQ=255)
                                                                                                 aCCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTc                                                     <  2:1049170/149‑1 (MQ=255)
                                                                                                  ccGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTcc                                                    >  1:1777169/1‑149 (MQ=255)
                                                                                                    gATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGACTGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCgg                                                  >  1:1427498/1‑149 (MQ=255)
                                                                                                         aGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTcaaa                                             <  2:1655916/149‑1 (MQ=255)
                                                                                                          gCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTcaaac                                            >  1:1076762/1‑149 (MQ=255)
                                                                                                          gCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTcaaac                                            >  1:1076515/1‑149 (MQ=255)
                                                                                                             caAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTcaaacaaa                                         >  1:3189960/1‑149 (MQ=255)
                                                                                                               aaTTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAg                                       >  2:2851885/1‑149 (MQ=255)
                                                                                                                aTTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAgg                                      >  1:2400030/1‑149 (MQ=255)
                                                                                                                  tttGAGAGATATCTTCTTCTGTCTTGTAACAGAATAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTa                                    <  1:1100087/149‑1 (MQ=255)
                                                                                                                   ttGAGAGATATCTTCTTCTGTCTTGTAACAGAATAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTaa                                   <  1:787224/149‑1 (MQ=255)
                                                                                                                    tGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAg                                  <  1:5242566/149‑1 (MQ=255)
                                                                                                                    tGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAg                                  >  1:3527423/1‑149 (MQ=255)
                                                                                                                     gagagaTATCTTCTTCTGTCTTGTAACAGAATAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGc                                 <  2:1741909/149‑1 (MQ=255)
                                                                                                                     gagagaTATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACACAAGGTAAGc                                 >  2:4917465/1‑149 (MQ=255)
                                                                                                                     gagagaTATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGc                                 <  2:1925906/149‑1 (MQ=255)
                                                                                                                     gagagaTATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGc                                 <  1:1379331/149‑1 (MQ=255)
                                                                                                                     gagagaTATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGc                                 >  2:2358498/1‑149 (MQ=255)
                                                                                                                     gagagaTATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAAACGACGCGTCCGGTCAAACAAAAGGTAAGc                                 >  1:1655076/1‑149 (MQ=255)
                                                                                                                      agagaTATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGcc                                >  1:3889688/1‑149 (MQ=255)
                                                                                                                         gaTATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCACAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCcgtc                             >  2:444857/1‑149 (MQ=255)
                                                                                                                         gaTATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCcgtc                             >  2:444926/1‑149 (MQ=255)
                                                                                                                         gaTATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAACcgtc                             <  1:939962/149‑1 (MQ=255)
                                                                                                                             tcttcttctGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCcgtcgtcg                         >  1:392181/1‑149 (MQ=255)
                                                                                                                              cttcttctGTCTTGTAACCGCAGAACAGCAAATCGGGCTTTTGCGCCGTAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCcgtcgtcgt                        <  1:1543030/149‑1 (MQ=255)
                                                                                                                              cttcttctGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAACcgtcgtcgt                        >  2:1824967/1‑149 (MQ=255)
                                                                                                                               ttcttctGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAACCGTCGTCGtt                       <  2:1662427/149‑1 (MQ=255)
                                                                                                                                tcttctGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTa                      >  1:1453738/1‑149 (MQ=255)
                                                                                                                                tcttctGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTa                      <  2:545262/149‑1 (MQ=255)
                                                                                                                                 cttctGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTAc                     >  1:2673008/1‑149 (MQ=255)
                                                                                                                                 cttctGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTAc                     <  1:4377437/149‑1 (MQ=255)
                                                                                                                                 cttctGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAACCGTCGTCGTTAc                     >  1:376303/1‑149 (MQ=255)
                                                                                                                                  ttctGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACg                    >  2:523706/1‑149 (MQ=255)
                                                                                                                                      gTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGa                >  1:2739567/1‑149 (MQ=255)
                                                                                                                                      gTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGa                >  1:113845/1‑149 (MQ=255)
                                                                                                                                      gTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGa                >  1:2368809/1‑149 (MQ=255)
                                                                                                                                         ttGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGa             <  2:121533/149‑1 (MQ=255)
                                                                                                                                          tGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCTTCGTCGTTACGAAGATGAc            >  1:703587/1‑149 (MQ=255)
                                                                                                                                           gTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACg           >  2:1940168/1‑149 (MQ=255)
                                                                                                                                             aaCAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACgat         >  1:298579/1‑149 (MQ=255)
                                                                                                                                             aaCAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACgat         >  1:3709607/1‑149 (MQ=255)
                                                                                                                                              aCAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACgatg        <  1:1917915/149‑1 (MQ=255)
                                                                                                                                              aCAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACgatg        <  1:4843063/149‑1 (MQ=255)
                                                                                                                                               cAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACgatga       >  2:3899547/1‑149 (MQ=255)
                                                                                                                                                agaagaACAGAAAATCGGGCTTTTGCGCCGGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCACTTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACgatgat      <  2:2144395/149‑1 (MQ=255)
                                                                                                                                                 gaagaaCAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGAtt     >  1:3031109/1‑149 (MQ=255)
                                                                                                                                                 gaagaaCAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGAtt     >  2:2587103/1‑149 (MQ=255)
                                                                                                                                                  aagaaCAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTa    >  1:499691/1‑149 (MQ=255)
                                                                                                                                                   agaaCAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTac   <  2:3970522/149‑1 (MQ=255)
                                                                                                                                                   agaaCAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTac   >  2:1400177/1‑149 (MQ=255)
                                                                                                                                                    gaaCAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTacg  >  2:4283847/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACG  >  CP009273/161023‑161319

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: