Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,689,817 C→A 15.5% intergenic (+66/‑35) manA → / → ydgA mannose‑6‑phosphate isomerase/DUF945 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,689,8170CA15.5% 267.9 / 28.3 100intergenic (+66/‑35)manA/ydgAmannose‑6‑phosphate isomerase/DUF945 family protein
Reads supporting (aligned to +/- strand):  ref base C (45/39);  new base A (10/6);  total (55/45)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 109 total aligned reads displayed.

GGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCG  >  NC_000913/1689719‑1689915
                                                                                                  |                                                                                                  
ggCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGcc                                                                                                   <  2:279949/100‑1 (MQ=255)
ggCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGcc                                                                                                   <  1:1748371/100‑1 (MQ=255)
  ccGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAg                                                                                                 <  2:410742/100‑1 (MQ=255)
   cGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGt                                                                                                >  2:1856810/1‑100 (MQ=255)
    gTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTc                                                                                               >  2:1024185/1‑100 (MQ=255)
     tttAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTcc                                                                                              <  2:2158893/100‑1 (MQ=255)
         gcgcGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCt                                                                                          <  1:1041008/100‑1 (MQ=255)
          cgcaTGTTTACTATACGCTGTAATATCTTACTGAAAAAATTAACATCTCTTTCTAATCTTATTAACGTCTTATATCTCCTTCTTGCGACCAGTCTGCCt                                                                                          <  2:1380041/99‑1 (MQ=255)
             gtgtTTACAACAAGCTGTAAGAGTTTACTTAAAAAATTAGCATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGAGAGTCCGCCTGAtt                                                                                      <  2:893904/100‑1 (MQ=255)
             gtgtTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGAtt                                                                                      <  1:217642/100‑1 (MQ=255)
             gtgtTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGGCCGCCTGaag                                                                                      <  1:1627869/100‑3 (MQ=255)
             gtgtTTACAACAAGCTGTAAGAGCTTACTGAAAAAAGTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGAtt                                                                                      <  1:1571516/100‑1 (MQ=255)
              tgtTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGGCCGCCTGaagg                                                                                     <  2:2221835/100‑4 (MQ=255)
              tgtTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGGCCGCCTGaagg                                                                                     <  2:2172406/100‑4 (MQ=255)
              tgtTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGGCCGCCTGaagg                                                                                     <  1:1820002/100‑4 (MQ=255)
                tttACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTc                                                                                    <  1:104446/99‑1 (MQ=255)
                tttACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTACGCCTGaaggtg                                                                                   >  1:531963/1‑95 (MQ=255)
                 ttACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAt                                                                                  >  1:1460624/1‑100 (MQ=255)
                  tACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAtt                                                                                 <  1:932324/100‑1 (MQ=255)
                  tACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAtt                                                                                 >  2:396754/1‑100 (MQ=255)
                    caacaaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAATGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAtttt                                                                               >  2:1428430/1‑100 (MQ=255)
                      acaaGCTGTCAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGAGGCCAGTCCGCCTGATTTCATTTTAt                                                                             >  2:1630470/1‑100 (MQ=255)
                      acaaGCTGTAAGAGCTTACTGATTAAATTAACATTTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGGCCGCCTGaaggtgtta                                                                                <  2:1228553/97‑9 (MQ=25)
                       caaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATg                                                                            <  1:1856822/100‑1 (MQ=255)
                         aGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAATCTTATTAAAGGCTTATAACTCCTTCAGTCGGCCAGTTTGCCTGATTTCATTTTATg                                                                            <  2:2179805/98‑1 (MQ=255)
                          gCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGAt                                                                         >  2:1256683/1‑100 (MQ=255)
                            tGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATaa                                                                       >  2:2095324/1‑100 (MQ=255)
                            tGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATaa                                                                       >  1:2182813/1‑100 (MQ=255)
                             gTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAAt                                                                      >  1:209470/1‑100 (MQ=255)
                              tAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATc                                                                     <  1:165305/100‑1 (MQ=255)
                               aaGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCa                                                                    <  1:316362/100‑1 (MQ=255)
                               aaGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCa                                                                    <  1:1838356/100‑1 (MQ=255)
                                   gCTGACTGAAAAAATTAACAAAAAATACTAAGCATATTAAAGACTTATAACACCTTCAGACAAACAGTCCGCATAATTTCATTTTATGGAAAATCAtaaa                                                                >  2:260895/1‑97 (MQ=255)
                                    cttacttaAAAAATTAACATCTCTTGCTAAGCTTATTAAATTCTTATAACACCTTCATGCGGCCAGTCTGCCTGTTTTCATTTTATGGATAATAATTATg                                                               <  2:557690/100‑1 (MQ=255)
                                       aCTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAAt                                                            >  1:837255/1‑100 (MQ=255)
                                        cTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCATGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATa                                                           <  2:1948685/100‑1 (MQ=255)
                                        cTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATa                                                           >  1:978539/1‑100 (MQ=255)
                                          gAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATaaa                                                         >  1:349809/1‑100 (MQ=255)
                                          gAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGCCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATaaa                                                         >  2:792139/1‑100 (MQ=255)
                                          gAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTAAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAAt                                                            >  2:741451/1‑97 (MQ=255)
                                          gAAAAAATTAACATCTCTTGCTAAGCTCATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATTGCTAATCATTATGAATa                                                           >  1:269294/1‑98 (MQ=255)
                                            aaaaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACCCCTTCAGGCGGCCAGACCGCCTGATTTAATTTTATGTATAATCATTATGAATaaa                                                         >  2:1370633/1‑98 (MQ=255)
                                            aaaaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATc                                                       <  2:459610/100‑1 (MQ=255)
                                              aaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGc                                                     <  1:1563201/100‑1 (MQ=255)
                                              aaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGc                                                     <  1:367304/100‑1 (MQ=255)
                                              aaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGc                                                     <  1:1222591/100‑1 (MQ=255)
                                               aaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCt                                                    <  1:741451/100‑1 (MQ=255)
                                                  tAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGt                                                 >  2:1490299/1‑100 (MQ=255)
                                                  tAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGt                                                 >  2:1490093/1‑100 (MQ=255)
                                                    aCATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAg                                               >  2:943281/1‑100 (MQ=255)
                                                     cATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGTTTTCATTTTATGGATAATCATTATGAATGAATCGCTGGTacc                                              <  1:1514348/100‑3 (MQ=255)
                                                       tctcTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCgg                                            >  2:422949/1‑100 (MQ=255)
                                                          ctTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAg                                         <  1:43456/100‑1 (MQ=255)
                                                          ctTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAg                                         >  2:1797962/1‑100 (MQ=255)
                                                           ttGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAgg                                        <  2:869621/100‑1 (MQ=255)
                                                            tGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCAGGTAGCGGt                                           >  1:1303411/1‑96 (MQ=255)
                                                              cTAAGCTTATGAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGt                                     <  2:871802/100‑1 (MQ=255)
                                                                  gCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGACTAAATCGCTGGTAGCGGTAGGCGTCAtt                                 >  2:2000287/1‑100 (MQ=255)
                                                                  gCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAtt                                 >  2:740004/1‑100 (MQ=255)
                                                                  gCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAtt                                 >  2:112753/1‑100 (MQ=255)
                                                                  gCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAtt                                 >  2:1289659/1‑100 (MQ=255)
                                                                     tattaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAttgtt                              >  1:590416/1‑100 (MQ=255)
                                                                      attaAAAGATTATAACAACATCAGCCGGCCAGTCCGCCTGATCTCATTTTATGGATACTCATTATGCATCACTCGCTGGTAGCGGTAGGCGTCActgttg                             <  2:449023/100‑1 (MQ=255)
                                                                       ttaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTgc                            <  2:2205284/100‑1 (MQ=255)
                                                                         aaaGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTgcgc                          >  1:1282646/1‑100 (MQ=255)
                                                                         aaaGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTgcgc                          <  1:2235702/100‑1 (MQ=255)
                                                                         aaaGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTgcgc                          <  1:408911/100‑1 (MQ=255)
                                                                         aaaGGCTTATAACACCTTCAGGCGACCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTgcg                           >  1:143127/1‑99 (MQ=255)
                                                                          aaGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCt                         <  1:1884289/100‑1 (MQ=255)
                                                                            ggCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAg                       >  2:612219/1‑100 (MQ=255)
                                                                            ggCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAg                       >  2:444444/1‑100 (MQ=255)
                                                                            ggCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAg                       <  2:1766154/100‑1 (MQ=255)
                                                                              cTTATAACACCTTCAGGCGACCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTACAGGTAGGCGTCATTGTTGCGATAGGc                     >  2:434948/1‑100 (MQ=255)
                                                                               ttATAACACCTTCAGGCGGCCAGTCCGACTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGc                     <  2:412257/99‑1 (MQ=255)
                                                                                tataACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGt                   >  2:1985500/1‑100 (MQ=255)
                                                                                tataACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGt                   >  1:241726/1‑100 (MQ=255)
                                                                                  aaatgaaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGTCGTCATTGTTGCGCTAGGc                     >  2:1818965/9‑96 (MQ=37)
                                                                                  aaatgaaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAg                 >  1:834810/9‑100 (MQ=38)
                                                                                   aaCACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCACTAGGCGTAGt                <  1:1490103/100‑1 (MQ=255)
                                                                                    atgaaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTc               >  2:1137308/7‑100 (MQ=38)
                                                                                    acacCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTc               >  2:1498000/1‑100 (MQ=255)
                                                                                     tgaaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTTCGCTAGGCGTAGTCt              >  2:1471394/6‑100 (MQ=38)
                                                                                      gaaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTg             >  2:598504/5‑100 (MQ=255)
                                                                                      gaaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGAGTCATTGTTACGATAGGCGTAGTCTg             >  2:1971282/5‑100 (MQ=255)
                                                                                       tcTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTgg            <  1:834989/99‑1 (MQ=255)
                                                                                       ccTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTgg            <  1:1788632/100‑1 (MQ=255)
                                                                                       ccTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTgg            >  2:685174/1‑100 (MQ=255)
                                                                                       aaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTg             >  1:292182/4‑99 (MQ=255)
                                                                                        cTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGa           >  1:1826768/1‑100 (MQ=255)
                                                                                         ttCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGAc          <  2:499699/100‑1 (MQ=255)
                                                                                         ttCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGAc          <  1:1985514/100‑1 (MQ=255)
                                                                                         ttCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTACGCTAGGCGTAGTATGGAc          >  2:2060541/1‑100 (MQ=255)
                                                                                          tCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCTGTAGACGTCATTGTTGCGATAGATGTAGTTTgg            >  2:937435/1‑97 (MQ=255)
                                                                                          tCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGaga         >  2:834989/1‑98 (MQ=255)
                                                                                          tCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACa         >  1:1344675/1‑100 (MQ=255)
                                                                                           cAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCCTTGTTGCGCTAGGCGTAGTCTGGACAg        >  2:162478/1‑100 (MQ=255)
                                                                                             ggcggcCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAggc      >  2:285225/1‑100 (MQ=255)
                                                                                                ggcCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAggcggc   <  2:1617327/100‑1 (MQ=255)
                                                                                                ggcCAATCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAggcggc   >  2:1072085/1‑100 (MQ=255)
                                                                                                 gcCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGgcggcg  <  2:2182798/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
GGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCG  >  NC_000913/1689719‑1689915

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: