Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,347,069 G→A *164* (TAG→TAA ptsN → sugar‑specific enzyme IIA component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,347,0690GA98.8% 255.6 / NA 80*164* (TAG→TAAptsNsugar‑specific enzyme IIA component of PTS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (39/40);  minor base T (1/0);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTGCATACTCTGTCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAGTTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGATATGGGTTTTTACTGCGTGGA  >  NC_000913/3346936‑3347203
                                                                                                                                     |                                                                                                                                      
aCCTGCATACTCTGTCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTa                                                                                                                                     <  2:1274926/137‑1 (MQ=255)
aCCTGCATACTCTGTCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTa                                                                                                                                     >  1:1274926/1‑137 (MQ=255)
         ctctGTCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGt                                                                                                                          <  1:1136/139‑1 (MQ=255)
            tGTCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGtctc                                                                                                                       >  1:771337/1‑139 (MQ=255)
            tGTCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGtctc                                                                                                                       >  1:1287057/1‑139 (MQ=255)
              tCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCtt                                                                                                                     >  1:46309/1‑139 (MQ=255)
              tCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCtt                                                                                                                     >  1:349246/1‑139 (MQ=255)
               cGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCttt                                                                                                                    >  2:968780/1‑139 (MQ=255)
                      ggCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGt                                                                                                             <  2:771337/139‑1 (MQ=255)
                      ggCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGt                                                                                                             >  2:1334954/1‑139 (MQ=255)
                      ggCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGt                                                                                                             <  2:1287057/139‑1 (MQ=255)
                        cGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTtg                                                                                                           >  1:1055332/1‑139 (MQ=255)
                        cGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTtg                                                                                                           <  2:465052/139‑1 (MQ=255)
                                  ggCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGaaa                                                                                                 <  1:410612/139‑1 (MQ=255)
                                    cGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACa                                                                                               <  2:700147/139‑1 (MQ=255)
                                       aCAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTa                                                                                            <  2:1505401/139‑1 (MQ=255)
                                       aCAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTa                                                                                            <  1:811374/139‑1 (MQ=255)
                                         aaaaCCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAg                                                                                                                 <  2:1227522/116‑1 (MQ=255)
                                         aaaaCCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAg                                                                                                                 >  1:1227522/1‑116 (MQ=255)
                                         aaaaCCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACa                                                                                          <  1:1334954/139‑1 (MQ=255)
                                         aaaaCCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACa                                                                                          <  2:46309/139‑1 (MQ=255)
                                         aaaaCCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACa                                                                                          <  1:1176009/139‑1 (MQ=255)
                                         aaaaCCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACa                                                                                          <  1:968780/139‑1 (MQ=255)
                                                 cTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGaa                                                                                                  <  1:331905/123‑1 (MQ=255)
                                                 cTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGaa                                                                                                  >  2:331905/1‑123 (MQ=255)
                                                       ccgTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATc                                                                            >  1:341034/1‑139 (MQ=255)
                                                       ccgTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATc                                                                            >  2:1121733/1‑139 (MQ=255)
                                                       ccgTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATc                                                                            >  1:640342/1‑139 (MQ=255)
                                                              cgcgCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGt                                                                                                                          >  2:1376856/1‑86 (MQ=255)
                                                              cgcgCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGt                                                                                                                          <  1:1376856/86‑1 (MQ=255)
                                                                cgcAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGa                                                                   <  2:341034/139‑1 (MQ=255)
                                                                cgcAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGa                                                                   <  2:640342/139‑1 (MQ=255)
                                                                     cccAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGt                                                                                                                          >  1:475725/1‑79 (MQ=255)
                                                                     cccAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGt                                                                                                                          <  2:475725/79‑1 (MQ=255)
                                                                     cccAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACgttc                                                              >  1:1210438/1‑139 (MQ=255)
                                                                               tGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGtct                                                                                                                        <  1:458532/71‑1 (MQ=255)
                                                                               tGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGtct                                                                                                                        >  2:458532/1‑71 (MQ=255)
                                                                               tGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGt                                                                                                                          <  1:1382841/69‑1 (MQ=255)
                                                                               tGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGt                                                                                                                          >  2:1382841/1‑69 (MQ=255)
                                                                               tGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACa                                                                                               <  1:837534/96‑1 (MQ=255)
                                                                               tGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACa                                                                                               >  2:837534/1‑96 (MQ=255)
                                                                               tGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTtcaggttcaggt                                                    >  1:695480/1‑139 (MQ=255)
                                                                                gAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTa                                                   >  2:830722/1‑139 (MQ=255)
                                                                                    agCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGgaga                                                                                                   >  2:67706/1‑87 (MQ=255)
                                                                                    agCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGgaga                                                                                                   <  1:67706/87‑1 (MQ=255)
                                                                                    agCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACa                                                                                          >  2:1064081/1‑96 (MQ=255)
                                                                                    agCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACa                                                                                          <  1:1064081/96‑1 (MQ=255)
                                                                                         tATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACgttcaggt                                                          >  1:1314087/1‑123 (MQ=255)
                                                                                         tATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACgttcaggt                                                          <  2:1314087/123‑1 (MQ=255)
                                                                                         tATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCg                                          >  2:719943/1‑139 (MQ=255)
                                                                                               aTCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCtgc                                    <  1:830722/139‑1 (MQ=255)
                                                                                                   ttACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGt                                                                                                                          >  2:376753/1‑49 (MQ=255)
                                                                                                   ttACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGt                                                                                                                          <  1:376753/49‑1 (MQ=255)
                                                                                                   ttACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTctct                                                                                                                      <  1:1288003/53‑1 (MQ=255)
                                                                                                   ttACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTctct                                                                                                                      >  2:1288003/1‑53 (MQ=255)
                                                                                                   ttACGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATgg                                                                                       >  1:1014464/1‑84 (MQ=255)
                                                                                                   ttACGGATACCGAAGGTACTCCGGATGAAGCGTAATTACTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATgg                                                                                       <  2:1014464/84‑1 (MQ=255)
                                                                                                     aCGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTAc                                                                                    >  1:1273616/1‑85 (MQ=255)
                                                                                                     aCGGATACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTAc                                                                                    <  2:1273616/85‑1 (MQ=255)
                                                                                                         aTACCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGa                          >  1:472848/1‑139 (MQ=255)
                                                                                                           aCCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTtgtg                                                                                                         <  2:526142/58‑1 (MQ=255)
                                                                                                           aCCGAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTtgtg                                                                                                         >  1:526142/1‑58 (MQ=255)
                                                                                                              gAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGata                     <  2:1055332/139‑1 (MQ=255)
                                                                                                              gAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGata                     <  1:1121733/139‑1 (MQ=255)
                                                                                                              gAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGata                     <  2:472848/139‑1 (MQ=255)
                                                                                                              gAAGGTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGata                     <  2:695480/139‑1 (MQ=255)
                                                                                                                 ggTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGATATgg                  >  1:520241/1‑139 (MQ=255)
                                                                                                                  gTACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGATATggg                 >  2:651877/1‑139 (MQ=255)
                                                                                                                   tACTCCGGATGAAGCGTATTTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTa                                                                                            >  1:641347/1‑63 (MQ=255)
                                                                                                                   tACTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTa                                                                                            <  2:641347/63‑1 (MQ=255)
                                                                                                                     cTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGata                     >  2:158991/1‑132 (MQ=255)
                                                                                                                     cTCCGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGata                     <  1:158991/132‑1 (MQ=255)
                                                                                                                       ccGGATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGATATGGGttttt            >  2:856128/1‑139 (MQ=255)
                                                                                                                          gATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTtg                                                                                                           <  1:137070/41‑1 (MQ=255)
                                                                                                                          gATGAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTtg                                                                                                           >  2:137070/1‑41 (MQ=255)
                                                                                                                             gAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGaaa                                                                                                 <  2:1100746/48‑1 (MQ=255)
                                                                                                                             gAAGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGaaa                                                                                                 >  1:1100746/1‑48 (MQ=255)
                                                                                                                               aGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGaaa                                                                                                 <  1:1412034/46‑1 (MQ=255)
                                                                                                                               aGCGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGaaa                                                                                                 >  2:1412034/1‑46 (MQ=255)
                                                                                                                                 cGTAATTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGATATGGGTTTTTACTGCGTGGa  >  1:93092/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                      
ACCTGCATACTCTGTCGCTGGTGGCGAAACGTCTGGCGGACAAAACCATCTGCCGCCGTTTACGCGCAGCCCAGAGCGATGAAGAGCTGTATCAAATCATTACGGATACCGAAGGTACTCCGGATGAAGCGTAGTTATTCGGTAATGTCTCTTTTAGACGTTGTGAGGAGAAACAGTACATGGTACTGATGATCGTCAGCGGACGTTCAGGTTCAGGTAAATCTGTCGCCCTGCGTGCGCTGGAAGATATGGGTTTTTACTGCGTGGA  >  NC_000913/3346936‑3347203

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