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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,794,9730GC 6.9% 721.1 / 32.3 201M15I (ATG→ATCwaaFADP‑heptose:LPS heptosyltransferase II
Reads supporting (aligned to +/- strand):  ref base G (75/112);  new base C (8/6);  total (83/118)
Fisher's exact test for biased strand distribution p-value = 2.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

Only 100 of 220 total aligned reads displayed.

TGACAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTCACGGAGCGCTGG  >  NC_000913/3794825‑3795121
                                                                                                                                                    |                                                                                                                                                    
tGACAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatgatg                                                                                                                                                      >  1:152507/1‑149 (MQ=255)
tGACAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatgatg                                                                                                                                                      >  4:167733/1‑149 (MQ=255)
 gACAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatgatga                                                                                                                                                     >  1:154683/1‑149 (MQ=255)
 gACAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatgatga                                                                                                                                                     >  1:430925/1‑149 (MQ=255)
  aCAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGCGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatgatgat                                                                                                                                                    >  1:83018/1‑149 (MQ=255)
      aTCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCg                                                                                                                                                >  2:261447/1‑149 (MQ=255)
          gcgcgcGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGGCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCaaa                                                                                                                                            <  3:109179/149‑1 (MQ=255)
            gcgcgGCGGGTTACGACAAACCGTTCAACCCCGTTGCTGAAGGTGTAACGGAATACGTTGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGCTCGGCGACATGATGATGTCGCAAAGt                                                                                                                                          <  3:415568/149‑1 (MQ=255)
            gcgcgGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGt                                                                                                                                          >  2:83787/1‑149 (MQ=255)
               cggcggGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGtctc                                                                                                                                       <  4:47854/149‑1 (MQ=255)
                ggcggGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTctct                                                                                                                                      <  2:151776/149‑1 (MQ=255)
                          gACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATCATCATGTCGcc                                                                                                                                              <  1:440050/131‑2 (MQ=255)
                          gACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATCATCATGTCGCcaa                                                                                                                                            <  1:383202/133‑1 (MQ=255)
                             aaaCCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCa                                                                                                                         <  2:412513/149‑1 (MQ=255)
                              aaCCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAg                                                                                                                        <  4:71439/149‑1 (MQ=255)
                                   ttCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGgcgcg                                                                                                                   >  1:358297/1‑149 (MQ=255)
                                   ttCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGgcgcg                                                                                                                   >  2:351555/1‑149 (MQ=255)
                                   ttCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGgcgcg                                                                                                                   <  2:38142/149‑1 (MQ=255)
                                   ttCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGgcgcg                                                                                                                   >  2:476218/1‑149 (MQ=255)
                                      aaaaCCGTTGCTGAACGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGGTGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCCCGCTa                                                                                                                >  2:290390/1‑149 (MQ=255)
                                      aaaaCCGTGGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTa                                                                                                                >  3:437395/1‑149 (MQ=255)
                                       aaaCCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTAt                                                                                                               >  4:291204/1‑149 (MQ=255)
                                        aaCCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATc                                                                                                              <  4:280272/149‑1 (MQ=255)
                                             ttGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAg                                                                                                         <  4:184213/149‑1 (MQ=255)
                                              tGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGG                                                                                                        <  4:171956‑M1/149‑1 (MQ=255)
                                                 tGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGA                                                                                                     >  3:49768‑M1/1‑149 (MQ=255)
                                                  gAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGAT                                                                                                    >  3:260717‑M1/1‑149 (MQ=255)
                                                   aaGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGgcgc                                                                                                                    <  1:454064/132‑1 (MQ=255)
                                                   aaGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGgcgc                                                                                                                    >  2:454064/1‑132 (MQ=255)
                                                    aGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAA                                                                                                  <  1:476218‑M1/149‑1 (MQ=255)
                                                      gtgtAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATc                                                                                                <  4:429048/149‑1 (MQ=255)
                                                       tgtAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCg                                                                                               >  3:172354/1‑149 (MQ=255)
                                                         tAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGAt                                                                                             >  4:194975/1‑149 (MQ=255)
                                                         tAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGAt                                                                                             <  4:48283/149‑1 (MQ=255)
                                                           aCGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGAtgt                                                                                           <  1:127921/149‑1 (MQ=255)
                                                           aCGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGAtgt                                                                                           >  2:2789/1‑149 (MQ=255)
                                                           aCGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGC                                                                                                       >  2:409977‑M1/1‑137 (MQ=255)
                                                            cgcaaTACATGGCCTGGCGGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGACTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGAtgtg                                                                                          >  4:36572/4‑149 (MQ=255)
                                                             ggAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGa                                                                                         <  1:53757/149‑1 (MQ=255)
                                                                    aTGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGAt                                                                                        <  4:413160/143‑1 (MQ=255)
                                                                    aTGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGAt                                                                                        >  3:413160/1‑143 (MQ=255)
                                                                    aTGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCAcc                                                                                  >  4:191867/1‑149 (MQ=255)
                                                                         cTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCAt                                                                             >  4:333546/1‑149 (MQ=255)
                                                                            gCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATc                                                                                                                                   >  4:289454/1‑92 (MQ=255)
                                                                            gCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATc                                                                                                                                   <  3:289454/92‑1 (MQ=255)
                                                                            gCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGa                                                                                              <  2:141604/129‑1 (MQ=255)
                                                                               gAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCAt                                                                             >  4:122741/1‑143 (MQ=255)
                                                                               gAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCAt                                                                             <  3:122741/143‑1 (MQ=255)
                                                                               gAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGcc                                                                       <  1:394217/149‑1 (MQ=255)
                                                                                      gACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATTATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATcgatgtgacggcaccatgccggtgccat                                                                     <  3:448417‑M1/144‑29 (MQ=255)
                                                                                      gACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGcc                                                                       <  1:289105/142‑1 (MQ=255)
                                                                                      gACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCatt                                                                >  4:174251/1‑149 (MQ=255)
                                                                                       aCGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTa                                                               <  2:75311‑M2/149‑2 (MQ=255)
                                                                                         gCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATaatggacggcaccatgccggtgccatcacatcgatt                                                             <  1:376763‑M1/149‑37 (MQ=255)
                                                                                           aTAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTattat                                                           >  4:160487‑M2/1‑144 (MQ=255)
                                                                                            tAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATaatggacggcaccatgccggtgccatcacatcgattatc                                                          <  2:100574‑M1/149‑40 (MQ=255)
                                                                                            tAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTattatc                                                          >  2:370909‑M2/1‑143 (MQ=255)
                                                                                            tAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTattatc                                                          <  3:463143‑M2/149‑7 (MQ=255)
                                                                                              agagCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTattatcgc                                                        >  3:247796‑M2/1‑141 (MQ=255)
                                                                                               gagCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGa                                                                                         >  2:196125/1‑115 (MQ=255)
                                                                                               gagCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGa                                                                                         <  1:196125/115‑1 (MQ=255)
                                                                                                agCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATgg                                                                                      <  1:10180/117‑1 (MQ=255)
                                                                                                 gCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGa                                                     >  2:242648/1‑149 (MQ=255)
                                                                                                   tctGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatgatg                                                                                                                                                      <  2:309315/50‑1 (MQ=255)
                                                                                                   tctGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTAt                                                                                                                                    <  4:388063/68‑1 (MQ=255)
                                                                                                    ctGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGc                                                  <  4:463143/149‑1 (MQ=255)
                                                                                                       cATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCgg                                                                                >  1:414155/1‑116 (MQ=255)
                                                                                                        aTGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTg                                                                         >  1:348704/1‑122 (MQ=255)
                                                                                                        aTGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGc                                                  <  2:431130/145‑1 (MQ=255)
                                                                                                         tGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGAt                                                                                             >  4:223503/1‑101 (MQ=255)
                                                                                                         tGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGaa                                             <  3:194975/149‑1 (MQ=255)
                                                                                                         tGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGaa                                             <  3:96333/149‑1 (MQ=255)
                                                                                                            aaaTACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTAt                                                                                                                                    >  1:296656/1‑59 (MQ=255)
                                                                                                            aaaTACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGtt                                          >  1:323729/1‑149 (MQ=255)
                                                                                                              aTACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTaa                                        <  4:260717/149‑1 (MQ=255)
                                                                                                              aTACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTaa                                        <  3:333546/149‑1 (MQ=255)
                                                                                                              aTACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTaa                                        <  3:191867/149‑1 (MQ=255)
                                                                                                               tACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAAc                                       <  1:275595/149‑1 (MQ=255)
                                                                                                                  tGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGAt                                                                                        >  3:276133/1‑97 (MQ=255)
                                                                                                                    gTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGc                                  >  1:224055/1‑149 (MQ=255)
                                                                                                                       aTCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTCTCGCGGATGCCGGAAGTTAACGAAGCTAt                               <  3:387212/149‑1 (MQ=255)
                                                                                                                       aTCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTAt                               <  2:272841/149‑1 (MQ=255)
                                                                                                                          ggCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTcc                            >  3:73267/1‑149 (MQ=255)
                                                                                                                           gCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGGTGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCt                           <  2:383477/149‑1 (MQ=255)
                                                                                                                           gCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCt                           >  4:44312/1‑149 (MQ=255)
                                                                                                                             ccGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTAt                         <  2:53700/149‑1 (MQ=255)
                                                                                                                             ccGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTAt                         <  2:323729/149‑1 (MQ=255)
                                                                                                                             ccGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTAt                         <  3:163606/149‑1 (MQ=255)
                                                                                                                               gTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGc                       <  3:363564/149‑1 (MQ=255)
                                                                                                                                       tttgCGACATCATCATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTaa                                        <  1:404030/121‑1 (MQ=255)
                                                                                                                                       tttgCGACATCATCATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTa                                >  4:284024/4‑132 (MQ=255)
                                                                                                                                       tttgCGACATCATCATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTa                                <  3:284024/129‑1 (MQ=255)
                                                                                                                                       ttGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGt               >  2:358125/1‑149 (MQ=255)
                                                                                                                                             aCATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTCACGGa         <  3:211433/149‑1 (MQ=255)
                                                                                                                                             aCATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTCACGGa         <  2:341505/149‑1 (MQ=255)
                                                                                                                                               atgatgatgTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTCACGGAgc       >  2:169468/1‑149 (MQ=255)
                                                                                                                                               atgatgatgTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTCACGGAgc       <  2:440564/149‑1 (MQ=255)
                                                                                                                                                 gatgatgTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTAt                         <  2:63353/129‑1 (MQ=255)
                                                                                                                                                 gatgatgTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTAt                         >  1:63353/1‑129 (MQ=255)
                                                                                                                                                 gatgatgTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATg                        <  1:26776/130‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGACAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTCACGGAGCGCTGG  >  NC_000913/3794825‑3795121

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