Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,795,591 G→A S221S (TCG→TCA waaF → ADP‑heptose:LPS heptosyltransferase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,795,5910GA99.0% 338.0 / ‑4.7 97S221S (TCG→TCAwaaFADP‑heptose:LPS heptosyltransferase II
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (58/38);  total (59/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GGAGTTTGGTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGAAACACAGCTTGAT  >  NC_000913/3795495‑3795687
                                                                                                |                                                                                                
ggAGTTTGGTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCg                                                                                               <  1:1064418/100‑1 (MQ=255)
ggAGTTTGGTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCg                                                                                               >  1:1221523/1‑100 (MQ=255)
   gTTTGGTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGaaa                                                                                            >  1:753884/1‑100 (MQ=255)
    tttGGTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAg                                                                                           <  2:1674714/100‑1 (MQ=255)
       ggTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATc                                                                                        <  2:812548/100‑1 (MQ=255)
       ggTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGAAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATc                                                                                        >  2:702763/1‑100 (MQ=255)
        gTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCa                                                                                       >  1:1119747/1‑100 (MQ=255)
         tCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCAt                                                                                      <  2:1609/100‑1 (MQ=255)
          ccGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATg                                                                                     >  1:1546622/1‑100 (MQ=255)
              cAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGc                                                                                 <  2:210093/100‑1 (MQ=255)
               cgatctcTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCg                                                                                <  2:271699/94‑1 (MQ=255)
               aaaaCGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCg                                                                                >  2:801902/1‑100 (MQ=255)
               aaaaCGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCg                                                                                >  1:1162393/1‑100 (MQ=255)
                aaaCGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCgg                                                                               <  2:1119742/100‑1 (MQ=255)
                  aCGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGc                                                                             >  2:20616/1‑100 (MQ=255)
                  aCGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGc                                                                             <  2:1992802/100‑1 (MQ=255)
                     cTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCgagatcg                                                                          >  1:271700/1‑94 (MQ=255)
                      tGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATg                                                                         <  2:1162387/100‑1 (MQ=255)
                      tGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATg                                                                         <  2:340106/100‑1 (MQ=255)
                       ggCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATga                                                                        >  2:1469759/1‑100 (MQ=255)
                        gCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATgag                                                                       <  2:427016/100‑1 (MQ=255)
                          cacaCTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGAt                                                                     >  2:1576377/1‑100 (MQ=255)
                          cacaCTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGAt                                                                     >  2:1576448/1‑100 (MQ=255)
                           tcacTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGAtt                                                                    <  1:482262/99‑1 (MQ=255)
                            cacTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGAtta                                                                   >  2:482261/1‑99 (MQ=255)
                            cacTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTc                                                                   <  1:1477026/100‑1 (MQ=255)
                            cacTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTc                                                                   >  2:10522/1‑100 (MQ=255)
                            cacTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTc                                                                   >  2:1944663/1‑100 (MQ=255)
                            cacTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTAGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTc                                                                   >  1:1436083/1‑100 (MQ=255)
                               tACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTg                                                                <  2:281831/100‑1 (MQ=255)
                                 ccACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGcc                                                              >  1:1471381/1‑100 (MQ=255)
                                   aCTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCg                                                             >  2:1727715/1‑99 (MQ=255)
                                   aCTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGc                                                            >  2:1738268/1‑100 (MQ=255)
                                   aCTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGc                                                            >  2:1157922/1‑100 (MQ=255)
                                   aCTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGc                                                            >  1:1817354/1‑100 (MQ=255)
                                   aCTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGc                                                            >  2:1083400/1‑100 (MQ=255)
                                    gatctcGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCt                                                           <  2:1499731/95‑1 (MQ=255)
                                     tATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCtt                                                          <  2:527956/100‑1 (MQ=255)
                                        gCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGa                                                       >  2:564225/1‑100 (MQ=255)
                                         cGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCtagatc                                                      >  1:1499740/1‑95 (MQ=255)
                                         cGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGaa                                                      >  2:1711851/1‑100 (MQ=255)
                                          ggAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGaa                                                      <  1:692884/99‑1 (MQ=255)
                                           gcgCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTATGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATa                                                    >  1:1382755/3‑100 (MQ=255)
                                             gCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATAcc                                                  <  2:299051/100‑1 (MQ=255)
                                              cTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCg                                                 <  1:670477/100‑1 (MQ=255)
                                               tGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGa                                                >  1:1329210/1‑100 (MQ=255)
                                               tGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGa                                                >  2:1855911/1‑100 (MQ=255)
                                               tGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGa                                                >  2:1372257/1‑100 (MQ=255)
                                                ggCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGag                                               >  2:153169/1‑100 (MQ=255)
                                                ggCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGag                                               >  1:49926/1‑100 (MQ=255)
                                                 ccgatctagcTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGagc                                              <  2:1890484/93‑1 (MQ=255)
                                                  cAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGagca                                             >  2:440041/1‑100 (MQ=255)
                                                  cAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTGTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCGTTGAATACCGagca                                             >  2:440119/1‑100 (MQ=255)
                                                   aaaGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGagcag                                            <  1:1813780/100‑1 (MQ=255)
                                                    aaGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGagcagn                                           <  1:246597/100‑2 (MQ=255)
                                                    aaGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGagcagc                                           <  2:753881/100‑1 (MQ=255)
                                                        agcTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGagcagatcgg                                       >  1:1890502/1‑95 (MQ=255)
                                                        agcTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGc                                       >  1:100660/1‑100 (MQ=255)
                                                         gcTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGc                                       <  1:940126/99‑1 (MQ=255)
                                                          cTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTAAGCGAAAGATTATGGAGAGGGCAATGGGATTCTTGCCGCTTTGAATACCGAGCAGAAGGCAt                                     >  2:147535/1‑100 (MQ=255)
                                                           tgattgatGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCAt                                     >  2:1588109/1‑99 (MQ=255)
                                                           tgattgatGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATg                                    >  1:1394844/1‑100 (MQ=255)
                                                           tgattgatGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATg                                    >  2:1777859/1‑100 (MQ=255)
                                                           tgattgatGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATg                                    >  1:1227102/1‑100 (MQ=255)
                                                           tgattgatGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATg                                    <  2:459727/100‑1 (MQ=255)
                                                            gattgatGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATgg                                   >  1:1254313/1‑100 (MQ=255)
                                                            gattgatGAAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATgg                                   <  1:1348964/100‑1 (MQ=255)
                                                               tgatgaAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGgtgt                                >  2:979773/1‑100 (MQ=255)
                                                               tgatgaAGGTTATCAGGTGGTTCAATTTGGCTCAGCGAAAGATCATTAAGCGGTCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGTCATGgtg                                 >  2:167901/1‑99 (MQ=255)
                                                                gatgaAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTc                               >  2:705241/1‑100 (MQ=255)
                                                                 atgaAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCg                              >  1:957275/1‑100 (MQ=255)
                                                                  tgaAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCgg                             >  2:1049757/1‑100 (MQ=255)
                                                                   gaAGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGa                            >  1:92313/1‑100 (MQ=255)
                                                                    aaGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGaa                           >  2:781510/1‑100 (MQ=255)
                                                                    aaGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGaa                           <  1:1588120/100‑1 (MQ=255)
                                                                     aGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAAc                          <  1:20616/100‑1 (MQ=255)
                                                                     aGGTTATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAAc                          >  2:1686306/1‑100 (MQ=255)
                                                                        ttATCAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTg                       >  2:1919641/1‑100 (MQ=255)
                                                                            cAGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCgg                   <  1:822270/100‑1 (MQ=255)
                                                                             aGGTGGTTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCggg                  <  2:1382746/100‑1 (MQ=255)
                                                                              ggtggtTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCgggg                 <  2:1283506/100‑1 (MQ=255)
                                                                              ggtggtTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCgggg                 <  2:1817338/100‑1 (MQ=255)
                                                                               gtggtTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCggggg                <  2:1272067/100‑1 (MQ=255)
                                                                                tggtTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGa               <  2:49926/100‑1 (MQ=255)
                                                                                 ngtTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGaa              >  2:1631840/2‑100 (MQ=255)
                                                                                 ggtTCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGaa              >  1:986932/1‑100 (MQ=255)
                                                                                   ttCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGAAac            <  1:702765/100‑1 (MQ=255)
                                                                                   ttCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGAAac            <  1:564227/100‑1 (MQ=255)
                                                                                   ttCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGAAac            <  2:1675985/100‑1 (MQ=255)
                                                                                    tCTGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGTGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGAAaca           >  2:1258725/1‑100 (MQ=255)
                                                                                      tGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGAAAcaca         >  1:562583/1‑100 (MQ=255)
                                                                                      tGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGAAAcaca         >  1:534337/1‑100 (MQ=255)
                                                                                      tGTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGAAAcaca         >  1:1349145/1‑100 (MQ=255)
                                                                                       gTTTGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGAAACACAg        >  1:1412496/1‑100 (MQ=255)
                                                                                          tGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGAAACACAGCtt     >  1:1171393/1‑100 (MQ=255)
                                                                                          tGGCTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGAAACACAGCtt     <  1:290849/100‑1 (MQ=255)
                                                                                             cTCAGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGAAACACAGCTTGAt  >  1:1914425/1‑100 (MQ=255)
                                                                                                |                                                                                                
GGAGTTTGGTCCGGCAAAACGCTGGCCACACTACCACTATGCGGAGCTGGCAAAGCAGCTGATTGATGAAGGTTATCAGGTGGTTCTGTTTGGCTCGGCGAAAGATCATGAAGCGGGCAATGAGATTCTTGCCGCTTTGAATACCGAGCAGCAGGCATGGTGTCGGAACCTGGCGGGGGAAACACAGCTTGAT  >  NC_000913/3795495‑3795687

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