Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,596,878:1 +G 5.9% intergenic (‑141/+27) 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,8781.G 5.9% 350.3 / 15.5 85intergenic (‑141/+27)ypfJ/purCputative neutral zinc metallopeptidase/phosphoribosylaminoimidazole‑succinocarboxamide synthetase
Reads supporting (aligned to +/- strand):  ref base . (40/40);  new base G (3/2);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-01

CAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAGG  >  NC_000913/2596730‑2597022
                                                                                                                                                     |                                                                                                                                                 
cAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCt                                                                                                                                                    >  1:1128960/1‑149 (MQ=255)
cAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCt                                                                                                                                                    >  2:650075/1‑149 (MQ=255)
 aaCGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTg                                                                                                                                                  <  2:1013239/149‑1 (MQ=255)
    gctttgATACCTCAACAATATATTAATTATAGAAATGCTCGAAGGTGGTTCATAGGAAGATTACCACAGGAAAGCTGGATGAAATGGCGGAGGAATGAGGGTTCAATGCTTATATTCCGAGACCGGTTGTCTTGCCCTGCGAGCTGCAc                                                                                                                                               <  2:516122/144‑1 (MQ=255)
      aTAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGGGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACgg                                                                                                                                             <  2:1411985/149‑1 (MQ=255)
           tACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAga                                                                                                                                        <  2:589777/149‑1 (MQ=255)
             ccTCAACAATATATTCATTATAGAAATGATCGTTGGTGGTTCAGAGGAAGATTACCAAAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTTCGAGGTGCGCTGCAAgatt                                                                                                                                      >  2:785286/1‑148 (MQ=255)
             ccTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgatt                                                                                                                                      >  2:424804/1‑148 (MQ=255)
              cTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgatta                                                                                                                                     >  1:61994/1‑147 (MQ=255)
              cTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgatga                                                                                                                                     <  2:59014/149‑1 (MQ=255)
              cTCAACAATATATTCATAATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAgatga                                                                                                                                     <  2:119567/149‑1 (MQ=255)
               tCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTATCTTGCCGTGCAGCTGGCACGGCAAgacaac                                                                                                                                    <  2:1484600/149‑5 (MQ=255)
                 aacaaTATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGAtgatga                                                                                                                                  >  2:888337/1‑149 (MQ=255)
                 aacaaTATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGAtgatga                                                                                                                                  >  2:1347505/1‑149 (MQ=255)
                  acaaTATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGaa                                                                                                                                 <  1:180753/149‑1 (MQ=255)
                   caaTATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGaaa                                                                                                                                >  1:596179/1‑148 (MQ=255)
                   caaTATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAAc                                                                                                                                >  1:1340876/1‑149 (MQ=255)
                            cATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAAta                                                                                                                       >  2:338413/1‑148 (MQ=255)
                            cATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATc                                                                                                                       <  2:924720/149‑1 (MQ=255)
                             aTTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCa                                                                                                                      >  2:1531765/1‑149 (MQ=255)
                              ttATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAg                                                                                                                     <  1:762820/149‑1 (MQ=255)
                              ttATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAg                                                                                                                     <  2:1492034/149‑1 (MQ=255)
                              ttATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAg                                                                                                                     >  2:1333004/1‑149 (MQ=255)
                               tataGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCATTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCATGATGATGAACAGAAAAATCAGt                                                                                                                    <  1:911024/149‑1 (MQ=255)
                               tataGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGACCAGAAAAATCAGt                                                                                                                    >  2:1042798/1‑149 (MQ=255)
                                ataGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTc                                                                                                                   <  2:1538246/149‑1 (MQ=255)
                                ataGAAATGATCGTAGGTGGTTCAGAGGAAGAATACCACAGGAAAGCGGGATGAAATGGAGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTTCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTc                                                                                                                   >  1:1133071/1‑149 (MQ=255)
                                   gAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATTCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAACATGATGAACAGAAAAATCAGTCCAg                                                                                                                >  2:197904/1‑149 (MQ=255)
                                    aaaTGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGc                                                                                                               >  2:285751/1‑149 (MQ=255)
                                       tGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGACAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGt                                                                                                            <  1:251606/149‑1 (MQ=255)
                                           cGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTacac                                                                                                       >  1:478855/1‑149 (MQ=255)
                                            gTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACAcc                                                                                                      <  1:1212524/149‑1 (MQ=255)
                                             tAGGTGGTTCAGAGGAAGATTACCAAAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACAccc                                                                                                     >  2:1006183/1‑149 (MQ=255)
                                                    ttCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcggc                                                                                              >  1:935936/1‑149 (MQ=255)
                                                         agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcgggcg                                                                                         >  1:533269/1‑149 (MQ=255)
                                                         agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcgggcg                                                                                         >  2:613426/1‑149 (MQ=255)
                                                          ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGa                                                                                        <  1:557969/149‑1 (MQ=255)
                                                          ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGACGTGGCAGAATCACGGCAAAATGGAAAACAAAAAAATAAATCCACCGAGACACCCCAGGcgccggcga                                                                                         >  2:1212197/1‑146 (MQ=255)
                                                           gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGAc                                                                                       >  1:683309/1‑149 (MQ=255)
                                                           gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGAc                                                                                       <  2:1340876/149‑1 (MQ=255)
                                                                   ccACAGGCAAGCGGGAGGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGACTTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCAt                                                                               <  2:95647/149‑1 (MQ=255)
                                                                    cacaGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                                                                              >  2:919545/1‑149 (MQ=255)
                                                                      caGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAgg                                                                            >  1:1282911/1‑149 (MQ=255)
                                                                       aGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGc                                                                           <  1:1449717/149‑1 (MQ=255)
                                                                       aGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGc                                                                           <  2:1128960/149‑1 (MQ=255)
                                                                         gAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCtt                                                                         >  1:1099120/1‑149 (MQ=255)
                                                                         gAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCtt                                                                         <  2:1474582/149‑1 (MQ=255)
                                                                             gagggCTGAAATGGCGGAGGACTGAGCGTTCAATGCATATAATCCGCGAGCGGGTGGCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACCCCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGAt                                                                     <  2:65683/147‑1 (MQ=255)
                                                                             gCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACCCCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGAt                                                                     >  1:649481/1‑149 (MQ=255)
                                                                               gggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTTCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTACAGCTTTAAACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCa                                                                   >  2:504311/1‑149 (MQ=255)
                                                                               gggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCa                                                                   >  2:1049091/1‑149 (MQ=255)
                                                                                  aTGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGc                                                                <  2:935936/149‑1 (MQ=255)
                                                                                      aaTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAcc                                                            <  2:1282911/149‑1 (MQ=255)
                                                                                                 aaTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGc                                                 <  2:61994/149‑1 (MQ=255)
                                                                                                  aTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTTCCGTGCCAGCTGCCCGGCAAGATGATGAACAGAAAAATCAGTACAGCTGTACACCCCCGCGGCGGGCGCCGGGTTCATAGGCTTCGATCAGGCCACCCACGCTCTGGCg                                                >  2:1319193/1‑149 (MQ=255)
                                                                                                  aTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCg                                                >  2:300393/1‑149 (MQ=255)
                                                                                                  aTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCg                                                >  1:1595093/1‑149 (MQ=255)
                                                                                                        gTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACg                                          <  2:1043359/149‑1 (MQ=255)
                                                                                                        gTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACg                                          >  2:1042873/1‑149 (MQ=255)
                                                                                                           aaaTGCATATATTGCGAGAGCGGTTGTCTTGCCGTGCCAGCTGGACGGAAAGATGATGAACAGAAAAATCAGTCCATCTGTACACCCAGGGGGCGGGCGACGGCTTCATAGTCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTc                                       <  2:1508161/148‑1 (MQ=255)
                                                                                                               gCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTg                                   <  1:1515940/149‑1 (MQ=255)
                                                                                                                cATATATTCCGAGAGCGGTTGTCTTGCCTTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGt                                  <  1:650075/149‑1 (MQ=255)
                                                                                                                cATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACTGCAAGATGATGAACAGAAAATTCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGTCTCTGGCGGAAACGGTCTTTGt                                  <  2:608255/149‑1 (MQ=255)
                                                                                                                   atatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCa                               >  1:924286/1‑149 (MQ=255)
                                                                                                                   atatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCGTCATATGCTTCGATCAGGCCCCCGAGGCTCTGGCGGAAACGGTCTTTGTCCa                               >  2:480417/1‑149 (MQ=255)
                                                                                                                    aatTCCGAGAGCGGTGGTCTTGCCGTGCGAGCTGCACGGCAAGATGATGCACAGAAAAATCCGTCATGCTGTACTACCAGGCGGCTGGCGACGGCTTCATAGGCTTCGATCAGTCCACCGAGGCGCTGGCGGAAACTGTCTTTGTCCAt                              <  1:237657/148‑1 (MQ=255)
                                                                                                                     atTCCGAGAGCGGTTGTCTTGCCATGCCAGCAGCACGGCAAGATGATGAACAGAAAAATCAGTCCAACTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTACGATCAAGCCACCGAGGCTCTGGCGGAAACCGTCTTAGTCCaat                             >  2:448949/1‑147 (MQ=255)
                                                                                                                      ttCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGGGACGGCGTCATAGGCTTGGGTCAGGCCACCGAGGCTCTGGGGGAAACGGTCTTTGTCCAttt                            >  1:609694/1‑149 (MQ=255)
                                                                                                                        ccGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTtc                          >  1:198274/1‑149 (MQ=255)
                                                                                                                        ccGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTtc                          <  2:478855/149‑1 (MQ=255)
                                                                                                                        ccGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGGCTTTGTCCATTTtc                          <  1:285751/149‑1 (MQ=255)
                                                                                                                         cGAGAGCGGTTGTCTTGCCGTGCCAGCTGCCCGGGAAGATGATGCACAGAAAAATCAGTCCAGCTGTTCACCCAGGCGGCGGGCGCCGGGGTCATAGGGGTCGAGCAGGCCCCCGCGGCGCTGGCGGAAACGGTCTTTGTGCATTTtct                         >  1:647373/1‑149 (MQ=255)
                                                                                                                         cGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATAAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGACGCTCTGGCGGAAACGGTCTTTGTACAttttat                         >  2:1376112/1‑147 (MQ=255)
                                                                                                                             cgCGGTTGTCTTGCCGTGCCAGCGGCACGGCAAGATGATGCACAGAAAAATCAGTCCCGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAg                     <  2:1161405/148‑1 (MQ=255)
                                                                                                                             agCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAg                     <  2:1167752/149‑1 (MQ=255)
                                                                                                                              gCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGc                    <  1:1347505/149‑1 (MQ=255)
                                                                                                                                ggTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGt                  <  1:47273/149‑1 (MQ=255)
                                                                                                                                 gTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGtt                 <  1:1319193/149‑1 (MQ=255)
                                                                                                                                  tcatCTTGCCGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGttt                >  2:322493/4‑149 (MQ=255)
                                                                                                                                    gTCTTGCCGGTCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTAGACCCATGCGGCGGGCGACGGCTTCTTAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCt              <  1:1212197/149‑1 (MQ=255)
                                                                                                                                        tGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGt          >  2:1141801/1‑149 (MQ=255)
                                                                                                                                        tGCCGGGCCAGCTGCAAGGCAAGATGATGAACAGAAAAATCAGTCCCGCTGTCCACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGt          <  1:1006183/149‑1 (MQ=255)
                                                                                                                                         gCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGAGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTc         <  2:830284/149‑1 (MQ=255)
                                                                                                                                          ccGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCCCCGCGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTcc        >  1:461559/1‑149 (MQ=255)
                                                                                                                                          ccGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGaaa                                            >  1:1412899/1‑113 (MQ=255)
                                                                                                                                          ccGTGCAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGaaa                                            <  2:1412899/113‑1 (MQ=255)
                                                                                                                                               ccAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAg   >  2:679726/1‑149 (MQ=255)
                                                                                                                                                cAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAgg  <  1:1438202/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCTGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAGG  >  NC_000913/2596730‑2597022

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: