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.9% 823.0 / 1.2 283intergenic (+47/‑91)ytfQ/ytfRputative LACI‑type transcriptional regulator/YtfR
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base C (182/98);  total (184/99)
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.61e-01

Only 100 of 304 total aligned reads displayed.

AAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CG‑‑G‑‑CGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGG  >  NC_000913/4448461‑4448624
                                                                                      |                                                                                    
aagaaTTAGAGAAGAAGAAAAACATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCgg                                                                           >  2:5525875/1‑93 (MQ=37)
 agaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGc                                                                          >  2:3569382/1‑93 (MQ=37)
 agaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGc                                                                          >  1:3940437/1‑93 (MQ=37)
  gaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCg                                                                         >  2:318693/1‑93 (MQ=37)
   agttAGAGAAAAAGAAAAATAGGGGGTATTGAGGGTTGCTATGCTTGATGCAGCTTCGTAGGCAGGATAAGGCCCT‑‑CGCGCCGCATCCGGCGa                                                                        <  2:6569639/91‑1 (MQ=19)
   aaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGa                                                                        <  2:8510249/93‑1 (MQ=37)
   aaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGa                                                                        >  2:2781427/1‑93 (MQ=37)
    aTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGAt                                                                       >  2:4895387/1‑93 (MQ=37)
     ttAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATg                                                                      >  1:3486067/1‑93 (MQ=37)
      tAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATgg                                                                     >  2:262036/1‑93 (MQ=37)
       agagAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGt                                                                    >  1:7401540/1‑93 (MQ=37)
       agagAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGt                                                                    >  1:1676856/1‑93 (MQ=37)
        gagaAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTg                                                                   >  2:6755069/1‑93 (MQ=37)
          gaagaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCa                                                                 >  2:4644095/1‑93 (MQ=37)
           aagaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCAc                                                                <  2:5910265/93‑1 (MQ=37)
           aagaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCAc                                                                >  2:6336071/1‑93 (MQ=37)
             gaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTg                                                              >  2:1075839/1‑93 (MQ=37)
             gaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTg                                                              >  1:5255103/1‑93 (MQ=37)
             gaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTg                                                              >  1:2146213/1‑93 (MQ=37)
               agaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGaa                                                            >  1:5721958/1‑93 (MQ=37)
               agaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGaa                                                            <  1:548672/93‑1 (MQ=37)
                gaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAg                                                           >  2:2426577/1‑93 (MQ=37)
                aaaaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGGTAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAg                                                           <  1:6420096/92‑1 (MQ=37)
                 aaaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGc                                                          <  2:4076569/93‑1 (MQ=37)
                 aaaaaTTGGGATATTGGGGTTTGCTTTGCCTGGTGCCATTTCGTAGGCAGGATAGGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGcc                                                         <  2:6670500/93‑1 (MQ=12)
                  aaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGcc                                                         <  1:8165233/93‑1 (MQ=37)
                    aaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTg                                                       <  1:4754095/93‑1 (MQ=37)
                    aaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTg                                                       >  1:7023486/1‑93 (MQ=37)
                     atatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGa                                                      >  1:6215627/1‑93 (MQ=37)
                     atatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGa                                                      <  1:886804/93‑1 (MQ=37)
                      tatGTGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGAt                                                     >  1:6012642/1‑93 (MQ=25)
                      tatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCAC‑‑CGCATCCGGCGATGGTGCACTGAAGCCTGAt                                                     >  2:1488424/1‑93 (MQ=37)
                        tGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGc                                                   >  1:7536041/1‑93 (MQ=25)
                        tGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGc                                                   >  2:1697984/1‑93 (MQ=25)
                        tGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGc                                                   <  2:4133560/93‑1 (MQ=25)
                          gggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGa                                                 >  1:4977551/1‑93 (MQ=25)
                          gggTATGGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGa                                                 <  1:3346380/93‑1 (MQ=25)
                           ggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGag                                                >  2:3034822/1‑92 (MQ=25)
                           ggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGAc                                                >  1:3062627/1‑93 (MQ=25)
                             tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACgt                                              >  1:2504513/1‑92 (MQ=25)
                              aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCt                                             >  2:8710589/1‑93 (MQ=25)
                               ttGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCtt                                            >  2:8604866/1‑93 (MQ=25)
                               ttGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCtt                                            <  1:5585427/93‑1 (MQ=25)
                                tGAGGTTTGCTAGGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCG‑‑CGGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTa                                           <  2:3935510/93‑1 (MQ=25)
                                 gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAc                                          >  2:7737596/1‑93 (MQ=25)
                                  aGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAcc                                         >  1:4161228/1‑93 (MQ=25)
                                  aGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAcc                                         >  2:6939910/1‑93 (MQ=25)
                                  aGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAcc                                         <  1:3975537/93‑1 (MQ=25)
                                   gggTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcg                                        >  2:3647605/1‑93 (MQ=25)
                                    ggTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgc                                       <  1:2707275/93‑1 (MQ=25)
                                      ttGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGt                                     <  1:1449261/93‑1 (MQ=25)
                                         cTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCtt                                  >  1:6544935/1‑93 (MQ=25)
                                          tATGCCTGATGCCGATTCGTAGGCCGTATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTa                                 <  1:4078246/93‑1 (MQ=21)
                                          tATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGGGATGGTGCACTGAAGCCTGAGGGGACGCTTCCCGCGTCttt                                 >  2:4164535/1‑92 (MQ=21)
                                          tATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTa                                 >  2:7297823/1‑93 (MQ=25)
                                           ttGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTAt                                <  1:5091800/92‑1 (MQ=25)
                                            tGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGAGGCTTACCGCGTCTTAta                               >  1:7196502/1‑92 (MQ=21)
                                            tGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATc                               >  2:3838376/1‑93 (MQ=21)
                                             gCCTGATGCCGATTCGTAGGCCTGATAAGGCGCT‑‑CG‑‑GGCCGCATCCGGCGATGGGGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCa                              <  2:2106233/93‑1 (MQ=25)
                                             gCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCa                              >  1:77976/1‑93 (MQ=21)
                                             gCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCa                              >  2:7416376/1‑93 (MQ=21)
                                              ccTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCAt                             <  1:4611200/93‑1 (MQ=21)
                                              ccTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCAt                             >  2:2498373/1‑93 (MQ=21)
                                               cTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATg                            <  2:5264909/93‑1 (MQ=21)
                                               cTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATg                            >  1:1164548/1‑93 (MQ=21)
                                               cTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATg                            >  1:3427103/1‑93 (MQ=21)
                                                tGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCTACGCTTACCGCGTCTTATCATGc                           >  2:4762589/1‑93 (MQ=21)
                                                 gATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGcc                          >  2:6184649/1‑93 (MQ=21)
                                                 gATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGcc                          <  1:835248/93‑1 (MQ=21)
                                                  aTGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCt                         >  1:7945910/1‑93 (MQ=21)
                                                    gCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGATGCACTGAAGCCTGATGAGACGCTTACCGCGTCTTCTCATGCCTAc                       >  2:4313346/1‑93 (MQ=18)
                                                    gCCGATGCGTAGGCCGGATAAGGCGCG‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTAc                       <  1:2046844/93‑1 (MQ=18)
                                                     ccGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACt                      <  1:4238720/93‑1 (MQ=21)
                                                     ccGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACt                      >  2:2127818/1‑93 (MQ=21)
                                                      cGATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTg                     >  1:8688923/1‑93 (MQ=21)
                                                       gATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTgg                    >  1:7718637/1‑93 (MQ=21)
                                                       gATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTgg                    >  2:5196324/1‑93 (MQ=21)
                                                       gATTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTgg                    >  2:6417215/1‑93 (MQ=18)
                                                        aTTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTggg                   <  1:8524689/93‑1 (MQ=21)
                                                        aTTCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTggg                   <  2:2191645/93‑1 (MQ=18)
                                                         ttggTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGa                  <  1:1915117/90‑1 (MQ=18)
                                                         ttCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGa                  >  1:2810437/1‑93 (MQ=21)
                                                         ttCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGa                  <  2:6515198/93‑1 (MQ=21)
                                                          tCGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAg                 <  1:444971/93‑1 (MQ=18)
                                                           cGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGc                >  2:3320827/1‑93 (MQ=21)
                                                           cGTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGc                >  2:7116679/1‑93 (MQ=21)
                                                            gTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCa               >  1:3667557/1‑93 (MQ=21)
                                                            gTAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCa               >  1:1310128/1‑93 (MQ=21)
                                                             tAGTCCGTATAAGGCGCTCGCG‑‑G‑‑CGCATCCGGCGATGGTGCAGTGAACCCTGATGCGACGCTTACCTCGTCTTATCATGCCTACTGGGAGCAc              <  1:4304333/93‑1 (MQ=21)
                                                             tAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCAc              >  1:8564414/1‑93 (MQ=21)
                                                             tAGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCAc              <  1:1033167/93‑1 (MQ=21)
                                                              aGGCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACg             <  1:2925072/93‑1 (MQ=21)
                                                                gacGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCt           <  1:6836651/91‑1 (MQ=21)
                                                                gCCGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCt           >  2:1031794/1‑93 (MQ=21)
                                                                  cGGATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCttt         <  2:8274035/93‑1 (MQ=21)
                                                                    gATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTac       >  2:598337/1‑93 (MQ=21)
                                                                    gATAAGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTac       >  1:7321512/1‑93 (MQ=21)
                                                                      aaaGGCGCT‑‑CGCGCCGCATCGGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTacac     <  2:8002522/92‑1 (MQ=25)
                                                                       ggcgcgcg‑‑cg‑‑c‑‑CGCAGCCGGCGATGGGGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCgg  <  2:2231306/90‑1 (MQ=37)
                                                                        aGGCGCT‑‑CGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCg   >  2:1303527/1‑93 (MQ=25)
                                                                                      |                                                                                    
AAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CG‑‑G‑‑CGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGG  >  NC_000913/4448461‑4448624

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