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% 343.3 / ‑3.8 120intergenic (+47/‑91)ytfQ/ytfRputative LACI‑type transcriptional regulator/YtfR
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (64/55);  total (64/56)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 127 total aligned reads displayed.

AAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGGCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCG  >  NC_000913/4448461‑4448623
                                                                                   |                                                                                  
aagaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCgg                                                                          >  2:1356555/1‑93 (MQ=37)
  gaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCg                                                                        >  2:2989429/1‑93 (MQ=37)
  gaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCg                                                                        >  1:4074328/1‑93 (MQ=37)
   aaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGa                                                                       >  1:2896154/1‑93 (MQ=37)
   aaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGa                                                                       <  2:1856431/93‑1 (MQ=37)
    aTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGAt                                                                      >  2:3571215/1‑93 (MQ=37)
     ttAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATg                                                                     >  1:528330/1‑93 (MQ=37)
     ttAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATg                                                                     <  2:2180137/93‑1 (MQ=37)
     ttAGAGAAGAAGAAAAATATGGGGAATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATg                                                                     >  2:3267301/1‑93 (MQ=25)
      tAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATgg                                                                    <  1:2705302/93‑1 (MQ=37)
      tAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATgg                                                                    <  1:2453252/93‑1 (MQ=37)
       agcgaagtagaaAACTATGGGGAATTGAAGGGTGCTATGCCTGATGCCGATTAGTAGGCCGGATAAGGCGCTCGCGCCCAGTCCGGCGATggg                                                                   >  2:3758891/4‑92 (MQ=19)
       agagAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGt                                                                   >  2:2838956/1‑93 (MQ=37)
       agagAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGt                                                                   >  2:3758892/1‑93 (MQ=37)
       agagAAGAAGAAAAATATGGGGCATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGt                                                                   >  1:2025859/1‑93 (MQ=25)
        gagaAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTg                                                                  >  1:4369816/1‑93 (MQ=37)
        gagaAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTg                                                                  >  1:3037463/1‑93 (MQ=37)
           aagaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCAc                                                               >  1:2555288/1‑93 (MQ=37)
            agaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACt                                                              <  1:3966012/93‑1 (MQ=37)
            agaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACt                                                              >  2:2132841/1‑93 (MQ=37)
            agaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACt                                                              <  1:277201/93‑1 (MQ=37)
             gaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTg                                                             >  2:741458/1‑93 (MQ=37)
              aagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGa                                                            <  1:1438459/93‑1 (MQ=37)
               agaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGaa                                                           >  1:3242038/1‑93 (MQ=37)
               agaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGaa                                                           >  1:3242036/1‑93 (MQ=37)
                gaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAg                                                          <  2:2472126/93‑1 (MQ=37)
                 aaaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGc                                                         >  1:260879/1‑93 (MQ=37)
                 aaaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGc                                                         <  1:547422/93‑1 (MQ=37)
                 aaaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGc                                                         >  1:508425/1‑93 (MQ=37)
                  aaaaTATGGGGTATGGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGcc                                                        <  2:2292361/93‑1 (MQ=25)
                   aaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCt                                                       <  2:425608/93‑1 (MQ=37)
                    aaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTg                                                      >  1:3833138/1‑93 (MQ=37)
                     atatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGa                                                     >  1:3772709/1‑93 (MQ=37)
                      tatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGAt                                                    >  1:4170326/1‑93 (MQ=37)
                      tatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGAt                                                    <  2:786810/93‑1 (MQ=37)
                      tatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGAt                                                    >  1:215105/1‑93 (MQ=25)
                       atGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATg                                                   >  1:1105365/1‑93 (MQ=25)
                       atGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATg                                                   <  2:1995073/93‑1 (MQ=37)
                          gggTATTGAGTTTTGCTATTCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGa                                                <  1:2257442/93‑1 (MQ=21)
                          gggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGa                                                <  1:2037490/93‑1 (MQ=25)
                             tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGc                                             >  2:4309385/1‑93 (MQ=25)
                              aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCt                                            >  1:3986171/1‑93 (MQ=25)
                              aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCt                                            >  1:939264/1‑93 (MQ=25)
                               ttGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCtt                                           <  1:567366/93‑1 (MQ=25)
                                  aGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAcc                                        >  1:892464/1‑93 (MQ=25)
                                   gggTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcg                                       >  2:2621671/1‑93 (MQ=25)
                                      ttGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGt                                    <  2:4309331/93‑1 (MQ=25)
                                      ttGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGt                                    >  1:3709662/1‑93 (MQ=25)
                                      ttGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGt                                    >  1:3709663/1‑93 (MQ=25)
                                       tGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTc                                   <  1:2450460/93‑1 (MQ=25)
                                       tGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTc                                   >  1:2980682/1‑93 (MQ=25)
                                        tcTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCTCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCt                                  <  1:1906558/92‑1 (MQ=25)
                                        gCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCt                                  <  1:1561812/93‑1 (MQ=25)
                                         cTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCtt                                 <  2:1944178/93‑1 (MQ=25)
                                         cTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCtt                                 <  1:3204192/93‑1 (MQ=25)
                                          tATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTa                                >  1:5490/1‑93 (MQ=21)
                                          tATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTa                                >  1:253266/1‑93 (MQ=25)
                                          tATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTa                                >  1:4649264/1‑93 (MQ=21)
                                            tGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATc                              <  2:240992/93‑1 (MQ=25)
                                            tGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATc                              >  1:172199/1‑93 (MQ=21)
                                             gCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCa                             >  2:1189977/1‑93 (MQ=21)
                                             gCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCa                             >  2:2236947/1‑93 (MQ=21)
                                             gCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCa                             >  2:4581493/1‑93 (MQ=21)
                                              ccTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCAt                            <  1:2083303/93‑1 (MQ=21)
                                              ccTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCAt                            <  2:283195/93‑1 (MQ=21)
                                               cTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATg                           >  1:1585595/1‑93 (MQ=21)
                                               cTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATg                           >  1:2862924/1‑93 (MQ=21)
                                               cTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATg                           >  2:2912015/1‑93 (MQ=21)
                                                tGATGCCGATTCGTAGGCGGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGc                          <  2:1015942/93‑1 (MQ=21)
                                                tGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGc                          <  1:1773001/93‑1 (MQ=21)
                                                 gATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGcc                         <  2:3782634/93‑1 (MQ=21)
                                                  aTGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCt                        <  2:468634/93‑1 (MQ=21)
                                                  aTGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCt                        <  2:1259964/93‑1 (MQ=21)
                                                  aTGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCt                        >  2:3543074/1‑93 (MQ=21)
                                                   tGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTa                       <  1:4431106/93‑1 (MQ=21)
                                                    gCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTAc                      >  1:3403010/1‑93 (MQ=21)
                                                        aTTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTggg                  <  1:130751/93‑1 (MQ=21)
                                                        aTTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTggg                  >  2:1688134/1‑93 (MQ=21)
                                                           cGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGc               <  2:4004477/93‑1 (MQ=21)
                                                            gTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCa              >  1:4335397/1‑93 (MQ=21)
                                                             tAGGCGGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCAc             <  2:3555035/93‑1 (MQ=21)
                                                             tAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCAc             <  1:3996655/93‑1 (MQ=21)
                                                             tAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCAc             <  2:1644583/93‑1 (MQ=21)
                                                             tAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCAc             <  1:4079938/93‑1 (MQ=21)
                                                              aGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACg            >  2:780588/1‑93 (MQ=21)
                                                              aGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGGGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGTGAGCACg            >  1:946165/1‑93 (MQ=18)
                                                               ggCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGc           <  1:4531120/93‑1 (MQ=21)
                                                                gCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCt          <  2:3733495/93‑1 (MQ=21)
                                                                  cGGATAGGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCttt        <  2:1559/93‑1 (MQ=21)
                                                                  cGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCttt        >  1:4247790/1‑93 (MQ=21)
                                                                  cGGATAAGGCCCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCttt        <  2:3615565/93‑1 (MQ=21)
                                                                   ggATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTa       <  2:1172559/93‑1 (MQ=25)
                                                                   ggATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTa       <  1:126720/93‑1 (MQ=21)
                                                                    gATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTac      >  2:1717578/1‑93 (MQ=25)
                                                                    gATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTac      >  2:4042624/1‑93 (MQ=21)
                                                                     aTAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTaca     >  1:3571987/1‑93 (MQ=25)
                                                                      tAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTacac    <  2:1619910/93‑1 (MQ=25)
                                                                      tAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTacac    >  1:1299632/1‑93 (MQ=25)
                                                                       aaGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACAcc   <  2:323595/93‑1 (MQ=25)
                                                                        aGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCg  <  2:3124569/93‑1 (MQ=25)
                                                                                   |                                                                                  
AAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGGCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCG  >  NC_000913/4448461‑4448623

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