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.7% 1190.6 / NA 375intergenic (+46/‑92)ytfQ/ytfRputative LACI‑type transcriptional regulator/YtfR
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base C (237/137);  minor base A (1/0);  total (238/137)
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 436 total aligned reads displayed.

CTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCG‑‑GCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGG  >  NC_000913/4448446‑4448626
                                                                                                |                                                                                        
ctgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGAGAAGGCGCTCg                                                                                           >  2:1361313/1‑96 (MQ=255)
 tgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCgc                                                                                          <  2:10017675/96‑2 (MQ=255)
 tgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCgc                                                                                          <  2:12855253/96‑2 (MQ=255)
  gATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGTTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCgcg                                                                                         <  2:10948056/96‑3 (MQ=255)
  gATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCgcg                                                                                         >  1:7325707/1‑94 (MQ=255)
  gATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCgcg                                                                                         <  1:1200855/96‑3 (MQ=255)
  gATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCg‑‑gg                                                                                       >  1:14306435/1‑95 (MQ=255)
   aTACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCgcgc                                                                                        <  1:2381396/96‑4 (MQ=255)
    tACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCgcgcc                                                                                       <  1:67850/96‑5 (MQ=255)
    tACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCgcgcc                                                                                       <  2:6259647/96‑5 (MQ=255)
    tACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTCCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCgcgcc                                                                                       <  1:748820/96‑5 (MQ=255)
             agaagaATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGc                                                                            <  1:10373893/96‑1 (MQ=37)
             agaagaATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGc                                                                            <  1:2125342/96‑1 (MQ=37)
             agaagaATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGc                                                                            <  2:8838878/96‑1 (MQ=37)
               aagaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGa                                                                          >  1:1795223/1‑96 (MQ=37)
               aagaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGa                                                                          >  2:5469066/1‑96 (MQ=37)
                 gaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATg                                                                        >  1:4457167/1‑96 (MQ=37)
                  aaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATgg                                                                       >  1:1073264/1‑96 (MQ=37)
                   agtAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGt                                                                      >  1:10839540/3‑96 (MQ=37)
                   aTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGt                                                                      >  2:6086855/1‑96 (MQ=37)
                     tAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGc                                                                    <  2:6371265/96‑1 (MQ=37)
                      agagAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCa                                                                   >  2:2421225/1‑96 (MQ=37)
                           agaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGaa                                                              <  1:13698316/96‑1 (MQ=37)
                           agaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGGGCACTGaa                                                              >  1:11848505/1‑96 (MQ=37)
                              agaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGcc                                                           >  2:13240061/1‑96 (MQ=37)
                              agaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGcc                                                           >  1:5657322/1‑96 (MQ=37)
                              agaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGAGCCGCATCCGGCGATGGTGCACTGAAGcc                                                           >  1:5299550/1‑96 (MQ=37)
                               gaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCt                                                          <  2:4441811/96‑1 (MQ=37)
                                aaaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTg                                                         <  2:102023/96‑1 (MQ=37)
                                aaaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGCCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTg                                                         >  1:8652563/1‑96 (MQ=37)
                                 aaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGa                                                        >  2:5133264/1‑96 (MQ=37)
                                 aaaaTATGGGGTATGGAGGGTTGCTATGCCTGATGCCGATTCGTCGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGa                                                        <  2:12693323/96‑1 (MQ=25)
                                  aaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGAt                                                       >  2:8177090/1‑96 (MQ=37)
                                   aaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATg                                                      <  2:5334097/96‑1 (MQ=37)
                                   aaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATg                                                      >  2:9606454/1‑96 (MQ=37)
                                    atatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGc                                                     >  1:7681146/1‑96 (MQ=37)
                                    atatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGc                                                     >  2:2751125/1‑96 (MQ=25)
                                    atatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGc                                                     >  2:8666148/1‑96 (MQ=37)
                                     tatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCg                                                    <  2:8378489/96‑1 (MQ=37)
                                        ggggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGaga                                                 <  1:11310939/96‑3 (MQ=25)
                                          ggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCt                                               <  2:12142551/96‑1 (MQ=37)
                                           gTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCtt                                              >  1:11544585/1‑96 (MQ=25)
                                           gTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCtt                                              >  2:2889333/1‑96 (MQ=25)
                                           gTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCtt                                              >  1:10559070/1‑96 (MQ=25)
                                           gTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCtt                                              <  2:13080045/96‑1 (MQ=37)
                                           gTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCtt                                              >  2:11858917/1‑96 (MQ=25)
                                           gTATTGAGGGGTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCGTGATGCGACGCtt                                              >  1:6397319/1‑96 (MQ=25)
                                            tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTa                                             >  2:12327824/1‑96 (MQ=25)
                                            tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTa                                             <  2:9200239/96‑1 (MQ=25)
                                             aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAc                                            >  1:8437623/1‑96 (MQ=25)
                                              ttGAGGGTTGCTATCCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCTGCATCCGGCGATGGTGCACTGAAGCCTGCTGCGACGCTTAcc                                           <  1:10977503/96‑1 (MQ=25)
                                               tGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcg                                          >  1:14085747/1‑96 (MQ=25)
                                                taGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgc                                         <  1:3541009/95‑1 (MQ=25)
                                                gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgc                                         <  2:7485202/96‑1 (MQ=25)
                                                 aGGGTTCCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgcg                                        <  2:3919526/96‑1 (MQ=25)
                                                    gTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCt                                     >  2:4862591/1‑96 (MQ=25)
                                                    gTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCt                                     >  2:6751735/1‑96 (MQ=25)
                                                     ttGCTCTGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGCAGCCTGCTGCGACGCTTACCGCGTCtt                                    <  2:1479138/96‑1 (MQ=21)
                                                     ttGCTATCCCTGACGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCtt                                    <  1:3042354/96‑1 (MQ=21)
                                                      tGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTa                                   <  2:11489060/96‑1 (MQ=25)
                                                      tGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGGGCACTGAAGCCTGATGGGACGCTTACCGCGCCTTa                                   >  1:1573467/1‑96 (MQ=25)
                                                       gCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGTTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTAt                                  >  2:10525492/1‑96 (MQ=25)
                                                       gCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTAt                                  >  1:9198704/1‑96 (MQ=25)
                                                         tATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCa                                >  1:7881192/1‑96 (MQ=21)
                                                          aTGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCAt                               <  1:5170831/96‑1 (MQ=25)
                                                          aTGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCAt                               <  1:3195856/96‑1 (MQ=25)
                                                            gCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGc                             <  1:708037/96‑1 (MQ=21)
                                                            gCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGACTGATGCGACGCTTACCGCGTCTTATCATGc                             >  2:7148277/1‑96 (MQ=21)
                                                             ccTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGcc                            <  2:4244071/96‑1 (MQ=21)
                                                              cTGGTGCCGATCCGGAGGCCGAACAAGGCGCCCGCGCCGCATCTGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCt                           <  2:2582573/96‑1 (MQ=18)
                                                              cTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCt                           >  2:11323206/1‑96 (MQ=21)
                                                              cTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCt                           >  1:12484583/1‑96 (MQ=21)
                                                               tGATGCCGATTCGTAGGCCGGATAAGGCACTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTa                          >  2:11428579/1‑96 (MQ=21)
                                                                    ccGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTggg                     >  1:5563509/1‑96 (MQ=21)
                                                                     cGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGa                    >  2:3454781/1‑96 (MQ=21)
                                                                     cGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGa                    >  1:381911/1‑96 (MQ=21)
                                                                     cGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGa                    >  2:1656424/1‑96 (MQ=21)
                                                                     agATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGa                    >  1:8842766/2‑96 (MQ=21)
                                                                      gATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAg                   >  2:3630563/1‑96 (MQ=21)
                                                                      gATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAg                   >  1:6546185/1‑96 (MQ=21)
                                                                      gATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAg                   <  2:7640679/96‑1 (MQ=21)
                                                                       aTTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGc                  <  2:1982017/96‑1 (MQ=21)
                                                                       aTTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGc                  <  1:4122337/96‑1 (MQ=21)
                                                                        ttCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCa                 >  1:10086971/1‑96 (MQ=21)
                                                                          cGTAGGCCTGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCatg               >  1:12960715/1‑94 (MQ=18)
                                                                            tAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCt             >  1:9453601/1‑96 (MQ=21)
                                                                              ggCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCttt           <  2:8007733/96‑1 (MQ=21)
                                                                               gCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTa          <  2:5930737/96‑1 (MQ=21)
                                                                                ccGGATAAGGCGCTCGCCCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTac         <  1:3731504/96‑1 (MQ=25)
                                                                                ccGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTac         >  2:2132711/1‑96 (MQ=21)
                                                                                 cGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTacc        >  1:13217566/1‑95 (MQ=21)
                                                                                  ggATAAGGCGCTCGCGCCGCATCCGTCGATGGTGCACTGAAGCCTGATGCGAGGCTTACCGCGTCTTATCATGCCGACTGGGAGCTCTCTTTacag       >  1:6118565/1‑95 (MQ=21)
                                                                                  ggATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCCCGCTTTacac       <  1:4198809/96‑1 (MQ=21)
                                                                                   gATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACAcc      <  2:10797089/96‑1 (MQ=25)
                                                                                   gATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACAcc      <  2:4706347/96‑1 (MQ=25)
                                                                                      aaGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCCCGCTTTACACCggg   >  2:7974258/1‑96 (MQ=25)
                                                                                      aaGGCGCTCGCGCCGCATCCGGCGATGGCGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCggg   >  2:13916166/1‑96 (MQ=25)
                                                                                       aGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCgggg  <  2:7053219/96‑1 (MQ=25)
                                                                                       aGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCgggg  >  1:8176108/1‑96 (MQ=25)
                                                                                       aGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCgggg  >  1:11751248/1‑96 (MQ=25)
                                                                                                |                                                                                        
CTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCG‑‑GCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGG  >  NC_000913/4448446‑4448626

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: