Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,794,973 G→C 6.1% M15I (ATG→ATC waaF → ADP‑heptose‑‑LPS heptosyltransferase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,794,9730GC 6.1% 892.2 / 33.2 260M15I (ATG→ATCwaaFADP‑heptose‑‑LPS heptosyltransferase 2
Reads supporting (aligned to +/- strand):  ref base G (122/121);  new base C (8/8);  total (130/130)
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.54e-01

Only 100 of 276 total aligned reads displayed.

CTGACAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTCACGGAGCGCTGG  >  NC_000913/3794824‑3795121
                                                                                                                                                     |                                                                                                                                                    
cTGACAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatgatg                                                                                                                                                      <  2:1103159/150‑1 (MQ=255)
cTGACAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatgatg                                                                                                                                                      <  1:1234580/150‑1 (MQ=255)
      gaTCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCg                                                                                                                                                <  2:446979/149‑1 (MQ=255)
      aaTCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCg                                                                                                                                                >  1:376155/1‑150 (MQ=255)
      aaTCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCg                                                                                                                                                <  2:1017833/150‑1 (MQ=255)
      aaTCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCg                                                                                                                                                <  1:1423981/150‑1 (MQ=255)
        tCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCGCGTCTTGGGTTGGCGACATTATGATGTCGCa                                                                                                                                              <  2:396877/150‑1 (MQ=255)
        tCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATCATCATGTCGcc                                                                                                                                              <  2:1686405/150‑2 (MQ=255)
         cTGCGCGCGGCGGGTTAAGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACTGAATACATGGCCTGGCTTAATCGCGACGCATAAGAGCTCTGCATGAAAATACTTGTGATCGGCCCGTCTTGGGTTGGCGAAATCATCATGTCGcca                                                                                                                                             <  2:1065390/150‑3 (MQ=255)
           gcgcgcGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATCATCATGTCGCcaac                                                                                                                                           <  2:1536013/150‑2 (MQ=255)
           gcgcgcGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATCATCATGTCGCcaac                                                                                                                                           <  2:1630637/150‑2 (MQ=255)
               gcggcgGGTTACGACCAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGtctc                                                                                                                                       <  1:440039/150‑1 (MQ=255)
                 ggcggGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTa                                                                                                                                     >  1:1581840/1‑150 (MQ=255)
                   cggGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATc                                                                                                                                   <  2:1021376/150‑1 (MQ=255)
                    gggTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCg                                                                                                                                  <  1:1358635/150‑1 (MQ=255)
                     ggTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGc                                                                                                                                 <  2:320548/150‑1 (MQ=255)
                      gTTACGACAAACCGGTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGGTTGCGACATGATGATGTCGCAAAGTCTCTATCGCa                                                                                                                                <  2:182866/150‑1 (MQ=255)
                       ttACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCAc                                                                                                                               <  2:50486/150‑1 (MQ=255)
                         aCGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGc                                                                                                                             <  1:1638908/150‑1 (MQ=255)
                         aCGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGc                                                                                                                             >  1:170116/1‑150 (MQ=255)
                          cGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCt                                                                                                                            >  1:325343/1‑150 (MQ=255)
                            aCAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTcc                                                                                                                          <  1:1935435/150‑1 (MQ=255)
                             cAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCa                                                                                                                         >  1:284969/1‑150 (MQ=255)
                             cAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCa                                                                                                                         <  1:649680/150‑1 (MQ=255)
                              aaaCCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAg                                                                                                                        >  2:1769268/1‑150 (MQ=255)
                               aaCCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAgg                                                                                                                       <  1:1103970/150‑1 (MQ=255)
                                    ttCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGgcgcgc                                                                                                                  <  1:975255/150‑1 (MQ=255)
                                       aaaaCCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTAt                                                                                                               >  1:1492636/1‑150 (MQ=255)
                                             gTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAg                                                                                                         <  1:108927/150‑1 (MQ=255)
                                               tGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGC                                                                                                       <  2:1718212‑M1/150‑1 (MQ=255)
                                               tGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGC                                                                                                       <  1:1663770‑M1/150‑1 (MQ=255)
                                                 cTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGA                                                                                                     >  2:1164358‑M1/1‑150 (MQ=255)
                                                 cTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGA                                                                                                     >  1:9628‑M1/1‑150 (MQ=255)
                                                  tGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGAT                                                                                                    >  2:91246‑M1/1‑150 (MQ=255)
                                                      ggTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATc                                                                                                <  1:1552674/150‑1 (MQ=255)
                                                      ggTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATc                                                                                                <  1:1034413/150‑1 (MQ=255)
                                                          tAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGAtg                                                                                            <  2:728253/150‑1 (MQ=255)
                                                           aaCGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCAGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGgcgcg                                                                                                                   >  2:820877/1‑126 (MQ=255)
                                                           aaCGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCAGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGgcgcg                                                                                                                   <  1:820877/126‑1 (MQ=255)
                                                                    cATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTa                                                                                                                >  2:1435903/1‑120 (MQ=255)
                                                                    cATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCAcc                                                                                  >  2:977611/1‑150 (MQ=255)
                                                                       ggCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGc                                                                               >  2:514457/1‑150 (MQ=255)
                                                                             gCTGAATCGCGACGCATAAGAGCTCTGCATTAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTg                                                                         <  2:1297309/150‑1 (MQ=255)
                                                                              cTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGc                                                                        <  2:1134032/150‑1 (MQ=255)
                                                                              cTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGAGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGc                                                                        >  1:860465/1‑150 (MQ=255)
                                                                               tGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGcc                                                                       >  2:1694694/1‑150 (MQ=255)
                                                                                gAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCg                                                                      <  1:765732/150‑1 (MQ=255)
                                                                                 aaTCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGt                                                                     >  1:669318/1‑150 (MQ=255)
                                                                                      cgACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCatt                                                                >  1:850840/1‑150 (MQ=255)
                                                                                      cgACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCatt                                                                <  1:1697775/150‑1 (MQ=255)
                                                                                          gCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTatta                                                            >  2:1183426‑M2/1‑146 (MQ=255)
                                                                                               agagCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATcgcg                                                       <  2:669318/150‑1 (MQ=255)
                                                                                               agagCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATcgcg                                                       >  2:1351845/1‑150 (MQ=255)
                                                                                               agagCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATcgcg                                                       <  1:465743/150‑1 (MQ=255)
                                                                                                gagCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCgg                                                      <  1:349701/150‑1 (MQ=255)
                                                                                                 agCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGa                                                     <  1:91246/150‑1 (MQ=255)
                                                                                                      tGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATcgatgtgacggcaccatgccggtgccatcacatcgattatcgcctgggg                                                <  1:234809‑M1/150‑50 (MQ=255)
                                                                                                       gCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATg                                                   >  1:1755725/1‑146 (MQ=255)
                                                                                                       gCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCgg                                               >  2:531597/1‑150 (MQ=255)
                                                                                                        cATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGa                                              >  2:957015/1‑150 (MQ=255)
                                                                                                          tGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAg                                            >  1:1006224/1‑150 (MQ=255)
                                                                                                          tGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAg                                            >  1:18674/1‑150 (MQ=255)
                                                                                                          tGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAg                                            >  1:21435/1‑150 (MQ=255)
                                                                                                           gAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGt                                           <  1:1769268/150‑1 (MQ=255)
                                                                                                             aaaTACTGGTGATCGGCCCGTCTTGGGTGGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTa                                         <  2:1349378/150‑1 (MQ=255)
                                                                                                              aaTACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTaa                                        <  2:376155/150‑1 (MQ=255)
                                                                                                               aTACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAAc                                       <  2:192272/150‑1 (MQ=255)
                                                                                                                tACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCTCAAAGTCTCTATCGCACGCTCCAGTCTTGATATCACCAGGCGATAATCTATGTTATGGCACCGGAATGGTTTTGTCTATTATTATCGCTGATGCTTTAAGTTaagg                                      <  2:1343962/150‑3 (MQ=255)
                                                                                                                 aCTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCgg                                                      >  2:1867250/1‑133 (MQ=255)
                                                                                                                 aCTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGa                                     >  1:101993/1‑150 (MQ=255)
                                                                                                                 aCTGGTGATCGGACCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGa                                     <  1:1003602/150‑1 (MQ=255)
                                                                                                                    ggTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGc                                  <  1:723225/150‑1 (MQ=255)
                                                                                                                    ggTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGc                                  >  1:594973/1‑150 (MQ=255)
                                                                                                                    ggTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGc                                  <  2:1482100/150‑1 (MQ=255)
                                                                                                                    ggTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGc                                  <  2:1581840/150‑1 (MQ=255)
                                                                                                                      tGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTa                                >  1:1327358/1‑150 (MQ=255)
                                                                                                                      tGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTa                                >  2:288909/1‑150 (MQ=255)
                                                                                                                        aTCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTAtt                              >  1:1312625/1‑150 (MQ=255)
                                                                                                                          cGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTcc                            >  1:1107165/1‑150 (MQ=255)
                                                                                                                           ggCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCat                                                                 >  2:239617/1‑112 (MQ=255)
                                                                                                                           ggCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCt                           >  2:348741/1‑150 (MQ=255)
                                                                                                                            gCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGa                                     <  1:1475070/139‑1 (MQ=255)
                                                                                                                            gCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTa                          <  1:1213779/150‑1 (MQ=255)
                                                                                                                              ccGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATg                        <  1:1842287/150‑1 (MQ=255)
                                                                                                                              ccGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATg                        <  2:606477/150‑1 (MQ=255)
                                                                                                                                    tGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCtctc                  >  2:583552/1‑150 (MQ=255)
                                                                                                                                      ggTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCgg                <  1:777138/150‑1 (MQ=255)
                                                                                                                                       ctttgCGACATCATCATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGt               >  2:217911/5‑150 (MQ=255)
                                                                                                                                        tttgCGACATCATCATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTc              >  1:1886358/4‑150 (MQ=255)
                                                                                                                                        ttGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTc              <  2:170116/150‑1 (MQ=255)
                                                                                                                                          ggCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTCAc            <  1:1164358/150‑1 (MQ=255)
                                                                                                                                            cGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTCACgg          >  2:1911651/1‑150 (MQ=255)
                                                                                                                                            cGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTCACgg          >  1:1922486/1‑150 (MQ=255)
                                                                                                                                            cGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTCACgg          <  2:249493/150‑1 (MQ=255)
                                                                                                                                              aCATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTCACGGAg        <  2:1509710/150‑1 (MQ=255)
                                                                                                                                                 tgatgatgTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTCACGGAgcgc     >  2:1936550/1‑150 (MQ=255)
                                                                                                                                                 tgatgatgTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTCACGGAgcgc     >  2:875489/1‑150 (MQ=255)
                                                                                                                                                  gatgatgTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTCACGGAGCGCt    >  2:1629413/1‑150 (MQ=255)
                                                                                                                                                   atgatgTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTCACGGAGCGCTg   >  2:74104/1‑150 (MQ=255)
                                                                                                                                                    tgatgTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTCACGGAGCGCTgg  >  1:410486/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTGACAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTCACGGAGCGCTGG  >  NC_000913/3794824‑3795121

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