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.2% 711.1 / ‑1.9 246intergenic (+46/‑92)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.38e-01

Only 100 of 257 total aligned reads displayed.

TACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGGCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCG  >  NC_000913/4448450‑4448623
                                                                                                |                                                                                   
tACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑Cgc                                                                                     >  1:6722404/1‑92 (MQ=255)
tACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑Cgc                                                                                     >  1:3245431/1‑92 (MQ=255)
tACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑Cgc                                                                                     >  2:5395026/1‑92 (MQ=255)
 aCCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑Cgcg                                                                                    >  1:4391703/1‑91 (MQ=255)
  ccGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑Cgcgc                                                                                   <  2:3073286/93‑4 (MQ=255)
  ccGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑Cgcgc                                                                                   <  2:4089614/93‑4 (MQ=255)
   cGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑Cgcgcc                                                                                  <  2:332718/93‑5 (MQ=255)
           aagaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCgg                                                                          >  2:2707693/1‑93 (MQ=37)
            agaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGc                                                                         <  2:1318646/93‑1 (MQ=37)
            agaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGc                                                                         >  1:4348826/1‑93 (MQ=37)
             gaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCg                                                                        >  2:6063031/1‑93 (MQ=37)
             gaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCg                                                                        >  2:5019078/1‑93 (MQ=37)
              aaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGa                                                                       >  1:5622992/1‑93 (MQ=37)
              aaTTAGAGAAGAAGAAAAAGATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGa                                                                       <  2:3423477/93‑1 (MQ=25)
               aTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGAt                                                                      >  1:4928048/1‑93 (MQ=37)
               aTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGAt                                                                      >  1:5207220/1‑93 (MQ=37)
                ttAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATg                                                                     >  2:1515780/1‑93 (MQ=37)
                 tAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATgg                                                                    >  2:6259482/1‑93 (MQ=37)
                    agaagaagaAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGc                                                                 <  2:5721487/93‑1 (MQ=37)
                     gaagaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCa                                                                >  1:2526579/1‑93 (MQ=37)
                      aagaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCAc                                                               <  1:4012173/93‑1 (MQ=37)
                       agaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACt                                                              >  2:2920139/1‑93 (MQ=37)
                        gaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTg                                                             >  1:2955843/1‑93 (MQ=37)
                          agaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGaa                                                           >  2:4692212/1‑93 (MQ=37)
                           gaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAg                                                          >  1:1822832/1‑93 (MQ=37)
                             aaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGcc                                                        >  2:1955109/1‑93 (MQ=37)
                                atatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGa                                                     >  1:5325376/1‑93 (MQ=37)
                                 tatGGGGTATTGGGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGAt                                                    <  1:6473544/93‑1 (MQ=25)
                                 tatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGAt                                                    >  1:4091542/1‑93 (MQ=25)
                                 tatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGAt                                                    >  2:5913782/1‑93 (MQ=37)
                                 tatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGAt                                                    <  2:4385149/93‑1 (MQ=25)
                                 tatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGAt                                                    >  2:2840495/1‑93 (MQ=37)
                                 tatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGAt                                                    >  1:2612497/1‑93 (MQ=37)
                                  atGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAGG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATg                                                   <  2:323074/93‑1 (MQ=37)
                                   tGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGc                                                  >  1:4555668/1‑93 (MQ=25)
                                     gggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGa                                                <  2:1548424/93‑1 (MQ=25)
                                     gggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGa                                                <  1:1139022/93‑1 (MQ=25)
                                        tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGc                                             >  1:5713459/1‑93 (MQ=25)
                                        tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGc                                             >  1:5996184/1‑93 (MQ=25)
                                        tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGc                                             >  2:2740288/1‑93 (MQ=25)
                                        tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGc                                             >  2:4598861/1‑93 (MQ=25)
                                        tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGc                                             >  1:7317746/1‑93 (MQ=25)
                                        tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGc                                             >  1:1877079/1‑93 (MQ=25)
                                         aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCt                                            >  1:2242021/1‑93 (MQ=25)
                                         aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCt                                            <  1:7496355/93‑1 (MQ=25)
                                           tGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTa                                          <  2:2060734/93‑1 (MQ=25)
                                            gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAc                                         >  2:1157098/1‑93 (MQ=25)
                                             aGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAcc                                        <  1:7146641/93‑1 (MQ=25)
                                                gTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCGCTGAAGCCTGATGCGACGCTTACcgcg                                     >  2:2392762/1‑93 (MQ=25)
                                                gTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgcg                                     >  2:2314919/1‑93 (MQ=25)
                                                 ttGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgcgg                                    >  1:5847453/1‑92 (MQ=25)
                                                  tGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTc                                   <  1:3182820/93‑1 (MQ=25)
                                                   gCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCt                                  >  1:86725/1‑93 (MQ=25)
                                                   gCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCt                                  <  1:3429580/93‑1 (MQ=25)
                                                   gCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCt                                  <  1:2932914/93‑1 (MQ=25)
                                                   gCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCt                                  >  1:5375207/1‑93 (MQ=25)
                                                     tATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCTCCCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTa                                >  1:5086182/1‑93 (MQ=21)
                                                       tGCCTGATGCCGATTCGTAGGCCGGATAAGCCGCCCC‑‑CGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATc                              <  1:4059554/93‑1 (MQ=21)
                                                       tGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATc                              >  2:580805/1‑93 (MQ=21)
                                                        gCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCa                             >  2:4892112/1‑93 (MQ=21)
                                                        gCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCa                             >  2:4903233/1‑93 (MQ=21)
                                                         ccTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCAt                            <  2:4214896/93‑1 (MQ=21)
                                                         ccTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCAt                            <  1:6564613/93‑1 (MQ=21)
                                                          cTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATg                           >  1:6443873/1‑93 (MQ=21)
                                                          cTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATg                           >  1:3604922/1‑93 (MQ=21)
                                                           tGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGc                          >  2:6380196/1‑93 (MQ=21)
                                                            gATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGcc                         <  2:5188721/93‑1 (MQ=21)
                                                            gATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGcc                         <  1:3926248/93‑1 (MQ=21)
                                                             aTGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCt                        >  1:2843619/1‑93 (MQ=21)
                                                              tGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTa                       >  2:6825438/1‑93 (MQ=21)
                                                              tGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTa                       >  2:4842893/1‑93 (MQ=21)
                                                              tGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTa                       >  1:1112484/1‑93 (MQ=21)
                                                               gCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTAc                      <  1:1237718/93‑1 (MQ=21)
                                                               gCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTAc                      >  2:1124695/1‑93 (MQ=21)
                                                                 cGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTg                    >  1:5732432/1‑93 (MQ=21)
                                                                  gATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTgg                   <  1:3513714/93‑1 (MQ=21)
                                                                   aTTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTggg                  >  1:2077373/1‑93 (MQ=21)
                                                                    ttCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGa                 >  2:3793551/1‑93 (MQ=18)
                                                                    ttCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGa                 <  1:2736519/93‑1 (MQ=21)
                                                                     tCGTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAg                <  2:7432058/93‑1 (MQ=21)
                                                                       gTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCa              >  1:380861/1‑93 (MQ=21)
                                                                       gTAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCa              <  2:7149799/93‑1 (MQ=21)
                                                                        tAGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCAc             >  2:1285063/1‑93 (MQ=21)
                                                                         aGGCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACg            <  1:7230795/93‑1 (MQ=21)
                                                                          ggCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGc           <  1:1489832/93‑1 (MQ=21)
                                                                          ggCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGc           <  1:5748285/93‑1 (MQ=21)
                                                                           gCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCt          >  1:4579200/1‑93 (MQ=21)
                                                                           gCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCt          <  2:6991389/93‑1 (MQ=21)
                                                                           gCCGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCt          >  2:2640578/1‑93 (MQ=21)
                                                                             cGGATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCttt        >  2:6339686/1‑93 (MQ=21)
                                                                              ggATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTa       <  1:707388/93‑1 (MQ=21)
                                                                               gATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTac      >  2:3768332/1‑93 (MQ=21)
                                                                               gATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTac      <  2:5131447/93‑1 (MQ=25)
                                                                               gATAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTac      >  1:4698327/1‑93 (MQ=21)
                                                                                aTAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTaca     >  2:5612093/1‑93 (MQ=25)
                                                                                 tAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTacac    <  2:1715684/93‑1 (MQ=25)
                                                                                 tAAG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTacac    >  2:5318358/1‑93 (MQ=25)
                                                                                  aaG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACAcc   <  1:2675542/93‑1 (MQ=25)
                                                                                  aaG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACAcc   >  2:3664378/1‑93 (MQ=25)
                                                                                   aG‑‑GCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCg  <  2:6786168/93‑1 (MQ=25)
                                                                                                |                                                                                   
TACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAG‑‑GCGCT‑‑CGGCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCG  >  NC_000913/4448450‑4448623

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