Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,039,960 Δ1 bp intergenic (+42/‑36) yihD → / → yihE hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,039,9590A.99.6% 1137.6 / ‑2.5 243intergenic (+41/‑37)yihD/yihEhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (135/107);  total (136/107)
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.33e-01

Only 100 of 275 total aligned reads displayed.

AGAAGACTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATAATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATGG  >  NC_000913/4039867‑4040047
                                                                                            |                                                                                        
agaagaCTTCAAAACGGCGCTGTTACGCGTTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAAACATCGTGAAATGATa                                                                                          <  2:3812582/93‑1 (MQ=255)
agaagaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATa                                                                                          <  2:3596660/93‑1 (MQ=255)
agaagaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATa                                                                                          <  1:2979823/93‑1 (MQ=255)
 gaagaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATat                                                                                         >  1:5813811/1‑92 (MQ=255)
  aagaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATatc                                                                                        <  2:856883/93‑3 (MQ=255)
  aagaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGGAAAAGCTTGTAAGCGGCGCCACCAAAATCCTCGTGAAATGATatc                                                                                        >  1:2692169/1‑91 (MQ=255)
  aagaCTTCAAAACGGCGCTGTTACGCGCTCGCCGCGTAATTAAAGAGTAAAACCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATatc                                                                                        <  2:1258669/93‑3 (MQ=255)
   agaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATatcc                                                                                       >  2:5081459/1‑90 (MQ=255)
   agaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATatcc                                                                                       <  1:3840242/93‑4 (MQ=255)
    gaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATatcct                                                                                      >  2:3736084/1‑89 (MQ=255)
     aCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCtt                                                                                    >  2:878564/1‑93 (MQ=255)
      cTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTc                                                                                   <  2:2600570/93‑1 (MQ=255)
      cTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTc                                                                                   <  1:980793/93‑1 (MQ=255)
       ttCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCg                                                                                  >  1:1922468/1‑93 (MQ=255)
         cAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTc                                                                                >  1:5357483/1‑93 (MQ=255)
         cAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTc                                                                                >  1:5558709/1‑93 (MQ=255)
          aaaaCGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCa                                                                               >  1:5160589/1‑93 (MQ=255)
             aCGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTc                                                                            >  2:4788599/1‑93 (MQ=255)
               ggCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGt                                                                          <  2:6111204/93‑1 (MQ=255)
                gcgcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTa                                                                         <  1:72599/93‑1 (MQ=255)
                gagcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTa                                                                         <  1:3344335/91‑1 (MQ=255)
                  gcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAAt                                                                       >  1:1543592/1‑93 (MQ=255)
                  gcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAAt                                                                       >  2:5705002/1‑93 (MQ=255)
                   cTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATg                                                                      >  1:5418051/1‑93 (MQ=255)
                    tGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGt                                                                     >  2:5405115/1‑93 (MQ=255)
                     gTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGtt                                                                    >  1:2225341/1‑93 (MQ=255)
                     gTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCCTGAAATGATATCCTTCGTCATTCGTAATGtt                                                                    >  1:5509269/1‑93 (MQ=255)
                       tACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGtttt                                                                  <  2:4818218/93‑1 (MQ=255)
                       tACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGtttt                                                                  <  1:2392000/93‑1 (MQ=255)
                          gcgcTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCg                                                               <  2:3015054/93‑1 (MQ=255)
                            gcTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGga                                                             <  1:2886602/93‑1 (MQ=255)
                            gcTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGga                                                             <  1:1449650/93‑1 (MQ=255)
                             cTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgat                                                            >  2:4046240/1‑93 (MQ=255)
                               cgcgGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgatga                                                          >  2:746516/1‑93 (MQ=255)
                               cgcgGCGTAATTAAAGAGTAAAAGCGTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTACCGgatga                                                          <  2:1870896/93‑1 (MQ=255)
                                 cggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGatgatg                                                        <  1:834573/93‑1 (MQ=255)
                                   gcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATggg                                                      >  2:3122279/1‑93 (MQ=255)
                                    cgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGa                                                     <  1:4240779/93‑1 (MQ=255)
                                    cgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGa                                                     <  2:4696003/93‑1 (MQ=255)
                                     gTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGAt                                                    >  2:3675925/1‑93 (MQ=255)
                                     gTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGAt                                                    <  2:5021227/93‑1 (MQ=255)
                                      tAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATg                                                   >  1:4592019/1‑93 (MQ=255)
                                      tAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATg                                                   >  2:1768892/1‑93 (MQ=255)
                                      tAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATg                                                   >  2:1085295/1‑93 (MQ=255)
                                         ttAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAAt                                                >  1:1440994/1‑93 (MQ=255)
                                         ttAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAAt                                                >  2:2681183/1‑93 (MQ=255)
                                         ttAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAAt                                                >  2:1581520/1‑93 (MQ=255)
                                          tAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATa                                               <  1:5935930/93‑1 (MQ=255)
                                          tAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATa                                               <  2:4714584/93‑1 (MQ=255)
                                           aaaGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATaa                                              >  1:2165068/1‑93 (MQ=255)
                                            aaGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAAc                                             >  2:5129045/1‑93 (MQ=255)
                                              gagTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAg                                           <  2:3274771/93‑1 (MQ=255)
                                              gagTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAg                                           <  1:558150/93‑1 (MQ=255)
                                               agTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgc                                          <  2:1429328/93‑1 (MQ=255)
                                                gTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgcg                                         <  2:4154518/93‑1 (MQ=255)
                                                 tAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgcgc                                        <  1:320464/93‑1 (MQ=255)
                                                 tAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgcgc                                        <  1:5318558/93‑1 (MQ=255)
                                                      gCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTa                                   <  1:5181092/93‑1 (MQ=255)
                                                       cTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGTATAACAGCGCTTTTAc                                  >  2:5498864/1‑93 (MQ=255)
                                                       cTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTAc                                  >  1:3039416/1‑93 (MQ=255)
                                                       cTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTAc                                  >  2:221603/1‑93 (MQ=255)
                                                        ttGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACt                                 >  2:5282983/1‑93 (MQ=255)
                                                         tGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACtt                                <  2:3447356/93‑1 (MQ=255)
                                                         tGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACtt                                >  2:4391625/1‑93 (MQ=255)
                                                          gTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACttt                               >  2:5866791/1‑93 (MQ=255)
                                                          gTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACttt                               <  1:355180/93‑1 (MQ=255)
                                                          gTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACttt                               <  1:4440820/93‑1 (MQ=255)
                                                           tAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTc                              <  2:5962649/93‑1 (MQ=255)
                                                            aaGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTcc                             <  2:1372897/93‑1 (MQ=255)
                                                              gcggcgCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAg                           <  2:6025002/93‑1 (MQ=255)
                                                               cggcgCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGa                          <  2:205620/93‑1 (MQ=255)
                                                               cggcgCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGa                          <  1:1369606/93‑1 (MQ=255)
                                                                ggcgCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAGTAACAGCGCTTTTACTTTCCAGac                         >  2:447335/1‑93 (MQ=255)
                                                                ggcgCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCCGaa                         >  1:656545/1‑92 (MQ=255)
                                                                 gcgcCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGaca                        >  1:1025976/1‑93 (MQ=255)
                                                                 gcgcCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGaca                        >  1:4024010/1‑93 (MQ=255)
                                                                  cgcCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGagac                       <  1:4952953/93‑1 (MQ=255)
                                                                  cgcCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGacac                       <  2:4189043/93‑1 (MQ=255)
                                                                   gcCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACt                      >  1:2736006/1‑93 (MQ=255)
                                                                    ctaccaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTa                     <  1:5854963/91‑1 (MQ=255)
                                                                     caccaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTac                    <  2:5909715/93‑1 (MQ=255)
                                                                     caccaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTac                    >  2:4110630/1‑93 (MQ=255)
                                                                     caccaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTac                    >  1:1205556/1‑93 (MQ=255)
                                                                       ccaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTacac                  >  1:206289/1‑93 (MQ=255)
                                                                        caAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACAcc                 <  1:774010/93‑1 (MQ=255)
                                                                        caAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACAcc                 <  2:1466648/93‑1 (MQ=255)
                                                                        caAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACAcc                 <  2:704488/93‑1 (MQ=255)
                                                                        caAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACAcc                 >  2:3058986/1‑93 (MQ=255)
                                                                            atcatcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGa             >  1:275383/1‑93 (MQ=255)
                                                                            atcatcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGa             <  1:4121954/93‑1 (MQ=255)
                                                                                tcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATAcc         >  1:3795480/1‑93 (MQ=255)
                                                                                tcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATAcc         >  1:1875962/1‑93 (MQ=255)
                                                                                 cGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACca        <  2:189406/93‑1 (MQ=255)
                                                                                  gTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcat       <  2:2332222/93‑1 (MQ=255)
                                                                                   tGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcatc      >  2:3209909/1‑93 (MQ=255)
                                                                                   tGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcatc      <  1:105052/93‑1 (MQ=255)
                                                                                    gAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcatca     <  2:2388980/93‑1 (MQ=255)
                                                                                    gAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcatca     >  2:1169061/1‑93 (MQ=255)
                                                                                     aaaTGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcatcat    >  1:1962544/1‑93 (MQ=255)
                                                                                     aaaTGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcatcat    >  1:1765629/1‑93 (MQ=255)
                                                                                            |                                                                                        
AGAAGACTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATAATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATGG  >  NC_000913/4039867‑4040047

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: