Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,596,873 Δ1 bp 6.2% intergenic (‑136/+32) ypfJ ← / ← purC putative neutral zinc metallopeptidase/phosphoribosylaminoimidazole‑succinocarboxamide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,596,8730C. 6.2% 604.4 / 32.9 176intergenic (‑136/+32)ypfJ/purCputative neutral zinc metallopeptidase/phosphoribosylaminoimidazole‑succinocarboxamide synthetase
Reads supporting (aligned to +/- strand):  ref base C (70/94);  new base . (6/5);  total (76/100)
Fisher's exact test for biased strand distribution p-value = 5.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.96e-01

Only 100 of 178 total aligned reads displayed.

CTTGCCAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAG  >  NC_000913/2596725‑2597021
                                                                                                                                                    |                                                                                                                                                      
  tGCCAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCa                                                                                                                                                      <  2:1711841/149‑1 (MQ=255)
   gCCAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAg                                                                                                                                                     <  2:393530/149‑1 (MQ=255)
    ccAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGc                                                                                                                                                    <  1:280923/149‑1 (MQ=255)
    ccAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGc                                                                                                                                                    >  1:1754358/1‑149 (MQ=255)
    ccAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGc                                                                                                                                                    >  1:1549756/1‑149 (MQ=255)
      aaCGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑g                                                                                                                                                 >  1:1001775/1‑149 (MQ=255)
            tAGATACCTCTACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGc                                                                                                                                           <  1:1012155/149‑1 (MQ=255)
                tACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAga                                                                                                                                       <  1:800360/149‑1 (MQ=255)
                 aCCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAgat                                                                                                                                      >  3:269120/1‑149 (MQ=255)
                   cTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAgatga                                                                                                                                    >  4:1452142/1‑149 (MQ=255)
                   cTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAgatga                                                                                                                                    <  3:1159781/149‑1 (MQ=255)
                   cTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAgatga                                                                                                                                    >  1:335314/1‑149 (MQ=255)
                   cTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAgatga                                                                                                                                    >  3:94072/1‑149 (MQ=255)
                    tCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAgatgat                                                                                                                                   <  2:587952/149‑1 (MQ=255)
                    tCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAgatgat                                                                                                                                   >  1:307156/1‑149 (MQ=255)
                     caacaaTATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGatgatg                                                                                                                                  <  1:924404/149‑1 (MQ=255)
                     caacaaTATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGatgatg                                                                                                                                  <  1:858162/149‑1 (MQ=255)
                       acaaTATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGaa                                                                                                                                <  4:700764/149‑1 (MQ=255)
                        caaTATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCCGCT‑GCACGGCAAGATGATGAAc                                                                                                                               >  4:184734/1‑149 (MQ=255)
                          atatatTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAg                                                                                                                             <  1:69354/149‑1 (MQ=255)
                            atatTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑aa                                                                                                                          <  3:1215107/149‑1 (MQ=255)
                            atatTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑aa                                                                                                                          <  3:1590130/149‑1 (MQ=255)
                              atTCATTATGGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCTGCT‑GCACGGCAAGATGATGAACAG‑aaaa                                                                                                                        >  4:1482802/1‑149 (MQ=255)
                               tgcATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑aaaaa                                                                                                                       <  3:1376901/147‑1 (MQ=255)
                                  aTTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCa                                                                                                                    <  3:25607/149‑1 (MQ=255)
                                  aTTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCa                                                                                                                    <  1:434767/149‑1 (MQ=255)
                                  aTTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCa                                                                                                                    <  3:181623/149‑1 (MQ=255)
                                  aTTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCa                                                                                                                    <  1:1181107/149‑1 (MQ=255)
                                  aTTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTACAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCa                                                                                                                    <  3:1426836/149‑1 (MQ=255)
                                    tataGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGt                                                                                                                  <  2:1754358/149‑1 (MQ=255)
                                    tataGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGt                                                                                                                  <  2:1549756/149‑1 (MQ=255)
                                     ataGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTc                                                                                                                 <  1:100369/149‑1 (MQ=255)
                                      taGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCa                                                                                                                    <  3:430116/145‑1 (MQ=255)
                                      taGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCa                                                                                                                    >  4:430116/1‑145 (MQ=255)
                                      taGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTcc                                                                                                                >  3:789111/1‑149 (MQ=255)
                                          aaTGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCt                                                                                                            <  1:1231318/149‑1 (MQ=255)
                                          aaTGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCt                                                                                                            >  3:1047095/1‑149 (MQ=255)
                                          aaTGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCt                                                                                                            >  1:989710/1‑149 (MQ=255)
                                           aTGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTg                                                                                                           <  3:743206/149‑1 (MQ=255)
                                           aTGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTAAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTg                                                                                                           >  2:232705/1‑149 (MQ=255)
                                              aTCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTac                                                                                                        >  2:410641/1‑149 (MQ=255)
                                               tCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAg                                                                                                                   <  1:1576784/137‑1 (MQ=255)
                                               tCGTAGGTGGTACAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGCGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTaca                                                                                                       >  1:1431914/1‑149 (MQ=255)
                                                 gTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACAcc                                                                                                     <  1:16189/149‑1 (MQ=255)
                                                 gTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACAcc                                                                                                     >  2:958729/1‑149 (MQ=255)
                                                         ttCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAggcggc                                                                                             >  1:1660969/1‑149 (MQ=255)
                                                          tCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCAAg                                                                                                                                        >  4:758319/1‑106 (MQ=255)
                                                          tCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCAAg                                                                                                                                        <  3:758319/106‑1 (MQ=255)
                                                           cAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGcggcgg                                                                                           <  4:820307/149‑1 (MQ=255)
                                                            agagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCaa                                                                                                                                         >  4:387357/1‑103 (MQ=255)
                                                             gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCggcgggc                                                                                         >  1:904909/1‑149 (MQ=255)
                                                              agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCggcgggcg                                                                                        >  1:606308/1‑149 (MQ=255)
                                                               ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGa                                                                                       >  2:1183564/1‑149 (MQ=255)
                                                                gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGAc                                                                                      <  2:1272586/149‑1 (MQ=255)
                                                                   gATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATCAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGc                                                                                   <  4:136608/149‑1 (MQ=255)
                                                                       aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCa                                                                               <  2:1365917/149‑1 (MQ=255)
                                                                        ccACAGGAAAGCGGGATGAAATGGCGAAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCAt                                                                              <  3:551929/149‑1 (MQ=255)
                                                                         cacaGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCTGGCGGCGGGCGACGGCTTCATa                                                                             >  1:1506050/1‑149 (MQ=255)
                                                                         cacaGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                                                                             >  3:1444793/1‑149 (MQ=255)
                                                                          acaGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAg                                                                            <  4:1201082/149‑1 (MQ=255)
                                                                            aGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGt                                                                                                                  >  2:132559/1‑109 (MQ=255)
                                                                              gAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCtt                                                                        >  4:272636/1‑149 (MQ=255)
                                                                                  gCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGAt                                                                    <  2:1012048/149‑1 (MQ=255)
                                                                                    gggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCa                                                                  <  1:232705/149‑1 (MQ=255)
                                                                                     ggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAg                                                                 <  3:1151776/149‑1 (MQ=255)
                                                                                      gATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAAATCAGTCCACCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGACCAg                                                                 >  3:533340/1‑149 (MQ=255)
                                                                                      gATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGAt                                                                    >  4:1556821/1‑145 (MQ=255)
                                                                                      gATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAgg                                                                <  2:1005386/149‑1 (MQ=255)
                                                                                       aTGAAATGGCGGAGGAATGAGCGGTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGc                                                               >  4:10741/1‑149 (MQ=255)
                                                                                              ggCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCGGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCa                                                                                                                    >  2:522277/1‑89 (MQ=255)
                                                                                              ggCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCa                                                                                                                    <  1:522277/89‑1 (MQ=255)
                                                                                              ggCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                                                                             <  2:1089894/128‑1 (MQ=255)
                                                                                              ggCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                                                                             >  1:1089894/1‑128 (MQ=255)
                                                                                                   aggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGctct                                                   <  4:789111/149‑1 (MQ=255)
                                                                                                     gaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTgg                                                 <  3:442120/149‑1 (MQ=255)
                                                                                                      aaTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGc                                                >  2:810017/1‑149 (MQ=255)
                                                                                                       aTGAGCGTTCGATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACTGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCAGAGGCTTCGATCAGGCCACCGAGGCTCTGGCg                                               >  3:379436/1‑149 (MQ=255)
                                                                                                           gCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGaaa                                           <  2:921163/149‑1 (MQ=255)
                                                                                                           gCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGaaa                                           >  4:231359/1‑149 (MQ=255)
                                                                                                           gCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGaaa                                           <  2:1535840/149‑1 (MQ=255)
                                                                                                              ttCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACgg                                        >  4:1625611/1‑149 (MQ=255)
                                                                                                              ttCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACgg                                        >  2:1461599/1‑149 (MQ=255)
                                                                                                              aacAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCAGGGTAAGATGATGAAGAG‑AAAAATCTGTCCAGCTGGACGCCCAGGCGGCGGGCGACGGCCTCAAAGGCTTCGATCAGGCCAACGAGGCTCCGGCGGAGACga                                        >  2:1042617/3‑148 (MQ=255)
                                                                                                                cAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTc                                      <  1:593778/149‑1 (MQ=255)
                                                                                                                   tGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATTAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGc                                                                                   <  2:833027/101‑1 (MQ=255)
                                                                                                                    gCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTg                                  <  2:1579137/149‑1 (MQ=255)
                                                                                                                      atatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑aa                                                                                                                          <  2:903096/59‑1 (MQ=255)
                                                                                                                       tatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGa                                                         >  3:89711/1‑123 (MQ=255)
                                                                                                                        atatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCa                              >  4:957756/1‑149 (MQ=255)
                                                                                                                         tatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCAt                             <  1:1042617/149‑1 (MQ=255)
                                                                                                                              cGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTcc                               <  1:955173/142‑1 (MQ=255)
                                                                                                                                      gTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGtt                >  2:632964/1‑149 (MQ=255)
                                                                                                                                       ttGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGttt               >  1:1386653/1‑149 (MQ=255)
                                                                                                                                         gTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCt             <  3:957756/149‑1 (MQ=255)
                                                                                                                                         gTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCt             <  3:681126/149‑1 (MQ=255)
                                                                                                                                         gTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCt             <  3:1541521/149‑1 (MQ=255)
                                                                                                                                          tCTTGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                                                                             <  2:208369/84‑1 (MQ=38)
                                                                                                                                            ttGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTg          >  3:60787/1‑149 (MQ=255)
                                                                                                                                             tGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGt         <  2:286806/149‑1 (MQ=255)
                                                                                                                                              gCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTc        <  1:1461049/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CTTGCCAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAG‑AAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAG  >  NC_000913/2596725‑2597021

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑