Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,113,364 T→A 5.4% intergenic (‑303/+70) rfbB ← / ← wcaN dTDP‑glucose 4,6 dehydratase, NAD(P)‑binding/putative regulatory subunit for GalU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,113,3640TA 5.4% 426.6 / 7.3 139intergenic (‑303/+70)rfbB/wcaNdTDP‑glucose 4,6 dehydratase, NAD(P)‑binding/putative regulatory subunit for GalU
Reads supporting (aligned to +/- strand):  ref base T (64/67);  new base A (3/5);  total (67/72)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 140 total aligned reads displayed.

AAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGCCATACTTCACAAACG  >  NC_000913/2113217‑2113513
                                                                                                                                                    |                                                                                                                                                     
    ccGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACttt                                                                                                                                                 >  1:1532204/1‑150 (MQ=255)
     cGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACtttt                                                                                                                                                >  1:764794/1‑150 (MQ=255)
       aaaaaaTCAAACGCTAATTGTCTTAACAATCCGCGCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAATACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTtgc                                                                                                                                              <  2:339584/150‑1 (MQ=255)
       aaaaaTTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTtgc                                                                                                                                              <  2:2097047/150‑1 (MQ=255)
               aaaCGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGTAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAAtt                                                                                                                                      <  1:1412500/150‑1 (MQ=255)
                  cGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTc                                                                                                                                   <  1:1532156/150‑1 (MQ=255)
                   gCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAGCAAAAGTAATTTGTTGCGAATCTTcc                                                                                                                                  <  2:1725055/150‑2 (MQ=38)
                    cTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCAc                                                                                                                                 <  2:2042852/150‑1 (MQ=255)
                      aaTTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTg                                                                                                                               >  1:1888974/1‑150 (MQ=255)
                      aaTTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTCATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTAGCTGCGAATTTTCACTg                                                                                                                               >  2:299270/1‑150 (MQ=255)
                       aTTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGTGAAAATTCGCAGCAAAAGTAATTTGTTGCGAATCTTcctgc                                                                                                                               <  1:2025949/150‑5 (MQ=25)
                       aTTGTCTTAACAATCCGCTCTGGAAACAAGGATAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGTGAAACTTCGCAGCAAAAGTAATTTGTTGCGAATCTTcctgc                                                                                                                               <  1:317746/150‑5 (MQ=25)
                        ttGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGTGAAAATTCGCAGCAAAAGTAATTTGTTGCGAATCTTC‑CTGcc                                                                                                                             <  2:2208256/150‑1 (MQ=25)
                        ttGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGcc                                                                                                                             <  1:1913735/150‑1 (MQ=255)
                         tGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGTGAAAATTCGCAGCAAAAGTAATTTGTTGCGAATCTTC‑CTGCCg                                                                                                                            <  2:878492/150‑1 (MQ=25)
                         tGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCg                                                                                                                            >  1:69719/1‑150 (MQ=255)
                         tGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCg                                                                                                                            >  2:480438/1‑150 (MQ=255)
                         tGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCg                                                                                                                            >  2:1033497/1‑150 (MQ=255)
                            cTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGttt                                                                                                                         >  2:382572/1‑149 (MQ=255)
                                ccAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAAt                                                                                                                     >  1:1670446/1‑150 (MQ=255)
                                ccAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAAt                                                                                                                     >  1:456957/1‑150 (MQ=255)
                                  aaTCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAAttt                                                                                                                   >  1:2113166/1‑150 (MQ=255)
                                      cGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCtt                                                                                                               <  2:719972/150‑1 (MQ=255)
                                       gCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTa                                                                                                              <  1:805525/150‑1 (MQ=255)
                                       gCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTa                                                                                                              <  2:1003629/150‑1 (MQ=255)
                                         tctGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATc                                                                                                            >  2:111247/1‑150 (MQ=255)
                                              aagcAAGGAAAATCCTGGAAAACATTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACTGCAG‑GAAGATTCTCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCg                                                                                                       <  1:109450/147‑1 (MQ=255)
                                              aaaCAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCg                                                                                                       <  1:488670/150‑1 (MQ=255)
                                              aaaCAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCg                                                                                                       <  2:623430/150‑1 (MQ=255)
                                              aaaCAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCg                                                                                                       <  1:1857242/150‑1 (MQ=255)
                                              aaaCAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACAGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCg                                                                                                       <  2:223191/150‑1 (MQ=255)
                                                aCAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGtt                                                                                                     <  2:497386/150‑1 (MQ=255)
                                                 cAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTa                                                                                                    <  1:906498/150‑1 (MQ=255)
                                                 cAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTa                                                                                                    >  1:1056720/1‑150 (MQ=255)
                                                 aaaGGAAAATCCTTGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTa                                                                                                    <  2:491441/149‑1 (MQ=255)
                                                  aaGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTAc                                                                                                   >  1:561840/1‑150 (MQ=255)
                                                    ggAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACAt                                                                                                 >  1:447781/1‑150 (MQ=255)
                                                    ggAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACAt                                                                                                 >  1:1719949/1‑150 (MQ=255)
                                                    ggAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACAt                                                                                                 >  2:1413268/1‑150 (MQ=255)
                                                        aaTCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCgg                                                                                             <  1:517368/150‑1 (MQ=255)
                                                         aTCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGt                                                                                            <  2:878206/150‑1 (MQ=255)
                                                              ggAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCTGAt                                                                                       >  2:532289/1‑150 (MQ=255)
                                                              ggAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGAt                                                                                       >  1:2063449/1‑150 (MQ=255)
                                                                aaaaCTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGAttt                                                                                     >  1:1779475/1‑150 (MQ=255)
                                                                aaaaCTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATAAACCGTTACATCCGGTCAGAttt                                                                                     >  2:204083/1‑150 (MQ=255)
                                                                  aaCTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTc                                                                                   <  2:1601714/150‑1 (MQ=255)
                                                                   aCTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCa                                                                                  <  2:1803621/150‑1 (MQ=255)
                                                                   aCTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCa                                                                                  >  1:1466060/1‑150 (MQ=255)
                                                                    cTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCat                                                                                 >  2:955164/1‑150 (MQ=255)
                                                                    atttGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCat                                                                                 <  1:2028069/149‑1 (MQ=255)
                                                                     tttGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCatt                                                                                >  2:1444252/1‑150 (MQ=255)
                                                                      ttGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCatta                                                                               >  2:1082056/1‑150 (MQ=255)
                                                                            aaacCCCCACTCCTAACTCCTTGTTATTCTGATGGTCTATATAAAACAACGGCAG‑GAAGAATCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCt                                                                         <  2:1466060/150‑1 (MQ=255)
                                                                            aaaaCCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCt                                                                         <  2:398907/150‑1 (MQ=255)
                                                                            aaaaCCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCt                                                                         <  2:1425017/150‑1 (MQ=255)
                                                                            aaaaCCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCt                                                                         <  2:301575/150‑1 (MQ=255)
                                                                              aaCCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTa                                                                       <  1:1826277/150‑1 (MQ=255)
                                                                                 ccTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACa                                                                    <  2:1888974/150‑1 (MQ=255)
                                                                                 ccTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACa                                                                    <  2:1577351/150‑1 (MQ=255)
                                                                                     cTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCTGTCAGATTTTCATTATTCGCTTAACAGCtt                                                                >  1:100810/1‑150 (MQ=255)
                                                                                       gCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTtct                                                              >  1:848167/1‑150 (MQ=255)
                                                                                       gCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTtct                                                              >  1:1465092/1‑150 (MQ=255)
                                                                                       gCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATAAACCGTTACATACGGTCAGATTTTCATTATTCGCTTAACAGCTtct                                                              >  1:1988222/1‑150 (MQ=255)
                                                                                         tAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCa                                                            <  1:1533472/150‑1 (MQ=255)
                                                                                         tAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCa                                                            <  1:382572/150‑1 (MQ=255)
                                                                                         tAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCa                                                            <  1:1239167/150‑1 (MQ=255)
                                                                                         tAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCa                                                            <  1:1033497/150‑1 (MQ=255)
                                                                                           aCTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAAt                                                          >  2:616994/1‑150 (MQ=255)
                                                                                            cTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATa                                                         >  1:652357/1‑150 (MQ=255)
                                                                                             tCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATAc                                                        <  2:2121996/150‑1 (MQ=255)
                                                                                             tCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATAc                                                        <  1:1678719/150‑1 (MQ=255)
                                                                                              cGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATAcc                                                       <  1:621692/150‑1 (MQ=255)
                                                                                                ttgttATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTTCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCTGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCtt                                                     <  2:2097417/150‑1 (MQ=255)
                                                                                                  gttATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTa                                                   >  1:1012077/1‑150 (MQ=255)
                                                                                                     aTTCTGATGGTTTATATATAACAACTGCAG‑GAAGATTAGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTTTAATTTTCTTATCAACCGTTACATCCGGTCAGTTTTTGATTATTCGGTTAACAGGTTCTCAATACCTTTACgg                                                <  2:154518/150‑1 (MQ=255)
                                                                                                     aTTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACgg                                                <  1:1719576/150‑1 (MQ=255)
                                                                                                     aTTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACgg                                                <  1:2117852/150‑1 (MQ=255)
                                                                                                        cTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAAc                                             >  2:1862002/1‑150 (MQ=255)
                                                                                                         tGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTAc                                                  <  1:654481/144‑1 (MQ=255)
                                                                                                         tGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTAc                                                  >  2:654481/1‑144 (MQ=255)
                                                                                                           aTGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTc                                          <  1:1249617/150‑1 (MQ=255)
                                                                                                             ggTTTATATAAAACAACGGCAT‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCAGTGCCGTTATAATTTTCTTATCAACCGTAACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGc                                        <  2:1804861/150‑1 (MQ=255)
                                                                                                                   tataAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCcttg                                  >  1:1707418/1‑149 (MQ=255)
                                                                                                                   tataAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCcttc                                  >  1:695871/1‑150 (MQ=255)
                                                                                                                    ataAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCcttct                                 >  2:1967223/1‑150 (MQ=255)
                                                                                                                    ataAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCcttct                                 <  1:1790502/150‑1 (MQ=255)
                                                                                                                     taAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGATTTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGACCcttctt                                >  1:211660/1‑150 (MQ=255)
                                                                                                                     taAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAAGAGCTTCTCAATACCTTTACGGAACTTCGCCCcttctt                                >  2:624426/1‑150 (MQ=255)
                                                                                                                     taAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCcttctt                                >  2:680259/1‑150 (MQ=255)
                                                                                                                      aaaaCAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCttt                               <  1:2071060/150‑1 (MQ=255)
                                                                                                                       aaaCAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTc                              >  2:1756516/1‑150 (MQ=255)
                                                                                                                       aaaCAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTc                              >  1:629909/1‑150 (MQ=255)
                                                                                                                        tataacGGCAGTGAAAATTCGCAGCAAAAGTAATTTGTTGCGAATCTTC‑CTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCa                             >  2:1230330/4‑150 (MQ=21)
                                                                                                                                cAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCTTTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGt                     <  2:1599319/150‑1 (MQ=255)
                                                                                                                                     aaGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGc                 >  2:972621/1‑150 (MQ=255)
                                                                                                                                      aGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGcc                <  2:1206678/150‑1 (MQ=255)
                                                                                                                                         ttCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGCCATa             >  2:1737027/1‑150 (MQ=255)
                                                                                                                                            gCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGCCATACtt          >  2:795309/1‑150 (MQ=255)
                                                                                                                                            gCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGCAATACtt          >  1:204074/1‑150 (MQ=255)
                                                                                                                                              aacaaATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGCCATACTTca        <  1:437898/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGCCATACTTCACAAACG  >  NC_000913/2113217‑2113513

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