Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,897,050 +G coding (277/666 nt) yieH → phosphoenolpyruvate and 6‑phosphogluconate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,897,0461.G98.4% 219.3 / ‑3.8 63coding (273/666 nt)yieHphosphoenolpyruvate and 6‑phosphogluconate phosphatase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (32/30);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCAAAGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑A‑GGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCGG  >  NC_000913/3896909‑3897177
                                                                                                                                           |                                                                                                                                   
ttCAAAGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑A‑g                                                                                                                                    <  2:1048860/139‑1 (MQ=255)
          tAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGa‑a                                                                                                                                      >  2:150692/1‑128 (MQ=255)
          tAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGa‑a                                                                                                                                      <  1:150692/128‑1 (MQ=255)
          tAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAg                                                                                                                           >  2:616590/1‑139 (MQ=255)
              aCTGTACGAAATTATCGTTATTGTTGCCCTCGAACATGGTGTTACGTTAGCGAAAACAGCAGCTGAACACGTTTACCGTGCAGAAGCAGCTCGGCTGTTCGATTCAGAACTGGGAGGCACGGATG‑GGGGCTGGATCGCt                                                                                                                       <  2:739494/139‑1 (MQ=255)
                 gTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTcc                                                                                                                     >  2:1306162/1‑138 (MQ=255)
                 gTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTcc                                                                                                                     <  1:1306162/138‑1 (MQ=255)
                 gTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTcc                                                                                                                     >  1:1334792/1‑138 (MQ=255)
                 gTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTcc                                                                                                                     <  2:1334792/138‑1 (MQ=255)
                  tACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGa                                                                                                                            >  1:1142284/1‑130 (MQ=255)
                  tACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGa                                                                                                                            <  2:1142284/130‑1 (MQ=255)
                  tACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTg                                                                                                                   >  1:1493086/1‑139 (MQ=255)
                  tACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTg                                                                                                                   >  1:92925/1‑139 (MQ=255)
                  tACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTg                                                                                                                   <  2:942277/139‑1 (MQ=255)
                  tACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTg                                                                                                                   >  1:27227/1‑139 (MQ=255)
                  tACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTg                                                                                                                   >  2:1270702/1‑139 (MQ=255)
                          tATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGAt                                                                                                           <  2:27227/139‑1 (MQ=255)
                                 atTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTgcg                                                                                                    <  1:1496012/139‑1 (MQ=255)
                                   tGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGACGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAg                                                                                                               >  1:64266/1‑126 (MQ=255)
                                   tGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAg                                                                                                               <  2:64266/126‑1 (MQ=255)
                                     tttCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCaa                                                                                                >  1:200857/1‑139 (MQ=255)
                                     tttCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑A‑gggggc                                                                                                                               <  1:853169/107‑3 (MQ=255)
                                     tttCCCTTGAACATGGTGATACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑A‑gggggc                                                                                                                               >  2:853169/1‑105 (MQ=255)
                                      ttCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGAAGCGCTCCTGTCAGCGATCACTGCGCCAAt                                                                                               >  2:379192/1‑139 (MQ=255)
                                           ttGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                          <  2:1493086/139‑1 (MQ=255)
                                            tGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtgt                                                                                         >  1:497560/1‑139 (MQ=255)
                                                   ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCt                                                                                  >  2:823663/1‑139 (MQ=255)
                                                    gtgtTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑A‑g                                                                                                                                    <  1:1259913/87‑1 (MQ=255)
                                                    gtgtTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑A‑g                                                                                                                                    >  2:1259913/1‑87 (MQ=255)
                                                         aCGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGc                                                                            <  1:22263/139‑1 (MQ=255)
                                                              aGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTcc                                                                                                                     >  1:1296884/1‑93 (MQ=255)
                                                              aGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTcc                                                                                                                     <  2:1296884/93‑1 (MQ=255)
                                                              aGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCg                                                                                                             >  1:967429/1‑101 (MQ=255)
                                                              aGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCg                                                                                                             <  2:967429/101‑1 (MQ=255)
                                                              aGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACt                                                                                                       <  1:669810/107‑1 (MQ=255)
                                                              aGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACt                                                                                                       >  2:669810/1‑107 (MQ=255)
                                                              aGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCaaa                                                                       <  2:200857/139‑1 (MQ=255)
                                                                cGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCGCTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATa                                                                     >  1:964886/1‑139 (MQ=255)
                                                                     aCAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACaaa                                                                >  2:1242639/1‑139 (MQ=255)
                                                                                 cccGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCa                                                                                                         <  2:28644/84‑1 (MQ=255)
                                                                                 cacGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTGAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCa                                                                                                 <  2:1015848/94‑1 (MQ=255)
                                                                                 cacGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCa                                                                                                         >  1:28644/1‑86 (MQ=255)
                                                                                 cacGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCAGTGCGCCa                                                                                                 >  1:1015848/1‑94 (MQ=255)
                                                                                       tACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGc                                              <  2:497560/139‑1 (MQ=255)
                                                                                           gTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGc                                          <  2:179870/139‑1 (MQ=255)
                                                                                           gTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGc                                          <  2:699339/139‑1 (MQ=255)
                                                                                              cAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGa                                       <  1:616590/139‑1 (MQ=255)
                                                                                                 aaGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAata                                    <  1:1270702/139‑1 (MQ=255)
                                                                                                   gTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATggg                                               >  2:1234349/1‑126 (MQ=255)
                                                                                                   gTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATggg                                               <  1:1234349/126‑1 (MQ=255)
                                                                                                      gCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGAt                                                                                                           <  1:1239614/63‑1 (MQ=255)
                                                                                                      gCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGAt                                                                                                           >  2:1239614/1‑63 (MQ=255)
                                                                                                       cTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACaaa                                                                >  2:774618/1‑105 (MQ=255)
                                                                                                       cTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACaaa                                                                <  1:774618/105‑1 (MQ=255)
                                                                                                          ggCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCact                           <  2:1453966/139‑1 (MQ=255)
                                                                                                           gCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCacta                          >  2:351694/1‑139 (MQ=255)
                                                                                                           gCTGTTCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGAAGCATTCTATGGGCAAGCTGAATATGTTGCacta                          >  1:983589/1‑139 (MQ=255)
                                                                                                                tCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATgcagca                                                        <  1:1070241/104‑1 (MQ=255)
                                                                                                                tCGATTCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATgcagca                                                        >  2:1070241/1‑104 (MQ=255)
                                                                                                                    ttCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                          <  2:374380/66‑1 (MQ=255)
                                                                                                                    ttCAGAACTGGAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                          >  1:374380/1‑66 (MQ=255)
                                                                                                                              gAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCaaa                                                                       >  1:949975/1‑75 (MQ=255)
                                                                                                                              gAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCaaa                                                                       <  2:949975/75‑1 (MQ=255)
                                                                                                                              gAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATgca                                                           <  2:477746/87‑1 (MQ=255)
                                                                                                                              gAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATgca                                                           >  1:477746/1‑87 (MQ=255)
                                                                                                                              gAAGCCATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTc       >  1:4424/1‑139 (MQ=255)
                                                                                                                                   cATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCaa                                            >  2:1263988/1‑97 (MQ=255)
                                                                                                                                   cATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCaa                                            <  1:1263988/97‑1 (MQ=255)
                                                                                                                                   cATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCgg  >  2:743987/1‑139 (MQ=255)
                                                                                                                                   cATCGA‑AGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCgg  <  2:983589/139‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                   
TTCAAAGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGA‑A‑GGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCGG  >  NC_000913/3896909‑3897177

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑