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.C98.0% 237.0 / ‑5.1 101intergenic (+47/‑91)ytfQ/ytfRputative LACI‑type transcriptional regulator/YtfR
Reads supporting (aligned to +/- strand):  ref base . (0/2);  major base C (62/35);  minor base A (1/1);  total (63/38)
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 101 total aligned reads displayed.

AAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGGCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGG  >  NC_000913/4448478‑4448606
                                                                      |                                                                
aaaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCgg                                                         >  2:4011943/1‑76 (MQ=25)
aaaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCgg                                                         >  1:2513945/1‑76 (MQ=25)
aaaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCgg                                                         <  2:1500932/76‑1 (MQ=25)
aaaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCgg                                                         <  1:5039901/76‑1 (MQ=25)
 aaaaTATGGGGTCTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGc                                                        <  1:2865749/76‑1 (MQ=21)
    atatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGAt                                                     <  1:5515070/76‑1 (MQ=25)
     tatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATg                                                    >  1:3628557/1‑76 (MQ=25)
     tatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATg                                                    >  1:5428812/1‑76 (MQ=25)
      atGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATgg                                                   >  1:3801287/1‑76 (MQ=25)
       tGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGt                                                  <  1:780822/76‑1 (MQ=25)
        ggggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTg                                                 <  1:4116761/76‑1 (MQ=25)
        ggggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTg                                                 <  1:3937658/76‑1 (MQ=25)
         gggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGc                                                >  2:3894258/1‑76 (MQ=21)
         gggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCAGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTgg                                                >  1:4009636/1‑75 (MQ=21)
          ggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCa                                               <  2:3628557/76‑1 (MQ=25)
           gTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCAc                                              <  2:4833038/76‑1 (MQ=25)
            tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCAcc                                             >  2:2444402/1‑75 (MQ=21)
            tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACt                                             >  1:4791131/1‑76 (MQ=25)
            tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACt                                             >  1:1219108/1‑76 (MQ=25)
            tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACt                                             >  1:1933878/1‑76 (MQ=25)
            tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACt                                             <  2:2109806/76‑1 (MQ=21)
             tttGCGGGTTGCTCGGCCTGCTGCCGATTCGTAGTCCGGATAG‑‑GGCGCTGGCGCCGCATCCGTCGATGGTGCACTg                                            <  2:3067682/75‑1 (MQ=255)
             aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCCTCCGGCGATGGTGCACTg                                            >  2:1534859/1‑76 (MQ=21)
             aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTg                                            >  1:615344/1‑76 (MQ=25)
             aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTg                                            <  2:1149141/76‑1 (MQ=21)
             aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTg                                            >  1:4555458/1‑76 (MQ=21)
             aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTg                                            >  2:2581215/1‑76 (MQ=21)
              ttGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGGGCACTGa                                           >  2:3293210/1‑76 (MQ=21)
               tGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGaa                                          <  1:277222/76‑1 (MQ=25)
               tGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGaa                                          >  2:1921352/1‑76 (MQ=25)
               tGAGGGTGGCTATGCCTGATGCCGATCCGTAGGCCGGATAC‑‑GTCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGaa                                          <  2:4745543/76‑1 (MQ=18)
                gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAg                                         >  2:3782567/1‑76 (MQ=25)
                gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAg                                         >  1:2906156/1‑76 (MQ=25)
                gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAg                                         >  1:4816541/1‑76 (MQ=25)
                 aGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGc                                        <  1:3145587/76‑1 (MQ=25)
                 aGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACGGAAGc                                        >  1:91321/1‑76 (MQ=21)
                  gggTTGCTATGCCTGATGCCGATTCGTAGGCCTGATAA‑‑GGCGCT‑‑CGCGCCGCCTCCGGCGATGGTGCACTGAAGcc                                       <  2:5091588/76‑1 (MQ=21)
                  gggTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGcc                                       >  1:3953244/1‑76 (MQ=25)
                  gggTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGCGCACTGCAGcc                                       >  2:5171708/1‑76 (MQ=18)
                   ggTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCt                                      >  1:5051212/1‑76 (MQ=25)
                   ggTTGCCATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCt                                      >  2:5568063/1‑76 (MQ=21)
                    gTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTg                                     >  1:1846946/1‑76 (MQ=25)
                    gTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTg                                     <  1:1827895/76‑1 (MQ=25)
                    gTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTg                                     <  1:1478556/76‑1 (MQ=21)
                     ttGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGa                                    >  1:4439864/1‑76 (MQ=25)
                     ttGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAACCCTGa                                    >  2:1821866/1‑76 (MQ=21)
                      tGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGAt                                   <  1:5348817/76‑1 (MQ=25)
                      tGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGAt                                   >  1:567149/1‑76 (MQ=25)
                      tGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGAt                                   >  1:5433784/1‑76 (MQ=21)
                       gCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATg                                  >  1:5007556/1‑76 (MQ=21)
                       gCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATg                                  >  2:2421207/1‑76 (MQ=21)
                        cTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGc                                 <  1:719982/76‑1 (MQ=21)
                        cTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGc                                 >  1:943166/1‑76 (MQ=21)
                         tATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCg                                <  1:1828541/76‑1 (MQ=21)
                         tATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCg                                >  2:5037567/1‑76 (MQ=21)
                         tATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCg                                >  2:2236077/1‑76 (MQ=21)
                         tATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCg                                >  1:5375621/1‑76 (MQ=21)
                           tGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGAc                              >  2:3377452/1‑76 (MQ=21)
                           tGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGAc                              >  1:989053/1‑76 (MQ=21)
                           tGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACCGAAGCCTGATGCGAc                              <  1:2444402/76‑1 (MQ=18)
                            gCCTGATGCCGATTCGTAGGCCGGATAA‑‑TGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACg                             <  1:3644118/76‑1 (MQ=21)
                             ccTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCCTCCGGCGATGGTGCACTGAAGACTGATGCGACGc                            >  2:3219491/1‑76 (MQ=21)
                              cTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCt                           >  1:763957/1‑76 (MQ=21)
                               tGATGCCGATTGGTAGGCCGGATAG‑‑GTCGCT‑‑CGTGACGCATCCGGCGATGGTGCTCTGAAGCGTGATGCGCCGGtt                          <  2:2969181/76‑1 (MQ=16)
                               tGATGCCGATTCGTAGGCCGGATAA‑‑GTCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCtc                          >  2:2964933/1‑75 (MQ=21)
                               tGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGGGCACTGAAGCCTGATGCGACGCtt                          >  2:2924699/1‑76 (MQ=21)
                                  tGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAcc                       <  2:1475824/76‑1 (MQ=21)
                                    ccGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgc                     <  2:1376116/76‑1 (MQ=21)
                                    ccGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgc                     <  1:2380510/76‑1 (MQ=21)
                                    ccGAGTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgc                     <  1:197966/76‑1 (MQ=18)
                                     cGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGGCGCTTACcgcg                    <  1:4689754/76‑1 (MQ=21)
                                     cGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgcg                    >  1:1781314/1‑76 (MQ=21)
                                       aTTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTc                  >  1:3967232/1‑76 (MQ=21)
                                       aTTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTc                  <  1:5410439/76‑1 (MQ=21)
                                        ttagTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTCAAGCCTGATGCTACGCTTACCGCGTCt                 >  1:3063118/4‑76 (MQ=18)
                                        ttCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGCCGCATACCGCGGCt                 >  2:3014735/1‑76 (MQ=17)
                                        ttCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCt                 >  2:3703990/1‑76 (MQ=21)
                                         tCGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCtt                >  1:3697576/1‑76 (MQ=21)
                                          cGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTa               <  1:4621401/76‑1 (MQ=21)
                                          cGTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTa               >  1:4341491/1‑76 (MQ=21)
                                           gTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTAt              >  1:192332/1‑76 (MQ=21)
                                           gTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTAt              <  1:1146700/76‑1 (MQ=21)
                                           gTAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTAt              >  1:1795796/1‑76 (MQ=18)
                                            tAGGCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATc             >  1:1760812/1‑76 (MQ=18)
                                              ggCCGGATAAGGGGCGCG‑‑CGCCGCGTCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCAt           <  2:405406/76‑1 (MQ=18)
                                              ggCCGGATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCAt           >  2:3363195/1‑76 (MQ=18)
                                                ccGGATGA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGc         <  1:128038/76‑1 (MQ=18)
                                                ccGGATAA‑‑GGCGCT‑‑CGCGCCGCACCCGGCGATGGTGCACTGAAGCCTGATGCGGCGCTTACCGCGTCTTATCATGc         <  1:186818/76‑1 (MQ=18)
                                                  ggATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCt       >  1:4020962/1‑76 (MQ=21)
                                                  ggATAA‑‑GGAGCT‑‑CGCGACGCATCCGGCGATGGTGCACGGAAGCATGATGCGTCGCGTACCGCGTCGTATCATGCCt       >  2:3437730/1‑76 (MQ=19)
                                                   gATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTa      >  1:3197761/1‑76 (MQ=21)
                                                   gATAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTa      <  1:4094571/76‑1 (MQ=21)
                                                    gtAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTAc     <  2:616330/75‑1 (MQ=18)
                                                    aTAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTAc     >  1:83871/1‑76 (MQ=21)
                                                    aTAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGGGCACTGAAGCCTGATGCTACGCTTACCGCGTCCTATCATGCCTAc     >  1:859859/1‑76 (MQ=18)
                                                    aTAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGAGGGTGCACTGAAGACTGATGCGACGCTTACCGCGTCTTATCATGCCTAc     >  1:5478865/1‑76 (MQ=18)
                                                     tAA‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGGCGCTTACCGCGTCTTATCATGCCCACt    >  2:4097087/1‑76 (MQ=18)
                                                      aa‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTg   >  1:3044619/1‑76 (MQ=21)
                                                       a‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCCCTCAAGCCTGCTGCGACGCTTACCGCGACTTATCATGCCTACTgg  >  2:3819875/1‑76 (MQ=18)
                                                       a‑‑GGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTgg  <  2:2906156/76‑1 (MQ=21)
                                                                      |                                                                
AAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAA‑‑GGCGCT‑‑CGGCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGG  >  NC_000913/4448478‑4448606

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