Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,825,939 T→A 6.5% intergenic (‑314/+16) spy ← / ← astE periplasmic ATP‑independent protein refolding chaperone, stress‑induced/succinylglutamate desuccinylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,825,9390TA 6.5% 246.3 / 3.5 83intergenic (‑314/+16)spy/astEperiplasmic ATP‑independent protein refolding chaperone, stress‑induced/succinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base T (30/47);  new base A (3/3);  total (33/50)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCT  >  NC_000913/1825795‑1826088
                                                                                                                                                |                                                                                                                                                      
ggAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTg                                                                                                                                                   >  2:12740/1‑150 (MQ=255)
   gagTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCgg                                                                                                                                                <  1:1629105/150‑1 (MQ=255)
    agTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACTGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAAATCATTCTGCggg                                                                                                                                               <  2:85717/150‑1 (MQ=255)
    agTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCggg                                                                                                                                               <  2:1549299/150‑1 (MQ=255)
    agTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCggg                                                                                                                                               <  1:996977/150‑1 (MQ=255)
    agTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCggg                                                                                                                                               <  1:787531/150‑1 (MQ=255)
    agTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAAATCATTCTGCggg                                                                                                                                               <  1:1226697/150‑1 (MQ=255)
              taaTATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAg                                                                                                                                     <  2:661263/149‑1 (MQ=255)
                         cacGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTc                                                                                                                         >  1:1416448/1‑150 (MQ=255)
                         cacGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTc                                                                                                                         <  1:1493417/150‑1 (MQ=255)
                          acGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCg                                                                                                                        <  1:1695921/150‑1 (MQ=255)
                            gattatTTTGCAGAGATTATGAAGAACTATACCTGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCCAg                                                                                                                      <  1:1175308/150‑1 (MQ=255)
                            gatgatTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAg                                                                                                                      <  2:1471762/150‑1 (MQ=255)
                            gatgatTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGGGATACATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAg                                                                                                                      <  2:1349851/150‑1 (MQ=255)
                            gaagatTTTGCATACATTATTAAGAACTATACCGGAAGACTGGTGTTAAATAAAGCAATTAACCAGGATTATTCTGTATTAATTTATAAGAAAGCAACTTAATACGCGCAGAATGATTTCTGCGGGTATGTATTAGCTTATTTTTTCGAg                                                                                                                      <  1:1773448/147‑1 (MQ=255)
                                     gCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATcc                                                                                                             <  1:476620/150‑1 (MQ=255)
                                     gCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACAAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACACGAAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATcc                                                                                                             >  2:1031139/1‑150 (MQ=255)
                                         agaTTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCcgcg                                                                                                         <  2:834266/150‑1 (MQ=255)
                                         agaTTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCcgcg                                                                                                         <  2:604398/150‑1 (MQ=255)
                                            ttATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGt                                                                                                      >  2:959786/1‑150 (MQ=255)
                                               tGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAAt                                                                                                   >  2:1790692/1‑150 (MQ=255)
                                                  agaaCTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATccc                                                                                                <  2:721137/150‑1 (MQ=255)
                                                   gaaCTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCa                                                                                               >  2:286682/1‑150 (MQ=255)
                                                     aCTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCACAAc                                                                                             >  2:148945/1‑150 (MQ=255)
                                                       tataCCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGATTATTTTTTCGAGCATTAATCCCGAGCGTAATCACAACGc                                                                                           >  1:218249/1‑150 (MQ=255)
                                                         taCCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTa                                                                                         <  1:962444/150‑1 (MQ=255)
                                                           cctgATGACTGGTTATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTAcc                                                                                       <  1:132543/147‑1 (MQ=255)
                                                           ccGTATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCTCAGAATGATTTCTGTGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTAcc                                                                                       <  1:925354/150‑1 (MQ=255)
                                                           ccGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTAcc                                                                                       >  2:1287250/1‑150 (MQ=255)
                                                            cGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCACAATGATTTATGCGGGTAAGTATTAGGTTATTTTTTCGAGCATTAATACCCAGCGTAATCCCCACGCTACCa                                                                                      >  2:1128285/1‑150 (MQ=255)
                                                             ggATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCaa                                                                                     >  1:1362931/1‑150 (MQ=255)
                                                               aTGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACg                                                                                   <  2:1418824/150‑1 (MQ=255)
                                                                  aCTGTTGATAAATAAAGAAAATAACCATGATTAATCTGTATTAATTTAAAAGAAAGCAACTTAATACCCGAAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGAt                                                                                <  2:77265/150‑1 (MQ=255)
                                                                  aCTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGAt                                                                                >  2:1297878/1‑150 (MQ=255)
                                                                      gTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTTCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAgg                                                                            <  1:27111/150‑1 (MQ=255)
                                                                      gTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAgg                                                                            <  2:1681909/150‑1 (MQ=255)
                                                                      gTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAgg                                                                            >  1:904043/1‑150 (MQ=255)
                                                                       tGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTTATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAggg                                                                           <  2:1408416/150‑1 (MQ=255)
                                                                       tGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCTCGTAATCCCAACGCTACCAACTGATTAggg                                                                           <  2:72520/150‑1 (MQ=255)
                                                                        gATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGa                                                                          >  2:740536/1‑150 (MQ=255)
                                                                           aaataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATa                                                                       >  1:1644596/1‑150 (MQ=255)
                                                                           aaataaaGCAAATAACCAGGATTAATCTGAATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGAATTCTGCCGGTAAGTATTAGATTATTTTTTAGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATa                                                                       >  1:476622/1‑150 (MQ=255)
                                                                             ataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGAGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAac                                                                     <  1:760581/150‑1 (MQ=255)
                                                                              taaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAaca                                                                    <  1:513631/150‑1 (MQ=255)
                                                                                aaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAAcaca                                                                  <  1:515647/150‑1 (MQ=255)
                                                                                aaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAAcaca                                                                  <  1:1392440/150‑1 (MQ=255)
                                                                                aaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAAcaca                                                                  >  1:547392/1‑150 (MQ=255)
                                                                                  tcAAATACCCAGGATTAATCTGTATTAATTTATAAGAAATCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGAATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATa                                                                <  2:1291623/149‑1 (MQ=255)
                                                                                   cAAATAACCAGGAATAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCCGCGGGTAAGTATTAGCTTATTTTTTCTAGCATTAATCACTCGCGTAATCCCAACGCTACCAACGGATTATGGAATAACACATAc                                                               >  1:123562/1‑150 (MQ=255)
                                                                                      aTAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATActct                                                            <  2:226483/150‑1 (MQ=255)
                                                                                         aCCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTAca                                                         <  1:1277515/150‑1 (MQ=255)
                                                                                         aCCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTAca                                                         >  1:812551/1‑150 (MQ=255)
                                                                                          ccAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTAcat                                                        <  2:1788010/150‑1 (MQ=255)
                                                                                             ggATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTAcatcat                                                     <  1:1809710/150‑1 (MQ=255)
                                                                                             ggATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGAATTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTAcatcat                                                     <  1:516648/150‑1 (MQ=255)
                                                                                              gATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATg                                                    >  2:1531676/1‑150 (MQ=255)
                                                                                                    tCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTAc                                              <  2:1110088/150‑1 (MQ=255)
                                                                                                    tCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTAc                                              <  2:1782541/150‑1 (MQ=255)
                                                                                                     cTGTATTAATTTCTAAGCAAGCAACTTAATACCCGCAGCATGATTGCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACg                                             <  2:1268801/150‑1 (MQ=255)
                                                                                                       gTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGt                                           >  2:1564463/1‑150 (MQ=255)
                                                                                                         aTTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGATCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTaa                                         <  1:1394247/150‑1 (MQ=255)
                                                                                                          ttAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCTCGTAATCCCAACGCTACAAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTaaa                                        >  1:548064/1‑150 (MQ=255)
                                                                                                           tAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTaaaa                                       <  1:959786/150‑1 (MQ=255)
                                                                                                            aaTTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAAc                                      <  1:286682/150‑1 (MQ=255)
                                                                                                             aTTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACCCATACTCTACATCATGGGTTACGGTAAAACg                                     >  2:1509294/1‑150 (MQ=255)
                                                                                                                 ataAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGTTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTcc                                 <  2:338601/150‑1 (MQ=255)
                                                                                                                 ataAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTcc                                 >  2:441270/1‑150 (MQ=255)
                                                                                                                    aGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCctct                              >  2:1114388/1‑150 (MQ=255)
                                                                                                                      aaaGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTcc                            <  2:1196308/150‑1 (MQ=255)
                                                                                                                            aCTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTg                      <  1:1453283/150‑1 (MQ=255)
                                                                                                                            aCTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTACGTCCTg                      >  1:1349466/1‑150 (MQ=255)
                                                                                                                              ttAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTgcg                    <  1:362940/150‑1 (MQ=255)
                                                                                                                                aaTACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGcc                  <  2:702157/150‑1 (MQ=255)
                                                                                                                                    cccGCAGAAATCATTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCa              >  2:128103/1‑150 (MQ=255)
                                                                                                                                    cccGCAGAAATCATTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCAAGGGTTACGGTAAAACGTTCCTCTCCGTACTGCGCCAGCa              >  2:1595810/1‑150 (MQ=255)
                                                                                                                                      cGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAAt            >  1:1600645/1‑150 (MQ=255)
                                                                                                                                         aGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGt          <  1:12740/150‑1 (MQ=255)
                                                                                                                                          gAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTc        <  2:1416448/150‑1 (MQ=255)
                                                                                                                                           aaTGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTcc       <  1:1239336/150‑1 (MQ=255)
                                                                                                                                            aTGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCt      >  1:1816086/1‑150 (MQ=255)
                                                                                                                                            aTGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCt      >  1:395029/1‑150 (MQ=255)
                                                                                                                                              gATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCttt    >  2:799576/1‑150 (MQ=255)
                                                                                                                                               aTTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCGAAGGCTACCAACGGATTAGGGAATAACACATTCTCTACATCATGGGTTACGGTATAACGTTCCTCTCCGTCCTGAGCCAGCAATGTTCCttt  >  2:286790/1‑149 (MQ=255)
                                                                                                                                                tttCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCTTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATTTGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTtct  <  1:1031139/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
GGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCT  >  NC_000913/1825795‑1826088

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