Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,596,873 Δ1 bp 5.1% 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. 5.1% 333.3 / 10.0 98intergenic (‑136/+32)ypfJ/purCputative neutral zinc metallopeptidase/phosphoribosylaminoimidazole‑succinocarboxamide synthetase
Reads supporting (aligned to +/- strand):  ref base C (49/44);  new base . (3/2);  total (52/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

AACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACA  >  NC_000913/2596731‑2597020
                                                                                                                                              |                                                                                                                                                    
aaCGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑g                                                                                                                                               >  8:58075/1‑149 (MQ=255)
  cGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCa                                                                                                                                             <  2:360356/149‑1 (MQ=255)
    cATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACg                                                                                                                                           >  5:191117/1‑149 (MQ=255)
     aTAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACgg                                                                                                                                          >  8:51617/1‑149 (MQ=255)
     aTAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGCTTGTCTTGCCGTGCCACCT‑GCACgg                                                                                                                                          >  1:210373/1‑149 (MQ=255)
         aTACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGc                                                                                                                                         >  5:61033/1‑146 (MQ=255)
            ccTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAgat                                                                                                                                    >  8:161736/1‑148 (MQ=255)
             cTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGc                                                                                                                                                  <  1:273875/134‑1 (MQ=255)
             cTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGc                                                                                                                                                  >  2:273875/1‑134 (MQ=255)
             cTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGGTGGCTTGCCGTGCCAGCT‑GCACGGCAAgatga                                                                                                                                  <  3:169504/149‑1 (MQ=255)
             cTCAACAATATATTCATTATAGAAATGATCGGAGGTGGTTCAGAGGAAGATTACCACACGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATAAATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAgatga                                                                                                                                  <  3:310392/149‑1 (MQ=255)
              tCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCaa                                                                                                                                       >  8:309068/1‑143 (MQ=255)
              tCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCaa                                                                                                                                       >  8:314361/1‑143 (MQ=255)
              tCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCaa                                                                                                                                       <  7:309068/143‑1 (MQ=255)
              tCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCaa                                                                                                                                       <  7:314361/143‑1 (MQ=255)
                  caaTATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGaa                                                                                                                              >  7:300521/1‑148 (MQ=255)
                   aaTATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATTGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGATTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACa                                                                                                                            <  3:60390/149‑1 (MQ=255)
                         ttCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGaaaaa                                                                                                                      >  3:175292/1‑149 (MQ=255)
                            aTTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGATGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAAt                                                                                                                     >  1:177874/1‑147 (MQ=255)
                            aTTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCa                                                                                                                   <  3:119645/149‑1 (MQ=255)
                            aTTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCa                                                                                                                   <  6:61033/149‑1 (MQ=255)
                             ttATAGAAATGATCGTAGGTTGTTCAGACGAATATTACCACAGTAAAGCGGGATTAAATGGCGGAGGAATGATCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAAGCT‑GCACGGCAAGATGATGAACAGAAAAATAAg                                                                                                                  >  3:23781/1‑149 (MQ=255)
                                taGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGTAGGAATGAGCGTTCAATGCATATAATCCGAGAGCGGTTGTCTTTCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTc                                                                                                                >  7:331701/1‑148 (MQ=255)
                                  gAACTGATCGTTGGTGGTTCAGAGGAAGAGTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAg                                                                                                             <  2:210373/149‑1 (MQ=255)
                                   aaaTGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAg                                                                                                                                                   >  5:32875/1‑111 (MQ=255)
                                   aaaTGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAg                                                                                                                                                   <  6:32875/111‑1 (MQ=255)
                                   aaaTGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGc                                                                                                            >  5:143688/1‑149 (MQ=255)
                                     aTGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTc                                                                                                                <  3:176834/143‑1 (MQ=255)
                                     aTGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTc                                                                                                                >  4:176834/1‑143 (MQ=255)
                                     aTGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTg                                                                                                          <  1:273074/149‑1 (MQ=255)
                                       gATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTa                                                                                                        >  6:145208/1‑149 (MQ=255)
                                        aTCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTac                                                                                                       >  8:21226/1‑149 (MQ=255)
                                        aTCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTac                                                                                                       <  2:6208/149‑1 (MQ=255)
                                        aTCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTa                                                                                                        <  6:143688/148‑1 (MQ=255)
                                         tCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTaca                                                                                                      >  6:223141/1‑149 (MQ=255)
                                            tAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATATCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACAcc                                                                                                    >  5:157110/1‑148 (MQ=255)
                                             aGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACAcc                                                                                                    >  6:104736/1‑147 (MQ=255)
                                              ggtggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGTCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAg                                                                                                 >  6:101622/1‑149 (MQ=255)
                                              ggtggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAg                                                                                                 <  7:21226/149‑1 (MQ=255)
                                                  gtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTAGGCTTGCCGTGCCAGCT‑TCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcg                                                                                              <  1:262092/148‑1 (MQ=255)
                                                      agagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGcggcgg                                                                                          <  5:104736/148‑1 (MQ=255)
                                                      agagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcggg                                                                                         >  8:321227/1‑149 (MQ=255)
                                                      agagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAAAAGAAAAATCAGTCCAGCTGTACACCCAGGcggcgg                                                                                          <  7:51617/148‑1 (MQ=255)
                                                      agagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCAGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcggg                                                                                         >  3:107693/1‑149 (MQ=255)
                                                        agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcgggcg                                                                                       >  8:322671/1‑149 (MQ=255)
                                                         ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACAccc                                                                                                   >  8:159812/1‑136 (MQ=255)
                                                         ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACAccc                                                                                                   <  7:159812/136‑1 (MQ=255)
                                                            aGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACgg                                                                                   <  7:161736/149‑1 (MQ=255)
                                                                  ccACAGGAAAGCGGGATGAATTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGAAAGATGATGTCCAGAAAAATCAGTCCCGCGGTAATCCCAGTCGGCGGGCGACGCCTTCa                                                                              >  2:68668/1‑148 (MQ=255)
                                                                       ggAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCTGGCGACGGCTTCATAGGCt                                                                        <  8:84640/149‑1 (MQ=255)
                                                                           aGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGa                                                                    <  5:85196/149‑1 (MQ=255)
                                                                              gggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCa                                                                 <  7:158833/149‑1 (MQ=255)
                                                                              gggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCa                                                                 <  1:132564/149‑1 (MQ=255)
                                                                              gggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCa                                                                 <  5:16435/149‑1 (MQ=255)
                                                                               ggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAg                                                                <  5:223141/149‑1 (MQ=255)
                                                                               ggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAg                                                                <  8:48297/149‑1 (MQ=255)
                                                                                gATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCCCGGCAAGATGATGAACCGAAAACTCAGTCCAGCTGTACACCCAGGCGGCGCGCGACGGCTTCATAGGCTTCGATCAgg                                                               >  5:296958/1‑149 (MQ=255)
                                                                                gATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAgg                                                               >  1:256953/1‑149 (MQ=255)
                                                                                 aTGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACCCCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGc                                                              >  8:4062/1‑149 (MQ=255)
                                                                                       tGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGa                                                        <  8:300521/149‑1 (MQ=255)
                                                                                         gCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGTCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAg                                                       <  1:68668/148‑1 (MQ=255)
                                                                                           ggaggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGct                                                    >  1:109011/1‑149 (MQ=255)
                                                                                            gaggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGctc                                                   >  3:110107/1‑149 (MQ=255)
                                                                                             aggaATGAGCGGTCACGGAATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCCCAGCAAGATGATGAACAGAAGAATGAGTACAGCTGTACACCCAGGCGCCGGGCGACGGCTTCATAGGCTGCGATCAGGCCACCGAGGctct                                                  >  6:118089/1‑149 (MQ=255)
                                                                                               gaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCCATCAGGCCACCGAGGCTCTgg                                                >  2:97163/1‑149 (MQ=255)
                                                                                                 aTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCg                                              >  8:264411/1‑149 (MQ=255)
                                                                                                   gAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCgg                                             >  5:103467/1‑148 (MQ=255)
                                                                                                       gTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCCGGCGGCGGGCGCCGGCTTCATAGGCATCGATCAGGCCACCGCGGCTCTGGCGGCAACg                                        >  8:152312/1‑149 (MQ=255)
                                                                                                        ttCAATCCATATATTCCGAGATCGGGTGGCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAACATCAGTCCACCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGCTCAGGCCACCGAGGCTCTGGCGGAAACgg                                       <  3:34102/149‑1 (MQ=255)
                                                                                                           aaTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTc                                     <  4:175292/148‑1 (MQ=255)
                                                                                                           aaTGCATATATTCCGAGAGCGGTTGTCTTGCCGGGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCt                                    <  5:205458/149‑1 (MQ=255)
                                                                                                            aTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTCCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCtt                                   >  1:277411/1‑149 (MQ=255)
                                                                                                             ttcATATATTACGAGTGCGGTTGGCTTGCAGTGCCTGCT‑GCAAGGAAAGATGATTTACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGGTTACTAGGCTTCGATCAGGCCACGGATGATCTGGCGGAAACGGTCttt                                  <  4:23781/147‑1 (MQ=255)
                                                                                                             tGCATATATTCCGAGAGCTGTTGTCTTGCCGTGCGAGCT‑GCACGGCAAGATGATGAAGAGAAAAATCAGTCCAGGTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCttt                                  <  4:107693/149‑1 (MQ=255)
                                                                                                             tGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCttt                                  <  7:321227/149‑1 (MQ=255)
                                                                                                                atatatTCCGAGAGCGGTTGTCTTGCCGTGCCGGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTc                               >  8:225034/1‑149 (MQ=255)
                                                                                                                  atatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCttt                                  >  4:82251/1‑144 (MQ=255)
                                                                                                                  atatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCttt                                  <  3:82251/144‑1 (MQ=255)
                                                                                                                    atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCGCTGGCGGAAACGGTCt                                    <  7:279944/140‑1 (MQ=255)
                                                                                                                    atTCCGAGAGCGGTTCTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGACGCTCGGGCGGAAACGGACt                                    >  8:279944/1‑140 (MQ=255)
                                                                                                                     ttCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCAttt                          >  7:153698/1‑149 (MQ=255)
                                                                                                                        cGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTtct                       >  6:263216/1‑149 (MQ=255)
                                                                                                                        cGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTtc                        <  6:191117/148‑1 (MQ=255)
                                                                                                                          agagCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGCTGAACAGAAAAATCAGTCCCGGTGTACCCCCAGGCGGCGGGCGCCGGCTTCATAGGCTTCTCTCAGGCCACCGCGGCTCTGGCGGCAACGGGCTTTGGCCATTTtctc                      >  1:300188/1‑148 (MQ=255)
                                                                                                                          agagCGCTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGCCGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTtctc                      >  8:178494/1‑148 (MQ=255)
                                                                                                                           gagCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCTATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCa                    <  5:145208/149‑1 (MQ=255)
                                                                                                                           gaccGGGTGTCTTGCCGTGCCAGCT‑GCACTGCAAGATGATGAACAGAAAAATCAGTCCCGCTGAACACCCATTCGTCGGGCGACGGCTTCATACGCTTCTATCAGGCCTACGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCa                    <  2:98179/146‑1 (MQ=255)
                                                                                                                            agCGGTTGTCTTGCCGTGCCAGCA‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGCCGGCTTCATAGGCTTCGCTAAGGCCACCGCGGCTCTGGCGGAAACGGTCATTGTCCATTTTCTCCAg                   >  8:235439/1‑149 (MQ=255)
                                                                                                                               catcatCTTGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGt                >  4:202217/6‑149 (MQ=255)
                                                                                                                                gTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGtt               >  3:260816/1‑149 (MQ=255)
                                                                                                                                 ttGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTcc                     >  3:336832/1‑142 (MQ=255)
                                                                                                                                 ttGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTcc                     <  4:336832/142‑1 (MQ=255)
                                                                                                                                    tCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCGTTGTCCATTTTCTCCCGCGTTTCtt           >  4:191231/1‑149 (MQ=255)
                                                                                                                                        gCCGTGCCAGCT‑GCACGGCAAGATGATGAAGAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGAGGGCTTCATAGGCTTCTATCAGGCCACCGAGGCTCTGGCTGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCttt          <  6:157110/146‑1 (MQ=255)
                                                                                                                                        gCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGt        <  7:152312/148‑1 (MQ=255)
                                                                                                                                         ccGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGGCTTTGTCCATTTTCTCCAGCGTTTCTTTGTc       <  7:322671/148‑1 (MQ=255)
                                                                                                                                           gTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCca    <  4:260816/149‑1 (MQ=255)
                                                                                                                                            tGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCca    >  4:264154/1‑148 (MQ=255)
                                                                                                                                              ccAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCcaca  <  4:110107/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
AACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACA  >  NC_000913/2596731‑2597020

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