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

GGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCT  >  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)
                         cacGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTc                                                                                                                         >  1:1416448/1‑150 (MQ=255)
                         cacGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTc                                                                                                                         <  1:1493417/150‑1 (MQ=255)
                          acGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCg                                                                                                                        <  1:1695921/150‑1 (MQ=255)
                            gattatTTTGCAGAGATTATGAAGAACTATACCTGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCCAg                                                                                                                      <  1:1175308/150‑1 (MQ=255)
                            gatgatTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAg                                                                                                                      <  2:1471762/150‑1 (MQ=255)
                            gatgatTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGGGATACATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAg                                                                                                                      <  2:1349851/150‑1 (MQ=255)
                            gaagatTTTGCATACATTATTAAGAACTATACCGGAAGACTGGTGTTAAATAAAGCAATTAACCAGGATTATTCTGTATTAATTTATAAGAAAGCAACTTAATACGCGCAGAATGATTTCTGCGGGTATGTATTAGCTTA‑TTTTTTCGAg                                                                                                                      <  1:1773448/147‑1 (MQ=255)
                                     gCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATcc                                                                                                             <  1:476620/150‑1 (MQ=255)
                                     gCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACAAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACACGAAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATcc                                                                                                             >  2:1031139/1‑150 (MQ=255)
                                         agaTTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCcgcg                                                                                                         <  2:834266/150‑1 (MQ=255)
                                         agaTTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCcgcg                                                                                                         <  2:604398/150‑1 (MQ=255)
                                            ttATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGt                                                                                                      >  2:959786/1‑150 (MQ=255)
                                               tGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAAt                                                                                                   >  2:1790692/1‑150 (MQ=255)
                                                  agaaCTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATccc                                                                                                <  2:721137/150‑1 (MQ=255)
                                                   gaaCTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCa                                                                                               >  2:286682/1‑150 (MQ=255)
                                                     aCTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCACAAc                                                                                             >  2:148945/1‑150 (MQ=255)
                                                       tataCCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGATTA‑TTTTTTCGAGCATTAATCCCGAGCGTAATCACAACGc                                                                                           >  1:218249/1‑150 (MQ=255)
                                                         taCCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTa                                                                                         <  1:962444/150‑1 (MQ=255)
                                                           cctgATGACTGGTTATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTAcc                                                                                       <  1:132543/147‑1 (MQ=255)
                                                           ccGTATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCTCAGAATGATTTCTGTGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTAcc                                                                                       <  1:925354/150‑1 (MQ=255)
                                                           ccGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTAcc                                                                                       >  2:1287250/1‑150 (MQ=255)
                                                            cGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCACAATGATTTATGCGGGTAAGTATTAGGTTA‑TTTTTTCGAGCATTAATACCCAGCGTAATCCCCACGCTACCa                                                                                      >  2:1128285/1‑150 (MQ=255)
                                                             ggATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCaa                                                                                     >  1:1362931/1‑150 (MQ=255)
                                                               aTGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACg                                                                                   <  2:1418824/150‑1 (MQ=255)
                                                                  aCTGTTGATAAATAAAGAAAATAACCATGATTAATCTGTATTAATTTAAAAGAAAGCAACTTAATACCCGAAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGAt                                                                                <  2:77265/150‑1 (MQ=255)
                                                                  aCTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGAt                                                                                >  2:1297878/1‑150 (MQ=255)
                                                                      gTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTTCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAgg                                                                            <  1:27111/150‑1 (MQ=255)
                                                                      gTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAgg                                                                            <  2:1681909/150‑1 (MQ=255)
                                                                      gTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAgg                                                                            >  1:904043/1‑150 (MQ=255)
                                                                       tGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTTATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAggg                                                                           <  2:1408416/150‑1 (MQ=255)
                                                                       tGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCTCGTAATCCCAACGCTACCAACTGATTAggg                                                                           <  2:72520/150‑1 (MQ=255)
                                                                        gATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGa                                                                          >  2:740536/1‑150 (MQ=255)
                                                                           aaataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATa                                                                       >  1:1644596/1‑150 (MQ=255)
                                                                           aaataaaGCAAATAACCAGGATTAATCTGAATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGAATTCTGCCGGTAAGTATTAGATTA‑TTTTTTAGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATa                                                                       >  1:476622/1‑150 (MQ=255)
                                                                             ataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGAGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAac                                                                     <  1:760581/150‑1 (MQ=255)
                                                                              taaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAaca                                                                    <  1:513631/150‑1 (MQ=255)
                                                                                aaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAAcaca                                                                  <  1:515647/150‑1 (MQ=255)
                                                                                aaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAAcaca                                                                  <  1:1392440/150‑1 (MQ=255)
                                                                                aaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAAcaca                                                                  >  1:547392/1‑150 (MQ=255)
                                                                                  tcAAATACCCAGGATTAATCTGTATTAATTTATAAGAAATCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGAATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATa                                                                <  2:1291623/149‑1 (MQ=255)
                                                                                   cAAATAACCAGGAATAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCCGCGGGTAAGTATTAGCTTA‑TTTTTTCTAGCATTAATCACTCGCGTAATCCCAACGCTACCAACGGATTATGGAATAACACATAc                                                               >  1:123562/1‑150 (MQ=255)
                                                                                      aTAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATActct                                                            <  2:226483/150‑1 (MQ=255)
                                                                                         aCCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTAca                                                         <  1:1277515/150‑1 (MQ=255)
                                                                                         aCCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTAca                                                         >  1:812551/1‑150 (MQ=255)
                                                                                          ccAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTAcat                                                        <  2:1788010/150‑1 (MQ=255)
                                                                                             ggATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTAcatcat                                                     <  1:1809710/150‑1 (MQ=255)
                                                                                             ggATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGAATTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTAcatcat                                                     <  1:516648/150‑1 (MQ=255)
                                                                                              gATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATg                                                    >  2:1531676/1‑150 (MQ=255)
                                                                                                    tCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTAc                                              <  2:1110088/150‑1 (MQ=255)
                                                                                                    tCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTAc                                              <  2:1782541/150‑1 (MQ=255)
                                                                                                     cTGTATTAATTTCTAAGCAAGCAACTTAATACCCGCAGCATGATTGCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACg                                             <  2:1268801/150‑1 (MQ=255)
                                                                                                       gTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGt                                           >  2:1564463/1‑150 (MQ=255)
                                                                                                         aTTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGATCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTaa                                         <  1:1394247/150‑1 (MQ=255)
                                                                                                          ttAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCTCGTAATCCCAACGCTACAAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTaaa                                        >  1:548064/1‑150 (MQ=255)
                                                                                                           tAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTaaaa                                       <  1:959786/150‑1 (MQ=255)
                                                                                                            aaTTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAAc                                      <  1:286682/150‑1 (MQ=255)
                                                                                                             aTTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACCCATACTCTACATCATGGGTTACGGTAAAACg                                     >  2:1509294/1‑150 (MQ=255)
                                                                                                                 ataAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGTTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTcc                                 <  2:338601/150‑1 (MQ=255)
                                                                                                                 ataAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTcc                                 >  2:441270/1‑150 (MQ=255)
                                                                                                                    aGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCctct                              >  2:1114388/1‑150 (MQ=255)
                                                                                                                      aaaGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTcc                            <  2:1196308/150‑1 (MQ=255)
                                                                                                                            aCTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTg                      <  1:1453283/150‑1 (MQ=255)
                                                                                                                            aCTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTACGTCCTg                      >  1:1349466/1‑150 (MQ=255)
                                                                                                                              ttAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTgcg                    <  1:362940/150‑1 (MQ=255)
                                                                                                                                aaTACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGcc                  <  2:702157/150‑1 (MQ=255)
                                                                                                                                    cccGCAGAAATCATTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCa              >  2:128103/1‑150 (MQ=255)
                                                                                                                                    cccGCAGAAATCATTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCAAGGGTTACGGTAAAACGTTCCTCTCCGTACTGCGCCAGCa              >  2:1595810/1‑150 (MQ=255)
                                                                                                                                      cGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAAt            >  1:1600645/1‑150 (MQ=255)
                                                                                                                                         aGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGt          <  1:12740/150‑1 (MQ=255)
                                                                                                                                          gAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTc        <  2:1416448/150‑1 (MQ=255)
                                                                                                                                           aaTGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTcc       <  1:1239336/150‑1 (MQ=255)
                                                                                                                                            aTGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCt      >  1:1816086/1‑150 (MQ=255)
                                                                                                                                            aTGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCt      >  1:395029/1‑150 (MQ=255)
                                                                                                                                              gATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCttt    >  2:799576/1‑150 (MQ=255)
                                                                                                                                               aTTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCGAAGGCTACCAACGGATTAGGGAATAACACATTCTCTACATCATGGGTTACGGTATAACGTTCCTCTCCGTCCTGAGCCAGCAATGTTCCtttg   >  2:286790/1‑149 (MQ=255)
                                                                                                                                                tttCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCTTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATTTGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTtct  <  1:1031139/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
GGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTA‑TTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCT  >  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: ‑