Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,448,543 +C intergenic (+48/‑90) 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,5421.C99.2% 719.9 / ‑1.6 246intergenic (+47/‑91)ytfQ/ytfRputative LACI‑type transcriptional regulator/YtfR
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (143/101);  total (144/102)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 246 total aligned reads displayed.

AAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGGCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCG  >  NC_000913/4448461‑4448623
                                                                                     |                                                                                 
aagaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCgg                                                                          >  2:2707693/1‑93 (MQ=37)
aagaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCgg                                                                          >  1:5184553/1‑93 (MQ=37)
aagaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCgg                                                                          <  1:4296606/93‑1 (MQ=37)
aagaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCgg                                                                          <  2:2781884/93‑1 (MQ=37)
  gaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCg                                                                        >  2:5019078/1‑93 (MQ=37)
   aaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGa                                                                       >  1:5952261/1‑93 (MQ=37)
   aaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGa                                                                       >  1:5622992/1‑93 (MQ=37)
   aaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGa                                                                       <  2:2727869/93‑1 (MQ=37)
     ttAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATg                                                                     <  2:370847/93‑1 (MQ=37)
     ttAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATg                                                                     >  2:1515780/1‑93 (MQ=37)
         agaagaagaAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGc                                                                 <  2:1580558/93‑1 (MQ=37)
         agaagaagaAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGc                                                                 <  1:7123884/93‑1 (MQ=37)
           aagaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCAc                                                               <  1:4012173/93‑1 (MQ=37)
            agaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACt                                                              >  1:778534/1‑93 (MQ=37)
             gaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTg                                                             >  2:5325257/1‑93 (MQ=37)
               agaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGaa                                                           >  1:2645404/1‑93 (MQ=37)
                gaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAg                                                          >  1:5749213/1‑93 (MQ=37)
                 aaaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGc                                                         <  1:4044776/93‑1 (MQ=37)
                 aaaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGc                                                         <  1:5722209/93‑1 (MQ=37)
                  aaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGcc                                                        >  2:1955109/1‑93 (MQ=37)
                  aaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGcc                                                        <  2:6008745/93‑1 (MQ=37)
                  aaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGcc                                                        <  1:1989271/93‑1 (MQ=37)
                   aaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCt                                                       >  1:6446466/1‑93 (MQ=37)
                   aaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCt                                                       <  2:3959610/93‑1 (MQ=37)
                    aaTATGGGGTATTGAGGTTTGCTATGCCTGATGCCGATTCGTGGGCGGGAAAGGGCCCT‑‑CGCGCCCCACCCGGCGATGGTGCACTGAAGCCTg                                                      <  1:2520544/93‑1 (MQ=255)
                    aaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTg                                                      >  2:3028803/1‑93 (MQ=37)
                    aaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTg                                                      <  1:2799070/93‑1 (MQ=37)
                    aaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTg                                                      >  2:1215862/1‑93 (MQ=37)
                     atatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGa                                                     <  1:1200016/93‑1 (MQ=37)
                     atatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGa                                                     <  2:3385002/93‑1 (MQ=37)
                     atatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGa                                                     >  1:5325376/1‑93 (MQ=37)
                      tatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGAt                                                    >  2:2885801/1‑93 (MQ=37)
                      tatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGAt                                                    <  2:4289226/93‑1 (MQ=37)
                      tatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGAt                                                    >  1:4091542/1‑93 (MQ=25)
                       atGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAGGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATg                                                   <  2:323074/93‑1 (MQ=37)
                        tGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGc                                                  >  1:6728920/1‑93 (MQ=25)
                        tGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGc                                                  >  1:7277802/1‑93 (MQ=25)
                        tGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGc                                                  >  1:4555668/1‑93 (MQ=25)
                          gggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGa                                                <  1:1139022/93‑1 (MQ=25)
                          gggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCCCC‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGa                                                <  1:2940679/93‑1 (MQ=25)
                           ggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGAc                                               <  1:2881673/93‑1 (MQ=25)
                            gTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACg                                              >  1:6489060/1‑93 (MQ=25)
                             tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGc                                             >  2:4598861/1‑93 (MQ=25)
                             tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGc                                             >  1:1877079/1‑93 (MQ=25)
                             tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGc                                             >  2:2588681/1‑93 (MQ=25)
                             tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGc                                             >  1:7210117/1‑93 (MQ=25)
                              aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCt                                            <  2:5733596/93‑1 (MQ=25)
                              aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCt                                            >  1:2242021/1‑93 (MQ=25)
                                tGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTa                                          <  2:2060734/93‑1 (MQ=25)
                                 gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAc                                         <  1:2715246/93‑1 (MQ=25)
                                  aGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAcc                                        >  1:1687839/1‑93 (MQ=25)
                                  aGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAcc                                        <  1:6565452/93‑1 (MQ=25)
                                   gggTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcg                                       >  1:3063412/1‑93 (MQ=25)
                                   gggTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcg                                       >  2:5859940/1‑93 (MQ=25)
                                     gTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCGCTGAAGCCTGATGCGACGCTTACcgcg                                     >  2:2392762/1‑93 (MQ=25)
                                     gTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgcg                                     >  1:4391015/1‑93 (MQ=25)
                                     gTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgcg                                     >  2:2314919/1‑93 (MQ=25)
                                        gCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCt                                  <  1:3429580/93‑1 (MQ=25)
                                          tATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCCCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTa                                >  1:5086182/1‑93 (MQ=21)
                                             gCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCa                             >  2:4892112/1‑93 (MQ=21)
                                             gCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCa                             >  2:4903233/1‑93 (MQ=21)
                                              ccTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCAt                            <  2:5692630/93‑1 (MQ=25)
                                              ccTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCAt                            <  2:4214896/93‑1 (MQ=21)
                                               cTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCCCTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATg                           <  1:7447138/93‑1 (MQ=21)
                                               cTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATg                           >  1:3604922/1‑93 (MQ=21)
                                                tGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGc                          >  2:1329458/1‑93 (MQ=21)
                                                 gATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGcc                         <  1:3926248/93‑1 (MQ=21)
                                                 gATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGcc                         <  2:4737629/93‑1 (MQ=21)
                                                 gATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGcc                         <  2:5188721/93‑1 (MQ=21)
                                                 gATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGcc                         >  2:4624288/1‑93 (MQ=21)
                                                   tGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTa                       >  1:1112484/1‑93 (MQ=21)
                                                   tGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTa                       >  1:5373845/1‑93 (MQ=21)
                                                    gCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTAc                      >  2:1124695/1‑93 (MQ=21)
                                                     ccGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCAGCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACt                     <  2:3492978/93‑1 (MQ=18)
                                                      cGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTg                    >  2:6306460/1‑93 (MQ=21)
                                                      cGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTg                    <  2:6179814/93‑1 (MQ=21)
                                                       gATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTgg                   <  1:3513714/93‑1 (MQ=21)
                                                        aTTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTggg                  >  1:5746738/1‑93 (MQ=21)
                                                         ttCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGa                 >  2:6837837/1‑93 (MQ=21)
                                                         ttCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGa                 <  2:1031376/93‑1 (MQ=21)
                                                         ttCGGAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGa                 <  1:2588004/93‑1 (MQ=18)
                                                          tCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAg                <  2:7432058/93‑1 (MQ=21)
                                                            gTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCa              <  2:7244894/93‑1 (MQ=21)
                                                            gTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCa              <  1:3192565/93‑1 (MQ=21)
                                                            gTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCa              >  1:380861/1‑93 (MQ=21)
                                                             tAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCAc             >  2:1285063/1‑93 (MQ=21)
                                                               ggCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGc           <  1:1489832/93‑1 (MQ=21)
                                                               ggCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGc           >  1:6057435/1‑93 (MQ=21)
                                                                gCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCt          <  2:6991389/93‑1 (MQ=21)
                                                                gCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCt          >  1:4579200/1‑93 (MQ=21)
                                                                  cGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCttt        >  2:6339686/1‑93 (MQ=21)
                                                                   ggATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTa       <  2:6560443/93‑1 (MQ=21)
                                                                    gATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTac      <  2:5131447/93‑1 (MQ=25)
                                                                    gATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTac      >  1:4698327/1‑93 (MQ=21)
                                                                    gATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTac      >  2:3768332/1‑93 (MQ=21)
                                                                     aTAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTaca     <  1:485153/93‑1 (MQ=25)
                                                                      tAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTacac    <  2:1715684/93‑1 (MQ=25)
                                                                      tAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTacac    >  2:5318358/1‑93 (MQ=25)
                                                                       aaGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACAcc   <  1:5695083/93‑1 (MQ=25)
                                                                        aGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCg  <  2:6786168/93‑1 (MQ=25)
                                                                                     |                                                                                 
AAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGGCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCG  >  NC_000913/4448461‑4448623

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