Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,448,541 +C intergenic (+46/‑92) ytfQ → / → ytfR putative LACI‑type transcriptional regulator/YtfR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,448,5411.C99.5% 593.0 / ‑6.9 197intergenic (+46/‑92)ytfQ/ytfRputative LACI‑type transcriptional regulator/YtfR
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base C (131/65);  minor base A (1/0);  total (132/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 219 total aligned reads displayed.

ATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCG‑‑GCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGG  >  NC_000913/4448449‑4448625
                                                                                             |                                                                                      
aTACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCg                                                                                         >  2:1964592/1‑93 (MQ=255)
 tACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTGGGCCGGATAAGGCGCTCgc                                                                                        <  1:6493841/93‑2 (MQ=255)
 tACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCgc                                                                                        >  1:1067015/1‑92 (MQ=255)
            aagaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCgg                                                                            >  1:4654168/1‑93 (MQ=37)
               aaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGa                                                                         <  2:1366139/93‑1 (MQ=37)
               aaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGa                                                                         <  1:2993777/93‑1 (MQ=37)
               aaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGa                                                                         >  2:4818928/1‑93 (MQ=37)
                aTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGAt                                                                        <  2:5792552/93‑1 (MQ=37)
                  tAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATgg                                                                      <  2:2869261/93‑1 (MQ=37)
                    gagaAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTg                                                                    >  2:5432028/1‑93 (MQ=37)
                      gaagaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCa                                                                  >  1:5728057/1‑93 (MQ=37)
                      gaagaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCa                                                                  >  2:5574189/1‑93 (MQ=37)
                       aagaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCAc                                                                 >  1:111989/1‑93 (MQ=37)
                       aagaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCAc                                                                 >  2:5927447/1‑93 (MQ=37)
                        agaagaaAAATATGGGGTATTGAGGGTTTCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACt                                                                <  2:4455814/93‑1 (MQ=25)
                        agaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACt                                                                >  1:3230136/1‑93 (MQ=37)
                         gaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTg                                                               >  1:2700677/1‑93 (MQ=37)
                         gaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTg                                                               <  1:5262300/93‑1 (MQ=37)
                         gaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTg                                                               >  1:1876506/1‑93 (MQ=37)
                           agaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGaa                                                             >  1:4005472/1‑93 (MQ=37)
                            gaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAg                                                            <  2:5325908/93‑1 (MQ=37)
                             aaaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGc                                                           <  2:3579652/93‑1 (MQ=37)
                              aaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGcc                                                          <  2:4737876/93‑1 (MQ=37)
                              aaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGcc                                                          <  2:6229817/93‑1 (MQ=37)
                                aaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTg                                                        <  1:5043128/93‑1 (MQ=37)
                                 atatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGAGCCGCATCCGGCGATGGTGCACTGAAGCCTGa                                                       >  2:2788784/1‑93 (MQ=37)
                                   atGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATg                                                     >  1:4388087/1‑93 (MQ=25)
                                    tGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGc                                                    <  2:6362254/93‑1 (MQ=25)
                                      gggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGa                                                  >  2:1893232/1‑93 (MQ=25)
                                      gggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGa                                                  >  2:3468006/1‑93 (MQ=25)
                                       ggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGAc                                                 >  2:3765228/1‑93 (MQ=25)
                                        gTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGAATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACg                                                <  2:2834810/93‑1 (MQ=25)
                                         tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGc                                               <  1:5956459/93‑1 (MQ=25)
                                         tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGc                                               >  1:6080207/1‑93 (MQ=25)
                                          aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCt                                              >  1:954150/1‑93 (MQ=25)
                                           ttGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCtt                                             >  1:5859158/1‑93 (MQ=25)
                                            tGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTa                                            <  2:260184/93‑1 (MQ=25)
                                            tGAGGGTTGCTATGCCTGAGGCCGATTCGTAGGCGGGATAAGGCGCTCGCCCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTa                                            <  1:418916/93‑1 (MQ=25)
                                            tGAGGGTGGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTa                                            <  1:431686/93‑1 (MQ=25)
                                             gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAc                                           <  1:2482780/93‑1 (MQ=25)
                                             gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAc                                           >  1:1190281/1‑93 (MQ=25)
                                             gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAc                                           <  1:6589562/93‑1 (MQ=25)
                                             gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGCGATGGTGCACTGAAGCCTGATGCGACGCTTAcc                                          <  2:4608921/93‑1 (MQ=25)
                                               gggTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcg                                         >  1:4715986/1‑93 (MQ=25)
                                               gggTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcg                                         >  2:831137/1‑93 (MQ=25)
                                                ggTTGCTATGCCTGATGCCGATCCGTAGGCCGGATAGGGCGCTCGCCCCGCACCGGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgc                                        <  1:5624711/93‑1 (MQ=25)
                                                 gTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgcg                                       <  2:5089025/93‑1 (MQ=25)
                                                  ttGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGt                                      >  1:538391/1‑93 (MQ=25)
                                                  ttGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGt                                      <  2:4015615/93‑1 (MQ=25)
                                                    tcTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCt                                    <  2:4372954/92‑1 (MQ=25)
                                                    gCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCt                                    >  2:732511/1‑93 (MQ=25)
                                                    gCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCt                                    >  2:944746/1‑93 (MQ=25)
                                                    gCTATGCCTGATGCCGATTCGTAGGCAGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCt                                    <  1:2596472/93‑1 (MQ=25)
                                                     cTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCtt                                   <  2:6278625/93‑1 (MQ=25)
                                                     cTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCtt                                   <  1:5442475/93‑1 (MQ=25)
                                                     cTAGGCCTGATGCCGATTCGTGGGCCGAATAGGGCGCTCGCGCCGCATCCGGCGATGGTCCACTGAAGCCTGATGCGACGCTTACCGCGTCtt                                   <  1:1477784/93‑1 (MQ=21)
                                                      tATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTa                                  >  1:4931723/1‑93 (MQ=21)
                                                      tATGCCTGATGCCGATTCGGAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTa                                  <  2:389789/93‑1 (MQ=21)
                                                       aTGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTAt                                 >  1:279758/1‑93 (MQ=25)
                                                       aTGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTAt                                 >  1:1672227/1‑93 (MQ=21)
                                                       aTGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTAt                                 >  2:2358884/1‑93 (MQ=25)
                                                       aTGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTAt                                 <  2:6016269/93‑1 (MQ=25)
                                                        tGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATc                                >  2:1114309/1‑93 (MQ=21)
                                                        tGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATc                                >  2:1231260/1‑93 (MQ=21)
                                                         gCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCa                               <  2:6073233/93‑1 (MQ=21)
                                                         gCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCa                               >  1:6281152/1‑93 (MQ=21)
                                                          ccTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCAt                              >  1:1545220/1‑93 (MQ=21)
                                                          ccTGATGCCGATTCGTAGGCCGCATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCAt                              <  2:6151006/93‑1 (MQ=21)
                                                          ccTGATGCCGATTCGTAGCCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCAt                              <  1:307313/93‑1 (MQ=21)
                                                           cTGATGCCTATGCGTAGGCGGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATg                             <  1:1649289/93‑1 (MQ=18)
                                                           cTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATg                             >  1:3782990/1‑93 (MQ=21)
                                                            tGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGc                            >  2:4269976/1‑93 (MQ=21)
                                                            tGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGc                            >  2:1828001/1‑93 (MQ=21)
                                                             gATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGcc                           >  1:2032337/1‑93 (MQ=21)
                                                               tGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTa                         <  2:3299856/93‑1 (MQ=21)
                                                               tGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTa                         <  2:5659878/93‑1 (MQ=21)
                                                               tGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTa                         <  2:2189393/93‑1 (MQ=21)
                                                                 cggATTCGTAGGCCGGATAAGGCGCCCGCGCCGCATCGGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACt                       <  1:4535703/91‑1 (MQ=18)
                                                                  cGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTg                      >  1:5774863/1‑93 (MQ=21)
                                                                     ttCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGa                   >  2:2955456/1‑93 (MQ=21)
                                                                      tCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAg                  >  1:1269727/1‑93 (MQ=21)
                                                                       cGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGc                 >  1:2163547/1‑93 (MQ=21)
                                                                        gTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCa                >  1:2345975/1‑93 (MQ=21)
                                                                        gTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCa                >  1:4473730/1‑93 (MQ=21)
                                                                         tAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCAc               >  1:3370489/1‑93 (MQ=21)
                                                                           ggCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGc             >  2:1190899/1‑93 (MQ=21)
                                                                            gCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGGCGCTTACCGCGTCTTATCCTGGCTACTGGGATCACGCt            >  2:3146571/1‑93 (MQ=18)
                                                                            gCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCt            <  1:2884933/93‑1 (MQ=21)
                                                                               ggATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTa         >  1:4081847/1‑93 (MQ=21)
                                                                                gATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTac        <  2:6503851/93‑1 (MQ=25)
                                                                                gATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTac        >  2:3139372/1‑93 (MQ=21)
                                                                                gATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTac        <  2:879395/93‑1 (MQ=25)
                                                                                gATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGCCTTATCATGCCTACTGGGAGCACGCtttgc        >  2:2312996/1‑91 (MQ=21)
                                                                                  tAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTacac      >  1:3674706/1‑93 (MQ=25)
                                                                                  tAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTacac      <  1:2289198/93‑1 (MQ=25)
                                                                                  tAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTacac      <  2:5535684/93‑1 (MQ=25)
                                                                                  tAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTacac      <  2:1630599/93‑1 (MQ=25)
                                                                                    aGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCg    >  2:6631674/1‑93 (MQ=25)
                                                                                    aGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCg    >  1:5153345/1‑93 (MQ=25)
                                                                                    aGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCg    >  2:1774848/1‑93 (MQ=25)
                                                                                             |                                                                                      
ATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCG‑‑GCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGG  >  NC_000913/4448449‑4448625

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: