Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,080,532 A→G R957R (CGT→CGC)  fdoG ← formate dehydrogenase‑O, major subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,080,5320AG99.6% 817.9 / ‑4.9 257R957R (CGT→CGC) fdoGformate dehydrogenase‑O, major subunit
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (145/111);  total (145/112)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

Only 100 of 266 total aligned reads displayed.

AAAGCCTTTTTTAGCTACGCCTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATACGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTATCGCCCTGGGCAATGCCAAGTTTATTCGC  >  NC_000913/4080442‑4080621
                                                                                          |                                                                                         
 aaGCCTTTTTTAGCTACGCCTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGt                                                                                        <  1:516066/93‑1 (MQ=255)
  aGCCTTTTTTAGCTACGCCTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGtt                                                                                       <  1:2196388/93‑1 (MQ=255)
   gCCTTTTTTAGCTACGCCTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGttt                                                                                      >  1:1359824/1‑93 (MQ=255)
   gCCTTTTTTAGCTACGCCTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGttt                                                                                      >  2:1468358/1‑93 (MQ=255)
   gCCTTTTTTAGCTACGCCTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGttt                                                                                      >  2:1311917/1‑93 (MQ=255)
    ccTTTTTTAGCTACGCCTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTg                                                                                     >  1:3840170/1‑93 (MQ=255)
    ccTTTTTTAGCTACGCCTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTg                                                                                     >  2:525542/1‑93 (MQ=255)
     cTTTTTTAGCTACGCCTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTgg                                                                                    >  2:1076920/1‑93 (MQ=255)
     cTTTTTTAGCTACGCCTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTgg                                                                                    >  1:1887380/1‑93 (MQ=255)
       tttttAGCTACGCCTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTc                                                                                  >  1:2609894/1‑93 (MQ=255)
       tttttAGCTACGCCTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTc                                                                                  <  1:3864072/93‑1 (MQ=255)
        ttttAGCTACGCCTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTcc                                                                                 >  1:478030/1‑92 (MQ=255)
        ttttAGCTACGCCTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTca                                                                                 >  1:14248/1‑93 (MQ=255)
        ttttAGCTACGCCTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGGTTGGTca                                                                                 >  2:4270160/1‑93 (MQ=255)
         tttAGCTACGCCTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTcac                                                                                >  1:4125453/1‑93 (MQ=255)
         tttAGCTACGCCTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTcac                                                                                >  1:2978768/1‑93 (MQ=255)
            aGCTACGCCTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTcaccac                                                                             >  1:3236807/1‑93 (MQ=255)
             gCTACGCCTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCaccacc                                                                            >  2:3489482/1‑93 (MQ=255)
              cTACGCCTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCacaaccc                                                                           >  1:226741/1‑92 (MQ=255)
              cTACGCCTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCg                                                                           >  1:4325751/1‑93 (MQ=255)
               tACGCCTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCgc                                                                          >  1:4301170/1‑93 (MQ=255)
               tACGCCTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCgc                                                                          >  1:2846110/1‑93 (MQ=255)
                 cGCCTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCgctt                                                                        <  2:1102093/93‑1 (MQ=255)
                   ccTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCgctttg                                                                      >  1:309702/1‑93 (MQ=255)
                   ccTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCgctttg                                                                      <  1:3476004/93‑1 (MQ=255)
                    cTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCgctttgg                                                                     <  1:2456833/93‑1 (MQ=255)
                    cTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCgctttgg                                                                     >  2:323983/1‑93 (MQ=255)
                    cTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCgctttgg                                                                     >  2:2280636/1‑93 (MQ=255)
                     ttCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCgctttggc                                                                    <  1:4579794/93‑1 (MQ=255)
                     ttCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCgctttggc                                                                    >  2:2636722/1‑93 (MQ=255)
                       cATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCgctttggctt                                                                  <  1:861646/93‑1 (MQ=255)
                        aTAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCgctttggcttt                                                                 >  1:1659345/1‑93 (MQ=255)
                        aTAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCgctttggcttt                                                                 <  1:1507772/93‑1 (MQ=255)
                        aTAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCgctttggcttt                                                                 <  1:941179/93‑1 (MQ=255)
                        aTAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCCCTTTCAGCGTGCGAATGCGTTTGGTCACCACCgctttggcttt                                                                 <  2:747465/93‑1 (MQ=255)
                            ccccAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGata                                                             <  1:744515/93‑1 (MQ=255)
                             cccAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGatat                                                            >  1:1754894/1‑93 (MQ=255)
                             cccAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGatat                                                            <  2:4487131/93‑1 (MQ=255)
                             cccAGTGAATCGGAATACCGATGGGATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGatat                                                            >  2:1649084/1‑93 (MQ=255)
                              ccAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGAtata                                                           <  2:1380111/93‑1 (MQ=255)
                              ccAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGAtata                                                           >  2:3186239/1‑93 (MQ=255)
                               cAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAg                                                          >  1:599819/1‑93 (MQ=255)
                                aGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGc                                                         >  2:562501/1‑93 (MQ=255)
                                  tGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCa                                                       <  1:3696431/93‑1 (MQ=255)
                                  tGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCa                                                       <  2:2853248/93‑1 (MQ=255)
                                  tGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCa                                                       >  2:701490/1‑93 (MQ=255)
                                   gAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCAc                                                      >  2:2535612/1‑93 (MQ=255)
                                   gAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCAc                                                      <  1:1328598/93‑1 (MQ=255)
                                      tCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGt                                                   <  2:2776457/93‑1 (MQ=255)
                                      tCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGt                                                   <  2:4321937/93‑1 (MQ=255)
                                       cGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGtt                                                  <  2:1851561/93‑1 (MQ=255)
                                         gAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTgg                                                <  2:3221019/93‑1 (MQ=255)
                                         gAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTgg                                                <  2:80091/93‑1 (MQ=255)
                                          aaTACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTgga                                               >  2:890663/1‑93 (MQ=255)
                                           aTACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTggag                                              <  1:107554/93‑1 (MQ=255)
                                           aTACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTggag                                              >  2:211100/1‑93 (MQ=255)
                                            tACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTggagg                                             >  1:3445612/1‑93 (MQ=255)
                                             aCCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTggagga                                            <  2:2197051/93‑1 (MQ=255)
                                              ccGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGgaggag                                           <  2:301386/93‑1 (MQ=255)
                                                 aTGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACt                                        <  2:4686155/93‑1 (MQ=255)
                                                 aTGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACt                                        <  2:3653078/93‑1 (MQ=255)
                                                  tGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACtt                                       <  1:4238000/93‑1 (MQ=255)
                                                  tGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACtt                                       >  2:4184982/1‑93 (MQ=255)
                                                   ggTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACttt                                      >  2:1825685/1‑93 (MQ=255)
                                                    gTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGAGTTTc                                     >  2:2226682/1‑93 (MQ=255)
                                                       tCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGGGTTTGGTCACGACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACg                                  >  2:385834/1‑93 (MQ=255)
                                                         gATATCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGt                                >  1:2933996/1‑93 (MQ=255)
                                                          atatCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTa                               >  1:4009298/1‑93 (MQ=255)
                                                          atatCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTa                               >  1:4384591/1‑93 (MQ=255)
                                                           tatCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTAt                              >  2:4633251/1‑93 (MQ=255)
                                                           tatCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTAt                              >  1:1013050/1‑93 (MQ=255)
                                                            atCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTATc                             <  2:3367781/93‑1 (MQ=255)
                                                            atCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTATc                             >  1:632577/1‑93 (MQ=255)
                                                             tCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTATCg                            <  2:50816/93‑1 (MQ=255)
                                                             tCTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCAAGGTATCg                            >  2:739639/1‑93 (MQ=255)
                                                              cTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTATCGc                           <  2:3122747/93‑1 (MQ=255)
                                                              cTTTGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACAACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTATCGc                           >  1:1616314/1‑93 (MQ=255)
                                                               tttGCCGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTATCGcc                          >  2:3674449/1‑93 (MQ=255)
                                                                   ccGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTATCGCCCTgg                      <  2:3571091/93‑1 (MQ=255)
                                                                    cGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTATCGCCCTggg                     <  1:4650510/93‑1 (MQ=255)
                                                                    cGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTATCGCCCTggg                     >  1:1985786/1‑93 (MQ=255)
                                                                    cGTTCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTATCGCCCTggg                     <  1:188416/93‑1 (MQ=255)
                                                                      ttCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTATCGCCCTGGGCa                   >  2:191104/1‑93 (MQ=255)
                                                                      ttCGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTATCGCCCCGGGCa                   >  1:1251933/1‑93 (MQ=255)
                                                                        cGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTATCGCCCTGGGCAAt                 <  1:3557102/93‑1 (MQ=255)
                                                                        cGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTATCGCCCTGGGCAAt                 <  2:200839/93‑1 (MQ=255)
                                                                        cGCTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTATCGCCCTGGGCAAt                 <  1:319133/93‑1 (MQ=255)
                                                                          cTTTCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTATCGCCCTGGGCAATGc               <  1:3306395/93‑1 (MQ=255)
                                                                           tttCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTATCGCCCTGGGCAATGcc              <  2:2734126/93‑1 (MQ=255)
                                                                           tttCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTATCGCCCTGGGCAATGcc              >  1:3948761/1‑93 (MQ=255)
                                                                            ttCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTATCGCCCTGGGCAATGCCa             <  2:1287806/93‑1 (MQ=255)
                                                                             tCAGCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTATCGCCCTGGGCAATGCCaa            <  2:2822590/93‑1 (MQ=255)
                                                                                gCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTATCGCCCTGGGCAATGCCAAGtt         >  1:1431744/1‑93 (MQ=255)
                                                                                gCGTGCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTATCGCCCTGGGCAATGCCAAGtt         >  2:2970841/1‑93 (MQ=255)
                                                                                    gCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTATCGCCCTGGGCAATGCCAAGTTTAtt     <  2:4097819/93‑1 (MQ=255)
                                                                                    gCGAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTATCGCCCTGGGCAATGCCAAGTTTAtt     <  1:3813659/93‑1 (MQ=255)
                                                                                      gAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTATCGCCCTGGGCAATGCCAAGTTTATTCg   >  1:3886400/1‑93 (MQ=255)
                                                                                      gAATGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTATCGCCCTGGGCAATGCCAAGTTTATTCg   >  1:3144587/1‑93 (MQ=255)
                                                                                       aaTGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGTAGGAGACTTTCACGGTATCGCCCTGGGCAATGCCAAGTTTATTCGc  <  1:2164247/93‑1 (MQ=255)
                                                                                       aaTGCGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTATCGCCCTGGGCAATGCCAAGTTTATTCGc  >  1:3584280/1‑93 (MQ=255)
                                                                                          |                                                                                         
AAAGCCTTTTTTAGCTACGCCTTCATAGCCCCAGTGAATCGGAATACCGATGGTATCGATATCTTTGCCGTTCGCTTTCAGCGTGCGAATACGTTTGGTCACCACCGCTTTGGCTTTGATATAGCCACGGTTGGAGGAGACTTTCACGGTATCGCCCTGGGCAATGCCAAGTTTATTCGC  >  NC_000913/4080442‑4080621

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