Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,311,748 G→A 5.3% intergenic (‑210/+68) PP_1144 ← / ← rapA membrane protein/RNA polymerase‑associated protein RapA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,311,7480GA 5.3% 287.2 / 4.3 99intergenic (‑210/+68)PP_1144/rapAmembrane protein/RNA polymerase‑associated protein RapA
Reads supporting (aligned to +/- strand):  ref base G (33/59);  new base A (2/3);  total (37/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.81e-01

GGCAGTATAGGAAGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTT  >  NC_002947/1311600‑1311896
                                                                                                                                                    |                                                                                                                                                     
ggCAGTATAGGAAGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCgg                                                                                                                                                      <  2:1111607/150‑1 (MQ=255)
   aGTATAGGAAGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTgg                                                                                                                                                   <  2:1999367/150‑1 (MQ=255)
      ataGGAAGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATg                                                                                                                                               <  2:2547990/150‑1 (MQ=255)
      ataGGAAGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACCGCAGTGGATg                                                                                                                                               <  1:3600951/150‑1 (MQ=255)
      ataGGAAGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACCGCAGTGGATg                                                                                                                                               <  2:5356150/150‑1 (MQ=255)
       taGGAAGCACTACTTAAGCCTTCCTTAATAATGGAACAATTATTTATACTATTTTGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGCCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACCGCAGTGGATGc                                                                                                                                              <  1:2452343/150‑1 (MQ=255)
         ggaaggaGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTATTGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCa                                                                                                                                            <  1:1017674/150‑1 (MQ=255)
         ggaaggaGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCa                                                                                                                                            <  1:1017730/150‑1 (MQ=255)
            aggaGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGcttcgc                                                                                                                                         <  1:4909337/150‑6 (MQ=255)
                 tACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGTTTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCttt                                                                                                                                    <  1:3910755/150‑1 (MQ=255)
                   cTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGc                                                                                                                                  <  1:382729/150‑1 (MQ=255)
                    ttAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGCAGCAACCACAGCCGGGGGTGGCCCCCGCTTTGcc                                                                                                                                 >  1:2636565/1‑150 (MQ=255)
                    ttAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATTCCACCGGCGTTTcc                                                                                                                                 >  2:4806724/1‑150 (MQ=255)
                    ttAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCGAAGCCTCCCCAGCGGGGGATGCCCCCCGCGTTGcc                                                                                                                                 >  2:3130097/1‑150 (MQ=255)
                     tAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCg                                                                                                                                <  1:5440910/150‑1 (MQ=255)
                     tAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCg                                                                                                                                <  1:3876513/150‑1 (MQ=255)
                       aGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGTAATAATTTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGgt                                                                                                                              <  1:2006119/150‑1 (MQ=255)
                       aGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGgt                                                                                                                              <  2:3080794/150‑1 (MQ=255)
                        gCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATACCTCCAATCGCAGTGAGTAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCAAAGCCGgtg                                                                                                                             <  1:1362005/150‑1 (MQ=255)
                        gCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCAAAGCCGgtg                                                                                                                             <  1:1361999/150‑1 (MQ=255)
                        gCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCAAATCGCCGTGAGGAATAATTTTGGGGCCCTCCTGCCTTCCTGTGGGCGCGGGCGCGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCAAAGCCGgtg                                                                                                                             <  1:2456883/150‑1 (MQ=255)
                         ccTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGGCCGCCAACCATACACAGCGGTGGGTGCCACCCGCGGTGGCGgtgt                                                                                                                            >  2:3412657/1‑150 (MQ=255)
                           ttCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCCGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCAAAGCCGgtggca                                                                                                                          <  1:1983038/150‑4 (MQ=255)
                            tCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGCACCAACCACGGCGGGGGTTGCCACCCGCGGTGCCGGTGTTcg                                                                                                                         >  2:522601/1‑150 (MQ=255)
                             ccTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGTAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCAAAGCCGgtggcatc                                                                                                                        <  2:3175099/150‑6 (MQ=255)
                              cTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCAACCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTcgcg                                                                                                                       >  2:2664898/1‑150 (MQ=255)
                               tGACTATTGGCACAATGATTTCTACTATTTGTTTATAACGCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCgg                                                                                                                      <  1:237143/150‑1 (MQ=255)
                               tGACTAATGGCACAATTATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGTAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCgg                                                                                                                      <  1:2428590/150‑1 (MQ=255)
                                 aCTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAACAATCCACATAGGTGGTATGCAACCGCGTTTCCGGGGTTAGCgggg                                                                                                                    >  2:59314/1‑149 (MQ=255)
                                       ggCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgccc                                                                                                              >  1:1679751‑M2/1‑146 (MQ=255)
                                        gCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCCCAGCCACCCCAGCGGTGGGTGGCCCCCGCGTTGCCGGTGTTAGCGGGCGGgcccg                                                                                                             <  2:2749867‑M2/150‑6 (MQ=255)
                                            caggattTCTACGAGTGGGTTCTAACGCCAATCGCAGTGAGGAATAACTGTGGGGCCCTAATGCCGTCCTGGGGGAGCGGGCGAGCCCGCGACGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGGGTTCGCGGGCGTgcccgctcc                                                                                                         >  2:4717758‑M2/4‑141 (MQ=255)
                                              tGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCAACCACAGGGGGGGGAGGCCCCCGCGTTTCCGGGGTTCGCGGGCGGGCCCCCGcccc                                                                                                       >  2:3684185/1‑149 (MQ=255)
                                                aTTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCCAAGCAAACACAGCGGGGGATGCCACCCGCTTTTCCGGTGTTCGCGGGCGGgcccgcgcccaca                                                                                                     <  2:2307597‑M2/150‑14 (MQ=255)
                                                   tCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGCAGCAACCACTGCGGTGGGTGCCCCCCGCGTGGCCCGTGTTCGCGGGCGTgcccccgcccaccagg                                                                                                  <  1:4839563‑M2/150‑17 (MQ=255)
                                                   tCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAgg                                                                                                  <  2:752119/150‑1 (MQ=255)
                                                    cTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAggg                                                                                                 <  2:4467907/150‑1 (MQ=255)
                                                    cTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCGAAGCATACACAGCGGTGGATGCCACCGGCGGTGCCGGGGTTCGCCGGCGGGCCCGCTCCCACCAggg                                                                                                 >  2:4750366/1‑150 (MQ=255)
                                                     tACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCg                                                                                                                                >  1:4387704/1‑118 (MQ=255)
                                                     tACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCg                                                                                                                                <  2:4387704/118‑1 (MQ=255)
                                                     tACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCAAACAATGCGGGGGGTGCCACCCGCTTTTCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAgggg                                                                                                >  2:3289467/1‑150 (MQ=38)
                                                     tACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCGCAGCAACCACAGCGGGGGATGCCCCCCGCGTTGCCCGGGTTCGCCGGCGGGCCCCCTCCCCCAAcggg                                                                                                >  2:3586603/1‑150 (MQ=255)
                                                          tttGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCgg                                                                                           <  2:4763931/150‑1 (MQ=255)
                                                          tttGGTTATAACTCCAATCGCAGGGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCgg                                                                                           <  2:1979826/150‑1 (MQ=255)
                                                             ggTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCCAAGCAAACAAAGCGGGGGATTCCACCCGCGTTGGCCGGGTTCGCGGGCGTGCCCCCTCCCACAACGGGCGCGGGGc                                                                                        >  2:1841965/1‑150 (MQ=255)
                                                                tataACCCCAATCGCCGTGAGTAATAATTTTGGGGCCCTAAGGCCTTCGTGTGGGCGCGGGCGCGCCCGCGACGCATCCACTGCGGGGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTg                                                                                     <  1:907913/150‑1 (MQ=25)
                                                                     cTCCAATCGCAGTGCGGAATAAATTTGGGGCCCTAAGGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTa                                                                                <  1:3577010/150‑1 (MQ=255)
                                                                      tCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCCAAGCATACACAGCGGGGGTTGCCACCCGCGTTGCCGGTGTTCGCGGGCGTGCCCCCTCCCCAAACGGGCGCGGGGCCGGCGCTag                                                                               >  2:3068148/1‑150 (MQ=37)
                                                                      tCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCCGACCAACCACACCCGGGGCTGCCCCCCGCGTTGCCCGGGGTCGCCGGCCGGGCCCCCCCCCCAACGGGCGCCGGGCCAgccccac                                                                               >  2:1444373/1‑145 (MQ=255)
                                                                      tCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCCAAACAAACCCAGCGGGGGGTTGCACCCCCTTTTCCGGGGGTTGCGGGCGGGCCCCCGCCCCCAACGGGCGCCGGGCCAgcccaag                                                                               >  1:432941/1‑145 (MQ=255)
                                                                       ccAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCAACCACAGCGGGGGATTCCACCGGCGTTGCCGGGGGTCGCGGGCGTGCCCCCTCCCACAAGGGGCGCGGGGCCAGCGCTagc                                                                              >  2:1281256/1‑150 (MQ=38)
                                                                       ccAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCGAAGCATACACAGCGGTGGATGCCACCCGCGTTGCCCGTGTTAGCGGGCGTGCCCCCTCCCACAAGGGGCGCGGGGCCAGCGCTagc                                                                              >  2:1280811/1‑150 (MQ=255)
                                                                         aaTCGCAGTGAGGAATAAATTTGGGGCCCTCATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTagcca                                                                            <  2:4224767/150‑1 (MQ=255)
                                                                             gCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGCGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAgccagcca                                                                        <  2:4347960/150‑1 (MQ=255)
                                                                                gTGAGGAATAATTTTGGGGCCCTACGGCCTTCGTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGccacca                                                                     <  1:3220572/150‑1 (MQ=255)
                                                                                      aaTAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCAcgc                                                               <  2:2208840/150‑1 (MQ=255)
                                                                                        tAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctccc                                                                                                        <  2:353722‑M2/107‑11 (MQ=255)
                                                                                        tAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctccc                                                                                                        >  1:353722‑M2/1‑97 (MQ=255)
                                                                                          attttGGGGCCCTTTTGCCTTCGTGTGGGAGGGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAAt                                                           <  1:4750366/148‑1 (MQ=255)
                                                                                          aaTTTGGGGCCCTAATGCCTTCCTGTGGGAGGGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAAt                                                           <  1:2822558/150‑1 (MQ=255)
                                                                                          aaTTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAAt                                                           <  2:2117132/150‑1 (MQ=255)
                                                                                          aaTTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAAt                                                           <  2:5129936/150‑1 (MQ=255)
                                                                                                    ccTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCCCTGCGGGGGGTGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctcccacaggcaggccttagggccccaaatttattcctcacggcgattggcggtctaacc                                                 <  2:4152464‑M2/150‑66 (MQ=255)
                                                                                                             ttCCTGTGGGCGCGGGCGCGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCGAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGggcg                                        <  1:2740760/150‑1 (MQ=255)
                                                                                                             ttCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGggcg                                        <  2:1996645/150‑1 (MQ=255)
                                                                                                             ttCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGggcg                                        <  2:3936505/150‑1 (MQ=255)
                                                                                                               cgtgtgGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGggcggc                                      <  1:1354479/148‑1 (MQ=255)
                                                                                                               ccTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGggcggc                                      <  1:4267959/150‑1 (MQ=255)
                                                                                                               ccTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTGGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGggcggc                                      <  1:4268126/150‑1 (MQ=255)
                                                                                                                cggtgGGAGCGGGGGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGcc                                     <  2:3477877/148‑1 (MQ=255)
                                                                                                                   tgGGAGCGGGCGAGCCCGCCAAGCATACACTGCGGTGGATGCCCCCGGCATTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCCGGGCCAGCCCAAGCCACCCACCAACACGCGAATTGCATACAGGCGCCGGGCGGCCGtc                                  >  2:2235594/1‑150 (MQ=25)
                                                                                                                        gCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctcccacaggaaggcattagggccccaaatttattcctcactgcgattggagttataaccaaatcgtagaaatcattgtg                             <  2:4627715‑M2/150‑86 (MQ=255)
                                                                                                                            gCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctcccacaggaaggcattagggccccaaatttattcctcactgcgattggagttataaccaaatcgtagaaatcattgtgccat                         <  1:4648223‑M2/150‑90 (MQ=255)
                                                                                                                            gCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCAtt                         >  2:2952054/1‑150 (MQ=255)
                                                                                                                            gCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCGTCCAGGCGCAGGGCGGCCTTCTCCAGCAtt                         >  2:4696499/1‑150 (MQ=255)
                                                                                                                               aGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctcccacaggaaggcattagggccccaaatttattcctcactgcgattggagttataaccaaatcgtagaaatcattgtgccattag                      <  1:3185955‑M2/150‑93 (MQ=255)
                                                                                                                                 cccGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAg                    >  2:1306251/1‑150 (MQ=255)
                                                                                                                                 cccGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAg                    >  2:3748696/1‑150 (MQ=255)
                                                                                                                                  ccGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctcccacaaggggcgcggggcctcaaatttattcctcactgcgattggagttataaccaaatcgtagaaatcattgtgccattagtcagg                   <  1:381620‑M2/150‑96 (MQ=255)
                                                                                                                                  ccGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTGCCAGGCGCAGGGCGGCCTTCTCCCGCATTGCCAgg                   >  2:49804/1‑150 (MQ=255)
                                                                                                                                   cgcgAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGc                  <  1:3480590/150‑1 (MQ=255)
                                                                                                                                   cgcgAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGc                  >  2:3480590/1‑150 (MQ=255)
                                                                                                                                        aGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctcccacaggaaggcattagggccccaaatttattcctcactgcgattggagttataaccaaatcgtagaaatcattgtgccattagtcaggaagg             <  1:4995665‑M2/150‑102 (MQ=255)
                                                                                                                                         gCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCa                              <  2:4746595/132‑1 (MQ=255)
                                                                                                                                         gCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCa                              >  1:4746595/1‑132 (MQ=255)
                                                                                                                                         gCATCCACCGCAGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCtt            >  1:4219828/1‑150 (MQ=255)
                                                                                                                                          cATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTagcc                                                                             >  2:2390496/1‑84 (MQ=255)
                                                                                                                                          cATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTagcc                                                                             <  1:2390496/84‑1 (MQ=255)
                                                                                                                                          cATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTc           <  2:2561006/150‑1 (MQ=255)
                                                                                                                                          cATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTc           <  2:5127644/150‑1 (MQ=255)
                                                                                                                                          cATCCACGGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctcccacaggaaggcattagggccccaaatttattcctcactgcgattggagttataaccaaatcgtagaaatcattgtgccattagtcaggaaggct           >  1:4152464‑M2/1‑47 (MQ=255)
                                                                                                                                           aTCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTcg          >  2:1792957/1‑150 (MQ=255)
                                                                                                                                            tCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTcgc         <  1:2304984/150‑1 (MQ=255)
                                                                                                                                             ccACCGCAGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTcgcg        >  1:841259/1‑150 (MQ=255)
                                                                                                                                              cACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGcgcg                                                                                            >  1:4359932/1‑65 (MQ=255)
                                                                                                                                              cACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGcgcg                                                                                            <  2:4359944/65‑1 (MQ=255)
                                                                                                                                              cACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGcgcg                                                                                            >  1:4359944/1‑65 (MQ=255)
                                                                                                                                              cACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGcgcg                                                                                            <  2:4359932/65‑1 (MQ=255)
                                                                                                                                                 tGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctcccacaggaaggcattagggccccaaatttattcctcactgcgattggagttataaccaaatcgtagaaatcattgtgccattagtcaggaaggcttaagtac    <  1:5377029‑M2/150‑111 (MQ=255)
                                                                                                                                                   cGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCtt  >  1:1657228/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGCAGTATAGGAAGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTT  >  NC_002947/1311600‑1311896

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