Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 705,972 (T)9→10 intergenic (+82/‑121) nagE → / → glnS N‑acetyl glucosamine specific PTS enzyme IIC, IIB, and IIA components/glutamyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913705,9631.T93.4% 293.5 / 12.4 91intergenic (+73/‑130)nagE/glnSN‑acetyl glucosamine specific PTS enzyme IIC, IIB, and IIA components/glutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base . (3/3);  new base T (48/37);  total (51/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

Only 100 of 105 total aligned reads displayed.

GCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑TTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCA  >  NC_000913/705869‑706053
                                                                                               |                                                                                              
gCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑ttttt                                                                                         >  1:2174574/1‑100 (MQ=255)
 cTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑tttttt                                                                                       <  2:286261/100‑1 (MQ=255)
    tATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑ttttttttt                                                                                   <  1:2693698/100‑1 (MQ=255)
     aTGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑tttttttttt                                                                                  <  2:656472/100‑2 (MQ=255)
      tGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑ttttttttttc                                                                                 <  2:2337130/100‑3 (MQ=255)
       gAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑tttttttttttc                                                                                >  2:494790/1‑97 (MQ=255)
       gAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑ttttttttttca                                                                                <  2:789799/100‑4 (MQ=255)
        aaaTCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑ttttttttttcat                                                                               <  2:2692987/100‑5 (MQ=255)
        aaaTCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑ttttttttttcat                                                                               >  1:2340945/1‑96 (MQ=255)
         aaTCAAAAAGTAATCTGCTTTATGCCTGATGTGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑ttttttttttcatc                                                                              <  2:1812954/100‑6 (MQ=255)
         aaTCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑ttttttttttcatc                                                                              >  1:2244480/1‑95 (MQ=255)
         aaTCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑ttttttttttcatc                                                                              <  1:1332969/100‑6 (MQ=255)
          aTCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTtcatca                                                                              >  1:41699/1‑100 (MQ=255)
           tCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACTAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCatcaa                                                                             <  2:493131/100‑1 (MQ=255)
            cAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGTATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCatcaat                                                                            >  1:1338207/1‑100 (MQ=255)
            cAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCatcaat                                                                            >  1:1433914/1‑100 (MQ=255)
            cAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑TTTTTTTTttat                                                                                >  2:2506494/1‑95 (MQ=255)
             aaaaaGTAATCTGCTTTATGCCTGATTCGTCGTTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCatcaatc                                                                           <  1:997560/100‑1 (MQ=255)
             aaaaaGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCatcaatc                                                                           >  1:1204279/1‑100 (MQ=255)
             aaaaaGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCatcaatc                                                                           >  1:2632692/1‑100 (MQ=255)
              aaaaGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCttcaatca                                                                          >  2:2052840/1‑100 (MQ=255)
              aaaaGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑tttttttttttttttat                                                                           >  2:2283732/1‑90 (MQ=255)
                aaGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCTTCAatcatc                                                                        >  2:1040552/1‑100 (MQ=255)
                 aGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATcc                                                                       <  2:1594087/100‑1 (MQ=255)
                      tCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATa                                                                  >  1:1003448/1‑100 (MQ=255)
                       cTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATaa                                                                 >  1:422665/1‑100 (MQ=255)
                         gCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTTTTCTATTATCCCCATAAt                                                                >  2:2486448/1‑99 (MQ=255)
                           tttATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCt                                                             <  2:310639/100‑1 (MQ=255)
                             taggtatGATGGGACGCTGGAGAGTTGAATCCAAAAGTGACAAGCGGTGGAGATCTTCTCCGCCGCTTTTTTTTTTCATCAATCATCCCCATCATCCTTg                                                           <  2:1206021/94‑1 (MQ=255)
                             tATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGTCAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTg                                                           <  2:583118/100‑1 (MQ=255)
                             tATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTg                                                           <  2:1384964/100‑1 (MQ=255)
                              aggCCAGATACGACGCTTGAACATCGCATAAAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTg                                                           <  2:732604/97‑1 (MQ=255)
                              aTGCCTGATGCGGCGCGTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGt                                                          <  2:2527480/100‑1 (MQ=255)
                              aTGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGt                                                          <  2:1102115/100‑1 (MQ=255)
                               tGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTTCATCAATCATCCCCATAATCCtt                                                            >  2:2467372/1‑98 (MQ=255)
                               tGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGtt                                                         <  2:2174561/100‑1 (MQ=255)
                                 ccTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCCTAATCCTTGtt                                                         >  2:1483411/1‑98 (MQ=255)
                                 ccTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAg                                                       >  1:168274/1‑100 (MQ=255)
                                 ccTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑TTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGa                                                      >  2:165443/1‑100 (MQ=255)
                                  cTGGTGCGACGATTGAAAATCGTATCCAACAATGACCAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTa                                                        <  1:1290902/98‑1 (MQ=255)
                                  cTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGa                                                      <  2:2083498/100‑1 (MQ=255)
                                  cTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGa                                                      >  2:2007263/1‑100 (MQ=255)
                                  cTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGa                                                      <  2:2503638/100‑1 (MQ=255)
                                  cTGATGAGACGCTTGAGCGTCGCATCCAACAATGACTAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGa                                                      <  2:1051654/100‑1 (MQ=255)
                                   ggATGCGACGCTCAAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAAACATCCCCAAAATCCTTTTTa                                                        >  1:1337206/2‑97 (MQ=37)
                                    gATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCCATCATCCCCATAATCCTTGTTAGCtt                                                    >  1:238298/1‑100 (MQ=255)
                                    gATGCGACGCTCAAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGAt                                                     >  1:1489820/1‑100 (MQ=255)
                                    gATGCGACGCTCAAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTTATTTATCATCCCCATAATCCTTGTTAGAtt                                                    >  2:2828966/1‑100 (MQ=255)
                                    gATGCGACGCTCAAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGAtt                                                    >  1:1865674/1‑100 (MQ=255)
                                    gATGCGACGCTCAAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGAtt                                                    >  2:874211/1‑100 (MQ=255)
                                     aTGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTTATCAATCATCCCCATAATCCTTGTTAGAtt                                                    >  2:2491325/1‑99 (MQ=255)
                                     aTGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTa                                                   <  2:490047/100‑1 (MQ=255)
                                     aTGCGACGCTCAAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGAtt                                                    >  1:270742/1‑100 (MQ=255)
                                     aTGCGACGCTCAAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGAtt                                                    >  2:888367/1‑99 (MQ=255)
                                     aTGCGACGCTCAAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTa                                                   >  2:2581003/1‑100 (MQ=255)
                                      tGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTAATTAATTATTCTTTTTATCCTTTTTAGATTa                                                   >  2:620995/1‑99 (MQ=255)
                                       gCGACGCTTGAGCGTCGCATCAAAAAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTa                                                   <  2:2268445/98‑1 (MQ=255)
                                        gggcGCTTGGGCGTTGAATACAAAATTGATTGGCGGTGGGGATCTTCTCTGCGGTTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCa                                                <  2:2412766/97‑1 (MQ=255)
                                        cGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCa                                                >  1:377022/1‑100 (MQ=255)
                                          actattaaGCGTCTCATCCAACAATGTCAAGCGTCGGAGATCTTCTCTGCCGC‑‑TTTTTGTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCa                                                <  2:988874/91‑1 (MQ=255)
                                             cTTGAGCGTCGAATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑TTTTTTTTATCATCAATCATCCCCATAATCCTTGTTAGATTATCAAtttt                                           <  2:1577733/100‑1 (MQ=255)
                                              ttGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTTTTAGATTATCAATTTTa                                          >  1:1864338/1‑100 (MQ=255)
                                              ttGAGCCTCGCATCCAAAAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTa                                          <  2:2007655/100‑1 (MQ=255)
                                               tGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTaa                                         >  1:937948/1‑100 (MQ=255)
                                               tGAGCGTCGCATCCAAAAATGACGAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTaa                                         <  1:1040554/100‑1 (MQ=255)
                                                gAGCGTCGCATCCAACAATGACAATCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTaaa                                        <  2:1341631/100‑1 (MQ=255)
                                                gAGCGTCGCATCCAACAATGACAAGCAGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTaaa                                        <  2:1341513/100‑1 (MQ=255)
                                                    gTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAAc                                    >  1:339674/1‑100 (MQ=255)
                                                     tCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTTTTTTATTTTTTCCTTAATCCTTGTTTGATTATCAATTTTAAAAAACt                                   >  2:653535/1‑100 (MQ=255)
                                                        cATCCAACAATGACGAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAAc                                <  1:2491341/100‑1 (MQ=255)
                                                         aTCCAACAATGAAAAGCGGTAGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAAc                                <  2:1061299/99‑1 (MQ=255)
                                                          tCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAg                              >  1:1782064/1‑100 (MQ=255)
                                                           ccAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGt                             >  1:1238608/1‑100 (MQ=255)
                                                           ccAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGt                             >  1:1442507/1‑100 (MQ=255)
                                                           caaaaaaTGACAAGCGGTGGAGATCTTCTCTCCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAAc                                <  2:1864325/92‑1 (MQ=255)
                                                            caacaaTGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGtt                            >  1:1542631/1‑100 (MQ=255)
                                                             aacaaTGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCACTCATCCCCATAATCCTTGTTAGATTATCACTTTTAAAAAACTAACAGtt                            >  1:1686446/1‑99 (MQ=255)
                                                             aacaaTGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTg                           >  1:1556255/1‑100 (MQ=255)
                                                              acaaTGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGt                          >  1:2123311/1‑100 (MQ=255)
                                                              acaaAGCCAAGCGGTGGAGATTTTCTTTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAGAAACTAACAGtt                            <  2:1063403/98‑1 (MQ=255)
                                                                aaTGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCa                        >  1:1699419/1‑100 (MQ=255)
                                                                  tGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGc                      >  1:2761832/1‑100 (MQ=255)
                                                                   gACAAGCGGTGGAAATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTc                         <  2:1338202/96‑1 (MQ=255)
                                                                    aCAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGcc                     <  2:1436878/99‑1 (MQ=255)
                                                                     cAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTg                   >  1:1609425/1‑100 (MQ=255)
                                                                         cGGTGGAGATCTACTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTcc                <  1:757955/99‑1 (MQ=255)
                                                                           gTGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGc             <  2:2841271/100‑1 (MQ=255)
                                                                            tGGAGATCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCt            >  1:553866/1‑100 (MQ=255)
                                                                               agaTCTTCTCTGCCGCTTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCtt           <  2:1489813/99‑1 (MQ=255)
                                                                               agaTCTTCTCTGCCGCTTTTTTTTTTTATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTACCGCTtat         >  2:1226449/1‑100 (MQ=255)
                                                                                gaTCTTCTCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTtata        >  1:1447388/1‑100 (MQ=255)
                                                                                  tcttctCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAg      >  1:550426/1‑100 (MQ=255)
                                                                                   cttctCTGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGa     >  1:1656290/1‑100 (MQ=255)
                                                                                    ttctCTGCCGC‑‑TTTTTTTTTCATCAATAATACCCATAATCCTTGTTATCTTATCAATTTTAAAAAACTAACAGTTGTTAGTCTGTCACTCTTATAAGAt    >  2:716930/1‑99 (MQ=255)
                                                                                    tccccggCCGC‑‑TTATTTTATTCATCAATCATCCCCATAATCCTTGTTAGATTATCAAGTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGAt    <  1:588771/94‑1 (MQ=255)
                                                                                      ctctGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGAt    <  1:1500588/98‑1 (MQ=255)
                                                                                      ctctGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATc   <  1:2856622/99‑1 (MQ=255)
                                                                                      ctctGCCGCTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCa  <  2:1821603/100‑1 (MQ=255)
                                                                                       tctGCCGCTTTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGAt    <  2:1820953/98‑1 (MQ=255)
                                                                                        ccgCCGCTTTTTTTTTTCATCAATCCACCCCATCATCATTGTTAGATTATCAATCTTAAAAAACTAACAGTTGTCAGCCGGTCCCGCTTATAAGATc   <  2:238297/95‑1 (MQ=255)
                                                                                               |                                                                                              
GCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑‑TTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCA  >  NC_000913/705869‑706053

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