Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 266,105 A→G intergenic (‑328/+303) yafY ← / ← yafZ hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913266,1050AG98.8% 774.6 / ‑2.6 244intergenic (‑328/+303)yafY/yafZhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base G (131/110);  total (134/110)
Fisher's exact test for biased strand distribution p-value = 2.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 247 total aligned reads displayed.

ATACTCCAGGTCCAGATACATCAGCCGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACACTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCATGGCGGCGTTCTGCCTGTGTCATGAG  >  NC_000913/266015‑266194
                                                                                          |                                                                                         
aTACTCCAGGTCCAGATACATCAGCCGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCt                                                                                         <  1:3315696/93‑1 (MQ=255)
 tACTCCAGGTCCAGATACATCAGCCGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTc                                                                                        >  2:4485048/1‑93 (MQ=255)
  aCTCCAGGTCCAGATACATCAGCCGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCa                                                                                       <  1:1087605/93‑1 (MQ=255)
  aCTCCAGGTCCAGATACATCAGCCGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCa                                                                                       <  1:1470178/93‑1 (MQ=255)
  aCTCCAGGTCCAGATACATCAGCCGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCa                                                                                       >  2:5957297/1‑93 (MQ=255)
   cTCCAGGTCCAGATACATCAGCCGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACACTCAg                                                                                      >  2:5254060/1‑93 (MQ=255)
     ccAGGTCCAGATACATCAGCCGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCg                                                                                    >  2:5613885/1‑93 (MQ=255)
      cAGGTCCAGATACATCAGCCGCTCGCGAAAATCACGCCGCCGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGt                                                                                   >  1:341246/1‑93 (MQ=255)
      cAGGTCCAGATACATCAGCCGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGt                                                                                   >  2:707639/1‑93 (MQ=255)
        ggTCCAGATACATCAGCCGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGttt                                                                                 <  1:4845515/93‑1 (MQ=255)
        ggTCCAGATACATCAGCCGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGttt                                                                                 >  1:5936024/1‑93 (MQ=255)
            cacaTACATCAGCCGCTCAAGAAACTCACGCCGCAGCGTGCGCACCGACACCCCAACCTCAGCGGCAAGCTTACGCACGCTCAGCGTTTcccc                                                                             <  2:3431592/90‑1 (MQ=18)
            cAGATACATCAGCCGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTcccc                                                                             <  2:2113036/93‑1 (MQ=255)
            cAGATACATCAGCCGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTcccc                                                                             >  2:227504/1‑93 (MQ=255)
             aGATACATCAGCCGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCt                                                                            <  2:5818019/93‑1 (MQ=255)
              gATACATCAGCCGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCtc                                                                           >  1:1857873/1‑92 (MQ=255)
              gATACATCAGCCGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTg                                                                           >  1:6018847/1‑93 (MQ=255)
               aTACATCAGCCGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGc                                                                          >  2:4032096/1‑93 (MQ=255)
                tACATCAGCCGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGcc                                                                         <  1:5368792/93‑1 (MQ=255)
                  cATCAGCCGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGccac                                                                       >  1:4336764/1‑93 (MQ=255)
                  cATCAGCCGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGccac                                                                       <  2:4863211/93‑1 (MQ=255)
                   aTCAGCCGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGccacc                                                                      <  1:33642/93‑1 (MQ=255)
                   aTCAGCCGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGccacc                                                                      >  2:4942103/1‑93 (MQ=255)
                   aTCAGCCGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGccacc                                                                      <  2:4979214/93‑1 (MQ=255)
                    tCAGCCGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACTCCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGccacca                                                                     <  1:4478734/93‑1 (MQ=255)
                    tCAGCCGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGccacca                                                                     >  1:5006767/1‑93 (MQ=255)
                      aGCCGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGa                                                                   >  2:3445925/1‑93 (MQ=255)
                       gCCGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGAc                                                                  >  1:5068211/1‑93 (MQ=255)
                        ccGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACg                                                                 >  2:4793200/1‑93 (MQ=255)
                         cGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACgg                                                                <  2:1171875/93‑1 (MQ=255)
                         cGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACgg                                                                >  1:6492007/1‑93 (MQ=255)
                          gCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCTCTGCGACCAGACGGc                                                               >  1:3459651/1‑93 (MQ=255)
                           cTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCt                                                              >  2:842099/1‑93 (MQ=255)
                            tCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTg                                                             >  2:5528996/1‑93 (MQ=255)
                            tCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTg                                                             >  1:3093983/1‑93 (MQ=255)
                             cACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGa                                                            <  1:3613786/93‑1 (MQ=255)
                             cACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGATTGCCCTGCCACCAGACGGCTGa                                                            >  1:3182424/1‑93 (MQ=255)
                              aCGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGAt                                                           >  2:2996214/1‑93 (MQ=255)
                              aCGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGAt                                                           <  2:6135975/93‑1 (MQ=255)
                                 aaaaTCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTAt                                                        >  2:2135195/1‑93 (MQ=255)
                                 aaaaTCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTAt                                                        >  1:2346710/1‑93 (MQ=255)
                                 aaaaTCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTAt                                                        >  1:6639225/1‑93 (MQ=255)
                                   aaTCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCa                                                      <  1:1566408/93‑1 (MQ=255)
                                    aTCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAg                                                     >  2:246564/1‑93 (MQ=255)
                                    aTCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAg                                                     >  1:3950755/1‑93 (MQ=255)
                                      cACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGtg                                                   <  1:2934933/93‑1 (MQ=255)
                                       aCGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGtga                                                  <  2:4225767/93‑1 (MQ=255)
                                       aCGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGtga                                                  >  2:1763437/1‑93 (MQ=255)
                                       aCGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGtga                                                  >  1:627211/1‑93 (MQ=255)
                                        cgccgcAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTAAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGtgac                                                 <  2:825491/93‑1 (MQ=255)
                                          ccgcAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGtgacag                                               >  2:6369489/1‑93 (MQ=255)
                                            gcAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGtgacagcc                                             <  2:2021812/93‑1 (MQ=255)
                                            gcAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGtgacagcc                                             <  1:3996797/93‑1 (MQ=255)
                                             cAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGtgacagcct                                            <  1:2424023/93‑1 (MQ=255)
                                             cAGCGTGCGCACCGACACACCAAACTAAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGtgacagcct                                            <  1:6728593/93‑1 (MQ=255)
                                              aGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCCCCAGACGGCTGATTATCAGtgacagcctg                                           <  1:5327511/93‑1 (MQ=255)
                                               gCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGtgacagcctga                                          <  1:6335707/93‑1 (MQ=255)
                                               gCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGtgacagcctga                                          >  2:123106/1‑93 (MQ=255)
                                               gCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGtgacagcctga                                          >  2:1306900/1‑93 (MQ=255)
                                                cGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGtgacagcctgac                                         <  2:3606857/93‑1 (MQ=255)
                                                 gTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGtgacagcctgaca                                        <  1:5307053/93‑1 (MQ=255)
                                                   gcgcACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGtgacagcctgacagc                                      <  1:4691403/93‑1 (MQ=255)
                                                    cgcACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGtgacagcctgacagcc                                     <  1:349356/93‑1 (MQ=255)
                                                     gcACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGAcagcca                                    <  1:857042/93‑1 (MQ=255)
                                                      cACCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGAcagccag                                   >  2:6213118/1‑93 (MQ=255)
                                                       cccGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGAcagccagc                                  <  1:5794520/92‑1 (MQ=255)
                                                       aCCGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGAcagcccgc                                  >  1:810437/1‑93 (MQ=255)
                                                        ccGACACACCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACagccagcc                                 <  2:2094187/93‑1 (MQ=255)
                                                            cacacCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTc                             >  1:4144837/1‑93 (MQ=255)
                                                             ccacCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCa                            <  2:880897/92‑1 (MQ=255)
                                                             acacCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCa                            <  2:6003894/93‑1 (MQ=255)
                                                              cacCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCAt                           <  1:5635052/93‑1 (MQ=255)
                                                              cacCAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCAt                           <  1:4650251/93‑1 (MQ=255)
                                                                ccAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCATgg                         >  2:6561368/1‑93 (MQ=255)
                                                                ccAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCATgg                         >  1:1025467/1‑93 (MQ=255)
                                                                 cAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCATggc                        <  2:2771335/93‑1 (MQ=255)
                                                                 cAAACTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCATggc                        <  2:1145234/93‑1 (MQ=255)
                                                                  aaaCTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCATggcg                       >  1:6332627/1‑93 (MQ=255)
                                                                   aaCTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCATggcgg                      >  2:956996/1‑93 (MQ=255)
                                                                   aaCTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCATggcgg                      >  1:2921197/1‑93 (MQ=255)
                                                                   aaCTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCATggcgg                      <  1:3790834/93‑1 (MQ=255)
                                                                    aCTCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCATggcggg                     >  2:4094445/1‑92 (MQ=255)
                                                                      tCAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACATCCTGACAGCCAGCCGGTCATGGCGGCGt                   <  2:3175390/93‑1 (MQ=255)
                                                                      ccAGCGGCAAGCTTGCGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCATGGCGGCGt                   <  2:3307818/92‑1 (MQ=255)
                                                                       cAGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCAGTCATGGCGGCGtt                  <  2:3822701/93‑1 (MQ=255)
                                                                       cAGCGGCAAGCTCACGCACGCTCAGCGTTCCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCATGGCGGCGtt                  <  2:2330016/93‑1 (MQ=255)
                                                                        aGCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCATGGCGGCGGTc                 >  2:6453228/1‑93 (MQ=255)
                                                                         gCGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCATGGCGGCGTTCt                >  2:5483061/1‑93 (MQ=255)
                                                                          cGGCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCATGGCGGCGTTCTg               <  1:1919892/93‑1 (MQ=255)
                                                                           ggCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCATGGCGGCGGTCTGc              >  1:6669948/1‑93 (MQ=255)
                                                                            gCAAGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCATGGCGGCGTTCTGcc             <  1:5537697/93‑1 (MQ=255)
                                                                              aaGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCATGGCGGCGTTCTGCCtg           <  2:2392942/93‑1 (MQ=255)
                                                                              aaGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCATGGCGGCGTTCTGCCtg           <  2:3849047/93‑1 (MQ=255)
                                                                               aGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCATGGCGGCGTTCTGCCtgt          >  1:1675628/1‑93 (MQ=255)
                                                                               aGCTTACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCATGGCGGCGTTCTGCCtgt          >  2:1028701/1‑93 (MQ=255)
                                                                                ggttACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCATGGCGGCGTTCTGCCtgtg         <  1:6410159/91‑1 (MQ=255)
                                                                                  ttACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCATGGCGGCGTTCTGCCTGTGTc       >  2:2918605/1‑93 (MQ=255)
                                                                                   tACGCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCAAGCCGGTCATGGCGGCGTTCTGCCTGTGTCa      >  2:5090119/1‑93 (MQ=255)
                                                                                      gCACGCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCATGGCGGCGTTCTGCCTGTGTCATga   <  2:2841022/93‑1 (MQ=255)
                                                                                       cacgCTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCATGGCGGCGTTCTGCCTGTGTCATgag  >  2:5225125/1‑93 (MQ=255)
                                                                                          |                                                                                         
ATACTCCAGGTCCAGATACATCAGCCGCTCACGAAAATCACGCCGCAGCGTGCGCACCGACACACCAAACTCAGCGGCAAGCTTACGCACACTCAGCGTTTCCCCTGCCACCAGACGGCTGATTATCAGTGACAGCCTGACAGCCAGCCGGTCATGGCGGCGTTCTGCCTGTGTCATGAG  >  NC_000913/266015‑266194

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