Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,794,973 G→C 5.7% 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 5.7% 1001.7 / 36.5 282M15I (ATG→ATCwaaFADP‑heptose:LPS heptosyltransferase II
Reads supporting (aligned to +/- strand):  ref base G (123/143);  new base C (7/9);  total (130/152)
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.58e-01

Only 100 of 305 total aligned reads displayed.

TGACAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTCACGGAGCGCTGG  >  NC_000913/3794825‑3795121
                                                                                                                                                    |                                                                                                                                                    
  aCAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatgatgat                                                                                                                                                    <  4:581322/149‑1 (MQ=255)
  aCAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatgatgat                                                                                                                                                    <  4:388240/149‑1 (MQ=255)
  aCAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatgatgat                                                                                                                                                    >  1:297818/1‑149 (MQ=255)
     aaTCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTc                                                                                                                                                 <  3:277914/149‑1 (MQ=255)
     aaTCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTc                                                                                                                                                 <  2:566184/149‑1 (MQ=255)
       tCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATCATCATGTCGc                                                                                                                                               <  2:347942/149‑1 (MQ=255)
       tCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATCATCATGTCGc                                                                                                                                               <  3:316294/149‑1 (MQ=255)
        cTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATCATCATGTCGcc                                                                                                                                              <  1:647477/149‑2 (MQ=255)
          gcgcgcGGCTGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCaaa                                                                                                                                            <  1:418435/149‑1 (MQ=255)
          gcgcgcGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCaaa                                                                                                                                            <  2:619829/149‑1 (MQ=255)
          gcgcgcGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCaaa                                                                                                                                            <  2:92419/149‑1 (MQ=255)
          gcgcgcGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCaaa                                                                                                                                            <  3:593939/149‑1 (MQ=255)
                  cggGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTAt                                                                                                                                    <  2:420439/149‑1 (MQ=255)
                   gggTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATc                                                                                                                                   <  3:437221/149‑1 (MQ=255)
                          gACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCt                                                                                                                            >  4:193935/1‑149 (MQ=255)
                          gACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCt                                                                                                                            >  2:175581/1‑149 (MQ=255)
                          gACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGAAATACATAGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGAAATGATGATGTCGCAAAGTCTCTATCGCACGCt                                                                                                                            >  3:638300/1‑149 (MQ=255)
                            cAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTcc                                                                                                                          <  3:53733/149‑1 (MQ=255)
                              aaCCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAg                                                                                                                        <  3:166413/149‑1 (MQ=255)
                              aaCCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAg                                                                                                                        <  3:553649/149‑1 (MQ=255)
                               aCCGATCAAAACCGTTGATGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAATAGCTCTGCTTGAAAATAAGGGTGATCGGCCCGTCTTGCGTTGGCGACATGATGATCTCGCAAAGTCTCTATCGCACGCTCCAgg                                                                                                                       >  2:647598/1‑149 (MQ=255)
                                  gTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGgcgc                                                                                                                    >  4:446402/1‑149 (MQ=255)
                                   ttCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCaaa                                                                                                                                            <  2:444882/124‑1 (MQ=255)
                                   ttCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGgcgcg                                                                                                                   >  2:365416/1‑149 (MQ=255)
                                        aaCCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATCATCATGTCGCcaac                                                                                                                                           <  1:360965/120‑2 (MQ=255)
                                           cGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATcccc                                                                                                           <  2:528596/149‑1 (MQ=255)
                                            gTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCa                                                                                                          >  4:702650/1‑149 (MQ=255)
                                             ttGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAg                                                                                                         <  1:260251/149‑1 (MQ=255)
                                                 tGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGA                                                                                                     >  3:371210‑M1/1‑149 (MQ=255)
                                                  gAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGAT                                                                                                    >  1:238324‑M1/1‑149 (MQ=255)
                                                       tgtAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCg                                                                                               >  2:312904/1‑149 (MQ=255)
                                                       tgtAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCg                                                                                               >  1:438571/1‑149 (MQ=255)
                                                       tgtAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAAT                                                                                                 <  2:299887‑M1/147‑1 (MQ=255)
                                                         tAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGAt                                                                                             <  1:555479/149‑1 (MQ=255)
                                                         tAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGAt                                                                                             <  4:558583/149‑1 (MQ=255)
                                                           aCGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGAtgt                                                                                           <  1:190374/149‑1 (MQ=255)
                                                           aCGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGACCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGAtgt                                                                                           >  4:215964/1‑149 (MQ=255)
                                                             ggAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGa                                                                                         >  3:620574/1‑149 (MQ=255)
                                                              gAATACGTGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGAt                                                                                        <  3:218274/149‑1 (MQ=255)
                                                              gAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGAt                                                                                        <  2:389615/149‑1 (MQ=255)
                                                              gAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGAt                                                                                        >  2:584449/1‑149 (MQ=255)
                                                               aaTACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATg                                                                                       <  3:379478/149‑1 (MQ=255)
                                                                 tACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGC                                                                                                       >  1:451640‑M1/1‑131 (MQ=255)
                                                                   cATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCAc                                                                                   >  3:734959/1‑149 (MQ=255)
                                                                    aTGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCAcc                                                                                  >  4:133762/1‑149 (MQ=255)
                                                                           ggCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGAT                                                                                                    >  1:535655‑M1/1‑124 (MQ=255)
                                                                            gCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGt                                                                          <  4:286354/149‑1 (MQ=255)
                                                                            gCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGt                                                                          <  1:151883/149‑1 (MQ=255)
                                                                              tGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGc                                                                        >  3:482983/1‑149 (MQ=255)
                                                                              tGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGc                                                                        >  4:64576/1‑149 (MQ=255)
                                                                               gAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGcc                                                                       >  4:107621/1‑149 (MQ=255)
                                                                               gAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGcc                                                                       <  2:238324/149‑1 (MQ=255)
                                                                                  tCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGAt                                                                                             >  4:749377/1‑124 (MQ=255)
                                                                                    gcgACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATcgatgtgacggcaccatgccggtgccatcac                                                                  <  1:482316‑M1/149‑32 (MQ=255)
                                                                                    gcgACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCa                                                                  <  4:497143/149‑1 (MQ=255)
                                                                                     cgACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCat                                                                 <  1:172348/149‑1 (MQ=255)
                                                                                      taCGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCAGCGGCATGGTGCCGTCTatt                                                                <  4:152479/148‑1 (MQ=255)
                                                                                      gACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCatt                                                                >  4:705562/1‑149 (MQ=255)
                                                                                         gCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATaatggacggcaccatgccggtgccatcacatcgatt                                                             <  3:195261‑M1/149‑37 (MQ=255)
                                                                                         gCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATaatggacggcaccatgccggtgccatcacatcgatt                                                             <  3:585212‑M1/149‑37 (MQ=255)
                                                                                         gCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTatt                                                             <  2:220692‑M2/149‑4 (MQ=255)
                                                                                          cATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTatta                                                            <  2:746806‑M2/149‑5 (MQ=255)
                                                                                            tAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTattatc                                                          <  3:469400‑M2/149‑7 (MQ=255)
                                                                                            tAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTattatc                                                          <  1:212365‑M2/149‑7 (MQ=255)
                                                                                            tAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTattatc                                                          <  2:306471‑M2/149‑7 (MQ=255)
                                                                                              agagCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCGAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCAACAATCGATGTGATTGCGCCGGCATGGTGCCGTCCATTattatcgc                                                        >  4:465753‑M2/1‑141 (MQ=255)
                                                                                              agagCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTattatcgc                                                        >  4:606719‑M2/1‑141 (MQ=255)
                                                                                              agagCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTattatcgc                                                        <  4:587113‑M2/149‑9 (MQ=255)
                                                                                              agagCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTattatcgc                                                        >  2:143230‑M2/1‑141 (MQ=255)
                                                                                                     tGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGcc                                                 <  1:30181/149‑1 (MQ=255)
                                                                                                        aTGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGa                                              >  1:761680/1‑149 (MQ=255)
                                                                                                         tGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGaa                                             <  4:263404/149‑1 (MQ=255)
                                                                                                          gAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAg                                            >  1:60625/1‑149 (MQ=255)
                                                                                                             aaTACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGa                                                     <  3:709816/137‑1 (MQ=255)
                                                                                                             aaTACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGa                                                     >  4:709816/1‑137 (MQ=255)
                                                                                                               tACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGTCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAAc                                       <  3:684755/149‑1 (MQ=255)
                                                                                                                  tGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGaa                                    <  3:1719/149‑1 (MQ=255)
                                                                                                                   ggTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAg                                   >  4:727453/1‑149 (MQ=255)
                                                                                                                     tGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCt                                 <  2:187504/149‑1 (MQ=255)
                                                                                                                      gATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTa                                >  4:678030/1‑149 (MQ=255)
                                                                                                                       aTCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTAt                               <  3:77604/149‑1 (MQ=255)
                                                                                                                       aTCGGCCAGTCTGGAGTTGAAGAAATTATGATATCGCAAAGTATCTATCGAACGCTCCAATCGCGCCATCCCCAGGCGATAATCGTTGTAATGGCACCGGAATGGTGCCGTCCATTATTATCCCTGATGACGGAATTTAACGCAGGTAt                               <  1:732134/149‑1 (MQ=255)
                                                                                                                         cGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTg                                                                         >  4:426181/1‑105 (MQ=255)
                                                                                                                          ggCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTcc                            >  4:142129/1‑149 (MQ=255)
                                                                                                                            cccGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTa                          >  2:504066/1‑149 (MQ=255)
                                                                                                                             ccGTCTTGGGTTGGCGACATGATGATGTCGGAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTAt                         <  1:393272/149‑1 (MQ=255)
                                                                                                                             ccGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTAt                         >  1:122358/1‑149 (MQ=255)
                                                                                                                             ccGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTAt                         <  2:660450/149‑1 (MQ=255)
                                                                                                                             ccGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATAATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTAt                         <  3:107621/149‑1 (MQ=255)
                                                                                                                                tCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGcc                      >  4:322098/1‑149 (MQ=255)
                                                                                                                                       tttgCGACATCATCATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTc                             <  3:252840/132‑1 (MQ=255)
                                                                                                                                       tttgCGACATCATCATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGt               >  3:424592/4‑149 (MQ=255)
                                                                                                                                       tttgCGACATCATCATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTACGCCTCTCGGt               >  2:550236/4‑149 (MQ=255)
                                                                                                                                       ttGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGCTAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGt               >  3:143765/1‑149 (MQ=255)
                                                                                                                                       ttGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGt               <  1:115564/149‑1 (MQ=255)
                                                                                                                                       atGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGTGATAATCGATGTGATGGCACCGGCATGGTGCCGACCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGt               >  2:595149/2‑149 (MQ=255)
                                                                                                                                          gCGACATCATCATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTCAc            >  1:518345/1‑149 (MQ=255)
                                                                                                                                             aCATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTCACGGa         <  1:742681/149‑1 (MQ=255)
                                                                                                                                               atgatgatgTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTCACGGAgc       >  1:314064/1‑149 (MQ=255)
                                                                                                                                                 gatgatgTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGaa                                             >  3:577274/1‑109 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGACAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTCACGGAGCGCTGG  >  NC_000913/3794825‑3795121

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