Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,948,964 C→T A336T (GCC→ACC)  yifB ← magnesium chelatase family protein and putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,948,9640CT100.0% 520.4 / NA 136A336T (GCC→ACC) yifBmagnesium chelatase family protein and putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (56/80);  total (56/80)

Only 100 of 136 total aligned reads displayed.

ATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGG  >  NC_000913/3948816‑3949110
                                                                                                                                                    |                                                                                                                                                  
     cAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAg                                                                                                                                             <  2:336906/151‑1 (MQ=255)
       gTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGt                                                                                                                                           <  4:1197617/151‑1 (MQ=255)
         gagagaTCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTa                                                                                                                                         <  1:1065506/151‑1 (MQ=255)
          agagaTCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTAt                                                                                                                                        >  1:1044173/1‑151 (MQ=255)
            agaTCGAAGCGGTCGAGAAAGGGCCCCGAGCGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTAtt                                                                                                                                       <  4:1500207/150‑1 (MQ=255)
              aTCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTAttttt                                                                                                                                    <  1:338864/151‑1 (MQ=255)
                cGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTg                                                                                                                                   <  4:650569/150‑1 (MQ=255)
                cGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTg                                                                                                                                   <  1:1431627/150‑1 (MQ=255)
                cGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTg                                                                                                                                   <  1:803517/150‑1 (MQ=255)
                 gAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGtt                                                                                                                                          >  3:337804/1‑142 (MQ=255)
                 gAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCt                                                                                                                                 >  1:1906346/1‑151 (MQ=255)
                      gggcGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTg                                                                                                                                   <  3:388160/141‑1 (MQ=255)
                       gTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGtgtg                                                                                                                           <  4:950075/151‑1 (MQ=255)
                        tCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGc                                                                                                                          >  3:166305/1‑151 (MQ=255)
                         cGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGACTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCg                                                                                                                         <  2:326354/151‑1 (MQ=255)
                          gagaAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGt                                                                                                                        >  4:1210749/1‑151 (MQ=255)
                            gaAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGa                                                                                                                      <  1:692558/151‑1 (MQ=255)
                              aaGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGaa                                                                                                                     >  4:642702/1‑150 (MQ=255)
                               aGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGaa                                                                                                                     >  1:615100/1‑149 (MQ=255)
                                    cccGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATgg                                                                                                               >  2:513697/1‑150 (MQ=255)
                                    cccGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGa                                                                                                              <  4:1611922/151‑1 (MQ=255)
                                    cccGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGa                                                                                                              <  4:1741298/151‑1 (MQ=255)
                                         gagCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGc                                                                                                         >  1:1279234/1‑151 (MQ=255)
                                         gagCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGc                                                                                                         >  1:1480174/1‑151 (MQ=255)
                                         gagCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGc                                                                                                         <  4:154471/151‑1 (MQ=255)
                                           gCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGccc                                                                                                       >  1:755031/1‑151 (MQ=255)
                                           gCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGccc                                                                                                       >  1:902171/1‑151 (MQ=255)
                                               gTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGAt                                                                                                   <  4:1882681/151‑1 (MQ=255)
                                                ttGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGAtt                                                                                                  >  1:605608/1‑151 (MQ=255)
                                                    gaTAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAAt                                                                                              >  4:1794178/1‑151 (MQ=255)
                                                      tAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCg                                                                                            >  1:561376/1‑151 (MQ=255)
                                                      tAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCg                                                                                            >  2:516992/1‑151 (MQ=255)
                                                      tAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCg                                                                                            >  4:1861097/1‑151 (MQ=255)
                                                          gTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGctc                                                                                        >  1:1798785/1‑151 (MQ=255)
                                                           tAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGctct                                                                                       >  4:1672217/1‑151 (MQ=255)
                                                            aaTGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGctct                                                                                       <  1:1302540/150‑1 (MQ=255)
                                                             aTGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCg                                                                                     <  4:503299/151‑1 (MQ=255)
                                                              tGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGc                                                                                    >  3:306889/1‑151 (MQ=255)
                                                               gTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCa                                                                                   <  3:375900/151‑1 (MQ=255)
                                                                tCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCa                                                                                   >  3:545012/1‑150 (MQ=255)
                                                                 cTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCa                                                                                   >  1:863608/1‑149 (MQ=255)
                                                                   gTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGc                                                                               >  1:494974/1‑151 (MQ=255)
                                                                     tCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCAt                                                                             >  4:827859/1‑151 (MQ=255)
                                                                     tCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCAt                                                                             >  1:100793/1‑151 (MQ=255)
                                                                       tGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATcc                                                                           >  4:591901/1‑151 (MQ=255)
                                                                        ggCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATc                                                                            <  1:38071/149‑1 (MQ=255)
                                                                             gCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAgtgt                                                                      >  3:365355/1‑150 (MQ=255)
                                                                             gCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAgtgt                                                                      <  1:1207071/150‑1 (MQ=255)
                                                                             gCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTa                                                                     <  4:1433009/151‑1 (MQ=255)
                                                                                 cGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGcc                                                                 <  4:1561991/151‑1 (MQ=255)
                                                                                  ggTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCg                                                                >  1:1104827/1‑151 (MQ=255)
                                                                                   gTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGt                                                               >  1:343749/1‑151 (MQ=255)
                                                                                    ttATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGtt                                                              >  4:1772269/1‑151 (MQ=255)
                                                                                     tATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGtt                                                              <  3:340878/150‑1 (MQ=255)
                                                                                      aTGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTc                                                             >  1:311733/1‑150 (MQ=255)
                                                                                       aggTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCa                                                            >  3:213198/2‑150 (MQ=255)
                                                                                        ggTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCaa                                                           <  3:179889/150‑1 (MQ=255)
                                                                                        ggTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCaa                                                           <  4:1093873/150‑1 (MQ=255)
                                                                                         gTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAAt                                                         <  4:1863949/151‑1 (MQ=255)
                                                                                          tttCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAAt                                                         <  1:730435/150‑1 (MQ=255)
                                                                                              ccTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCATTGTACGCCGTTCAAATTCAgg                                                    <  1:301621/151‑1 (MQ=255)
                                                                                              ccTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCa                                                      <  4:365568/149‑1 (MQ=255)
                                                                                              ccTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAgg                                                    <  2:124816/151‑1 (MQ=255)
                                                                                              ccTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAgg                                                    <  1:70551/151‑1 (MQ=255)
                                                                                              ccTGAGAATGTCCGGTAGGGCTTGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAgg                                                    <  3:534109/151‑1 (MQ=255)
                                                                                                tGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGt                                                   >  2:330278/1‑150 (MQ=255)
                                                                                                 gATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTa                                                  <  1:589870/150‑1 (MQ=255)
                                                                                                       gTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCtcatc                                           >  1:1374351/1‑151 (MQ=255)
                                                                                                         ccGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCaa                                         <  4:520626/151‑1 (MQ=255)
                                                                                                         ccGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCaa                                         <  2:84721/151‑1 (MQ=255)
                                                                                                           ggTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGa                                       <  4:290058/151‑1 (MQ=255)
                                                                                                            gTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGaa                                      <  1:1706724/151‑1 (MQ=255)
                                                                                                             tAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGaa                                      >  4:21940/1‑150 (MQ=255)
                                                                                                               gggCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGaaaaa                                   <  1:1020448/151‑1 (MQ=255)
                                                                                                               gggCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGaaaaa                                   <  2:67356/151‑1 (MQ=255)
                                                                                                                ggCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAg                                  <  4:1044622/151‑1 (MQ=255)
                                                                                                                 gCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGc                                 <  2:167497/151‑1 (MQ=255)
                                                                                                                    gggATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCa                                <  4:1334423/149‑1 (MQ=255)
                                                                                                                      gATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGTACGcc                            <  4:764707/151‑1 (MQ=255)
                                                                                                                      gATTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGcc                            <  4:80926/151‑1 (MQ=255)
                                                                                                                       aTTCATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCg                           <  4:291208/151‑1 (MQ=255)
                                                                                                                          cATCGCCGCAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTa                        <  2:389285/151‑1 (MQ=255)
                                                                                                                             cgccgcAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATgc                     >  4:462248/1‑151 (MQ=255)
                                                                                                                                 gcAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCa                 <  2:350750/151‑1 (MQ=255)
                                                                                                                                 gcAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCa                 <  1:1210234/151‑1 (MQ=255)
                                                                                                                                 gcAACAAGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCa                 <  1:855105/151‑1 (MQ=255)
                                                                                                                                   aacaaGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAg                <  4:935069/150‑1 (MQ=255)
                                                                                                                                   aacaaGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAg                <  4:1838553/150‑1 (MQ=255)
                                                                                                                                   aacaaGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGc               <  2:542808/151‑1 (MQ=255)
                                                                                                                                     caaGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCg              >  4:679226/1‑150 (MQ=255)
                                                                                                                                      aaGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGaa            <  1:1771577/151‑1 (MQ=255)
                                                                                                                                      aaGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGa             >  1:298972/1‑150 (MQ=255)
                                                                                                                                       aGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCGTCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGaa            >  1:222802/1‑150 (MQ=255)
                                                                                                                                       aGCTGGAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGaa            >  4:125017/1‑150 (MQ=255)
                                                                                                                                           ggAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAAtt         <  4:5224/149‑1 (MQ=255)
                                                                                                                                           ggAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAAtt         <  2:350294/149‑1 (MQ=255)
                                                                                                                                            gAAACGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTc       >  4:1106624/1‑150 (MQ=255)
                                                                                                                                               aCGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCg   <  1:654090/151‑1 (MQ=255)
                                                                                                                                               aCGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCg   >  4:1918128/1‑151 (MQ=255)
                                                                                                                                                cGGGTTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCgg  <  4:215707/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCGCAACAAGCTGGAAACGGGCTGGATAGGTTATTTTTGCTCGTGTGCGTGAAAGATGGATCTGCCCGGATTCAATCGGCTCTCGCAAGGCATCCAGTGTACGCCGTTCAAATTCAGGTAGCTCATCAAGAAAAAGCACGCCGTTATGCGCCAGCGAAATTTCACCGG  >  NC_000913/3948816‑3949110

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