Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,794,970 G→C 6.9% M14I (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,9700GC 6.9% 718.3 / 31.9 201M14I (ATG→ATCwaaFADP‑heptose:LPS heptosyltransferase II
Reads supporting (aligned to +/- strand):  ref base G (75/111);  new base C (8/6);  total (83/118)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 226 total aligned reads displayed.

TCTGACAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTCACGGAGC  >  NC_000913/3794823‑3795116
                                                                                                                                                   |                                                                                                                                                  
tCTGACAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatga                                                                                                                                                   <  2:31118/149‑1 (MQ=255)
tCTGACAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatga                                                                                                                                                   <  3:90379/149‑1 (MQ=255)
tCTGACAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatga                                                                                                                                                   >  3:135910/1‑149 (MQ=255)
  tGACAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatgatg                                                                                                                                                 >  4:167733/1‑149 (MQ=255)
  tGACAAATATGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCCTAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGAcatcatc                                                                                                                                                 <  2:388564/149‑2 (MQ=255)
   gACAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatgatga                                                                                                                                                <  4:318961/149‑1 (MQ=255)
   gACAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatgatga                                                                                                                                                <  2:340924/149‑1 (MQ=255)
    aCAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGCGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatgatgat                                                                                                                                               >  1:83018/1‑149 (MQ=255)
    aCAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGAcatcatcat                                                                                                                                               <  1:153603/149‑1 (MQ=255)
    aCAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatgatgat                                                                                                                                               >  1:371676/1‑149 (MQ=255)
       aaTCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTc                                                                                                                                            <  1:377456/149‑1 (MQ=255)
         tCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATCATCATGTCGc                                                                                                                                          <  1:401141/149‑1 (MQ=255)
          cTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATCATCATGTCGcc                                                                                                                                         <  1:442581/149‑2 (MQ=255)
            gcgcgcGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCaaa                                                                                                                                       <  3:58306/149‑1 (MQ=255)
              gcgcgGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGt                                                                                                                                     >  2:83787/1‑149 (MQ=255)
                 cggcggGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGtctc                                                                                                                                  <  4:47854/149‑1 (MQ=255)
                     gggTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATc                                                                                                                              >  4:35157/1‑149 (MQ=255)
                     gggTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGACGTCGCAAAGTCTCTATc                                                                                                                              <  1:291820/149‑1 (MQ=255)
                           cGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGc                                                                                                                        <  2:1808/149‑1 (MQ=255)
                            gACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCt                                                                                                                       >  1:383477/1‑149 (MQ=255)
                               aaaCCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCa                                                                                                                    <  2:412513/149‑1 (MQ=255)
                                aaCCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAg                                                                                                                   <  4:71439/149‑1 (MQ=255)
                                    gTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatgat                                                                                                                                                  >  1:125090/1‑114 (MQ=255)
                                    gTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGgcgc                                                                                                               <  4:281459/149‑1 (MQ=255)
                                     ttCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGgcgcg                                                                                                              >  2:476218/1‑149 (MQ=255)
                                     ttCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGgcgcg                                                                                                              <  2:38142/149‑1 (MQ=255)
                                     ttCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGgcgcg                                                                                                              >  2:351555/1‑149 (MQ=255)
                                        aaaaCCGTGGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTa                                                                                                           >  3:437395/1‑149 (MQ=255)
                                         aaaCCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTAt                                                                                                          >  4:291204/1‑149 (MQ=255)
                                         aaaCCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTAt                                                                                                          >  2:275595/1‑149 (MQ=255)
                                               ttGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAg                                                                                                    <  4:184213/149‑1 (MQ=255)
                                               ttGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAg                                                                                                    <  2:154683/149‑1 (MQ=255)
                                                tGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGG                                                                                                   <  4:171956‑M1/149‑1 (MQ=255)
                                                 gCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGC                                                                                                  <  2:368862‑M1/149‑1 (MQ=255)
                                                   tGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGA                                                                                                >  4:2889‑M1/1‑149 (MQ=255)
                                                    gAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGAT                                                                                               >  3:260717‑M1/1‑149 (MQ=255)
                                                    gAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGAT                                                                                               <  2:430925‑M1/149‑1 (MQ=255)
                                                    gAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGAT                                                                                               <  3:429275‑M1/149‑1 (MQ=255)
                                                     aaGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGgcgc                                                                                                               >  2:454064/1‑132 (MQ=255)
                                                        gtgtAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATc                                                                                           <  4:429048/149‑1 (MQ=255)
                                                          gtAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGa                                                                                         <  1:415327/149‑1 (MQ=255)
                                                          gtAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGa                                                                                         <  2:149183/149‑1 (MQ=255)
                                                          gtAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGa                                                                                         >  2:53757/1‑149 (MQ=255)
                                                           tAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGAt                                                                                        <  4:48283/149‑1 (MQ=255)
                                                             aCGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGC                                                                                                  <  1:409977‑M1/137‑1 (MQ=255)
                                                             aCGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGC                                                                                                  >  2:409977‑M1/1‑137 (MQ=255)
                                                               ggAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGa                                                                                    <  3:171566/149‑1 (MQ=255)
                                                                      aTGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGAt                                                                                   >  3:413160/1‑143 (MQ=255)
                                                                      aTGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGAt                                                                                   <  4:413160/143‑1 (MQ=255)
                                                                      aTGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCAcc                                                                             >  4:191867/1‑149 (MQ=255)
                                                                              gCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATc                                                                                                                              <  3:289454/92‑1 (MQ=255)
                                                                              gCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGa                                                                                         >  1:141604/1‑129 (MQ=255)
                                                                              gCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGt                                                                     <  1:2789/149‑1 (MQ=255)
                                                                              gCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGt                                                                     <  4:251738/149‑1 (MQ=255)
                                                                                 gAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCAt                                                                        >  4:122741/1‑143 (MQ=255)
                                                                                    tCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATcgatgtgacggcaccatgccggtgccatc                                                               <  4:437395‑M1/149‑30 (MQ=255)
                                                                                       cgACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCat                                                            <  3:121066/149‑1 (MQ=255)
                                                                                        gACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATTATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATcgatgtgacggcaccatgccggtgccat                                                                <  3:448417‑M1/144‑29 (MQ=255)
                                                                                        gACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATcgatgtgacggcaccatgccggtgccat                                                                >  4:448417‑M1/1‑116 (MQ=255)
                                                                                        gACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGcc                                                                  >  2:289105/1‑142 (MQ=255)
                                                                                        gACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCatt                                                           >  4:174251/1‑149 (MQ=255)
                                                                                         aCGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTa                                                          <  2:75311‑M2/149‑2 (MQ=255)
                                                                                         aCGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTa                                                          >  2:259852‑M2/1‑148 (MQ=255)
                                                                                            cATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTatta                                                       >  3:147397‑M2/1‑145 (MQ=255)
                                                                                             aTAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTattat                                                      >  4:160487‑M2/1‑144 (MQ=255)
                                                                                             aTAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCCATGTGATGGCACCGGCATGGTGCCGTCCATTattat                                                      <  1:95153‑M2/149‑6 (MQ=255)
                                                                                              tAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTattatc                                                     <  3:463143‑M2/149‑7 (MQ=255)
                                                                                                 gagCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGa                                                                                    <  1:196125/115‑1 (MQ=255)
                                                                                                  agCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATgg                                                                                 <  1:10180/117‑1 (MQ=255)
                                                                                                  agCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGACGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATgg                                                                                 >  2:10180/1‑117 (MQ=255)
                                                                                                   gCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATaatggacggcaccatgccggtgccatcacatcgattatcgcctg                                                <  3:84540‑M1/149‑45 (MQ=255)
                                                                                                     tctGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatgatg                                                                                                                                                 >  1:309315/1‑50 (MQ=255)
                                                                                                     tctGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATg                                              >  1:53700/1‑149 (MQ=255)
                                                                                                         cATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGCCATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCgg                                          >  1:214272/1‑149 (MQ=255)
                                                                                                          aTGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGc                                             >  1:431130/1‑145 (MQ=255)
                                                                                                          aTGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGa                                         >  2:206310/1‑149 (MQ=255)
                                                                                                           tGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGAt                                                                                        <  3:223503/101‑1 (MQ=255)
                                                                                                           tGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGaa                                        <  3:194975/149‑1 (MQ=255)
                                                                                                              aaaTACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTAt                                                                                                                               <  2:296656/59‑1 (MQ=255)
                                                                                                              aaaTACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGtt                                     >  1:323729/1‑149 (MQ=255)
                                                                                                                aTACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTaa                                   <  3:333546/149‑1 (MQ=255)
                                                                                                                aTACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTaa                                   <  3:191867/149‑1 (MQ=255)
                                                                                                                    tGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGAt                                                                                   <  4:276133/97‑1 (MQ=255)
                                                                                                                    tGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGAt                                                                                   >  3:276133/1‑97 (MQ=255)
                                                                                                                    tGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGaa                               <  1:351555/149‑1 (MQ=255)
                                                                                                                      gTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGc                             >  1:224055/1‑149 (MQ=255)
                                                                                                                      gTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGc                             >  1:467388/1‑149 (MQ=255)
                                                                                                                         aTCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTCTCGCGGATGCCGGAAGTTAACGAAGCTAt                          <  3:387212/149‑1 (MQ=255)
                                                                                                                           cGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTc                        <  4:378200/149‑1 (MQ=255)
                                                                                                                           cGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTc                        <  1:5746/149‑1 (MQ=255)
                                                                                                                               ccGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTAt                    <  2:323729/149‑1 (MQ=255)
                                                                                                                                 gTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGc                  <  3:363564/149‑1 (MQ=255)
                                                                                                                                         tttgCGACATCATCATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTa                           <  3:284024/129‑1 (MQ=255)
                                                                                                                                         ttGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGt          >  2:358125/1‑149 (MQ=255)
                                                                                                                                              gACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTCACgg     >  4:197144/1‑149 (MQ=255)
                                                                                                                                               aCATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTCACGGa    <  3:211433/149‑1 (MQ=255)
                                                                                                                                                 atgatgatgTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTCACGGAgc  <  2:440564/149‑1 (MQ=255)
                                                                                                                                                   gatgatgTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTAt                    <  2:63353/129‑1 (MQ=255)
                                                                                                                                                   gatgatgTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATg                   >  2:26776/1‑130 (MQ=255)
                                                                                                                                                   gatgatgTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATg                   <  1:26776/130‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCTGACAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATGCCGGAAGTTAACGAAGCTATTCCTATGCCTCTCGGTCACGGAGC  >  NC_000913/3794823‑3795116

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