Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,389,812 Δ1 bp intergenic (‑149/+4151) tldD ← / ← cafA TldD/CafA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,389,8110C.99.3% 495.3 / ‑3.6 146intergenic (‑148/+4152)tldD/cafATldD/CafA
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (71/74);  total (72/74)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 170 total aligned reads displayed.

CGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC  >  NC_000913/3389719‑3389848
                                                                                            |                                     
cggcgTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGc                                       <  1:5814976/93‑1 (MQ=255)
cggcgTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGc                                       <  1:9791399/93‑1 (MQ=255)
 ggcgTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGcg                                      >  1:2165825/1‑92 (MQ=255)
 ggcgTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGcg                                      <  1:8285590/93‑2 (MQ=255)
  gcgTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGcgg                                     >  2:5383765/1‑91 (MQ=255)
  gcgTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGcgg                                     >  1:7984731/1‑91 (MQ=255)
  gcgTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGcgg                                     <  1:2038832/93‑3 (MQ=255)
  gcgTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGcgg                                     >  1:6369511/1‑91 (MQ=255)
   cgTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGcggc                                    <  2:6136460/93‑4 (MQ=255)
   cgTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCggg                                    >  2:2914547/1‑93 (MQ=255)
     tGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGt                                 >  1:6559405/1‑93 (MQ=255)
     tGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGt                                 >  1:8100246/1‑93 (MQ=255)
       aaCGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTaa                               <  2:1866565/93‑1 (MQ=255)
        aCGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTaaa                              >  1:6864292/1‑93 (MQ=255)
        aCGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTaaa                              <  2:354833/93‑1 (MQ=255)
        aCGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTaaa                              <  1:5437071/93‑1 (MQ=255)
        aCGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTaaa                              <  1:5437070/93‑1 (MQ=255)
        aCGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTaaa                              >  1:3083033/1‑93 (MQ=255)
        aCGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTaaa                              <  2:2644296/93‑1 (MQ=255)
         cGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAAc                             <  2:1769314/93‑1 (MQ=255)
          gCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACg                            >  1:2112146/1‑93 (MQ=255)
           cccTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGc                           >  1:558647/1‑93 (MQ=255)
           cccTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGc                           >  1:1905493/1‑93 (MQ=255)
           cccTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGc                           <  2:2260986/93‑1 (MQ=255)
           cacTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGc                           >  2:959362/3‑93 (MQ=255)
            ccTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGcc                          >  1:2422889/1‑93 (MQ=255)
            ccTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGcc                          >  1:4549600/1‑93 (MQ=255)
            ccTATCCAGCCTACGGTTATGTTCCGGTTGGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGcc                          <  2:6764134/93‑1 (MQ=255)
            ccTATCCAGCCTACAGTTATGTACCGGTAAGTAGGCCTGATAAGACGCACGGCGTCACATCAAACAACGGCAGTCGGATG‑CGGCGTAAACGcc                          <  2:9251466/93‑1 (MQ=14)
             cTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATACGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCt                         <  1:8949712/93‑1 (MQ=255)
             cTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCt                         <  2:2457015/93‑1 (MQ=255)
             cTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCt                         >  1:5978620/1‑93 (MQ=255)
             cTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCt                         <  2:6373607/93‑1 (MQ=255)
              tATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGccct                        >  1:3268567/1‑91 (MQ=255)
              tATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCtt                        <  2:2422012/93‑1 (MQ=255)
              tATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCtt                        >  2:4274940/1‑93 (MQ=255)
              tATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCtt                        <  1:6037384/93‑1 (MQ=255)
               aTCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTa                       <  2:6095240/93‑1 (MQ=255)
               aTCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTa                       <  2:52885/93‑1 (MQ=255)
                tCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTAt                      >  1:9151735/1‑93 (MQ=255)
                tCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTAt                      >  1:157028/1‑93 (MQ=255)
                 ccAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATc                     <  1:8332989/93‑1 (MQ=255)
                  cAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATcc                    <  1:4278846/93‑1 (MQ=255)
                  cAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATcc                    >  1:9502727/1‑93 (MQ=255)
                  cAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATcc                    <  2:8380262/93‑1 (MQ=255)
                  cAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATcc                    <  1:5208175/93‑1 (MQ=255)
                  cAGCCGACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATcc                    <  1:5208176/93‑1 (MQ=255)
                   aGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCg                   >  1:3706310/1‑93 (MQ=255)
                   aGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCg                   <  2:6384151/93‑1 (MQ=255)
                    gcgtACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGa                  <  2:7203631/90‑1 (MQ=255)
                    gCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGa                  >  1:5993750/1‑93 (MQ=255)
                     ccTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGAc                 <  1:358462/93‑1 (MQ=19)
                     ccTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGAc                 >  1:3568107/1‑93 (MQ=19)
                     ccTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGAc                 <  2:2134075/93‑1 (MQ=255)
                     ccTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGAc                 <  1:3434798/93‑1 (MQ=255)
                     ccTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGAc                 <  2:9243079/93‑1 (MQ=19)
                      cTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGAcc                <  2:6044155/93‑1 (MQ=255)
                      cTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGAcc                <  2:6044156/93‑1 (MQ=19)
                       tACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑GGGCGTAAACGCCTTATCCGGCCt               >  1:5317164/1‑93 (MQ=11)
                       tACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCt               >  1:8819849/1‑93 (MQ=19)
                       tACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCt               >  2:463741/1‑93 (MQ=19)
                       tACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCt               >  2:4989104/1‑93 (MQ=19)
                       tACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCt               >  2:8770905/1‑93 (MQ=19)
                       tACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCt               >  2:7563463/1‑93 (MQ=19)
                       tACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCt               >  1:7597774/1‑93 (MQ=19)
                       tACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAAACCCTTATCCGACCt               >  1:5009261/1‑93 (MQ=11)
                        aCGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCtc              <  2:4943739/93‑2 (MQ=255)
                        aCGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTa              <  1:708383/93‑1 (MQ=255)
                         cGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTAc             <  1:6804525/93‑1 (MQ=19)
                         cGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTAc             <  1:6804526/93‑1 (MQ=19)
                         cGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTAc             <  1:1189835/93‑1 (MQ=19)
                         cGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTAc             <  2:5646144/93‑1 (MQ=19)
                          ggTTATGTTCCGGTTTGTATGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTACg            <  1:7819067/93‑1 (MQ=255)
                          ggTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTACg            >  2:8142221/1‑93 (MQ=19)
                           gTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTACgg           >  2:1564069/1‑93 (MQ=19)
                           gTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTACgg           >  1:7336287/1‑93 (MQ=19)
                           gTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTACgg           >  2:6657971/1‑93 (MQ=19)
                           gTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTACgg           >  2:3075556/1‑93 (MQ=19)
                            ttATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTACGGt          >  2:5369356/1‑93 (MQ=19)
                            ttATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTACGGt          >  2:787908/1‑93 (MQ=19)
                            ttATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTACGGt          <  1:5010540/93‑1 (MQ=19)
                            ttATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTACGGt          <  1:9024573/93‑1 (MQ=19)
                            ttATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTACGGt          >  2:8198860/1‑93 (MQ=19)
                            ttATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTACGGt          >  1:3945467/1‑93 (MQ=19)
                             tATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTACGGtt         >  2:1075836/1‑93 (MQ=19)
                             tATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTACGGtt         <  1:7485320/93‑1 (MQ=19)
                              aTGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTACGGTTa        >  2:5744668/1‑93 (MQ=19)
                              aTGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTACGGTTa        >  1:7773646/1‑93 (MQ=19)
                              aTGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTACGGTTa        >  1:4262327/1‑93 (MQ=19)
                              aTGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTACGGTTa        >  1:7773647/1‑93 (MQ=19)
                               tGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTACGGTTAt       >  1:6056150/1‑93 (MQ=19)
                               tGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTACGGTTAt       <  2:638957/93‑1 (MQ=19)
                                gTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCTACCTACGGTTATg      <  2:10027465/93‑1 (MQ=255)
                                gTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTACGGTTATg      <  1:3808630/93‑1 (MQ=19)
                                gTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTACGGTTATg      >  2:9458669/1‑93 (MQ=19)
                                  tCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTACGGTTATGtt    >  1:3236457/1‑93 (MQ=12)
                                   ccGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTACGGTTATGTTc   >  1:918627/1‑93 (MQ=9)
                                   ccGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTACGGTTATGTTc   >  1:2479340/1‑93 (MQ=9)
                                   ccGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTACGGTTATGTTc   <  1:149425/93‑1 (MQ=9)
                                    cGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATG‑CGGCGTAAACGCCTTATCCGACCTACGGTTATGTTcc  <  2:3939273/93‑1 (MQ=0)
                                                                                            |                                     
CGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCC  >  NC_000913/3389719‑3389848

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