Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,113,359 A→G 5.1% intergenic (‑298/+75) 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,3590AG 5.1% 431.6 / 9.0 137intergenic (‑298/+75)rfbB/wcaNdTDP‑glucose 4,6 dehydratase, NAD(P)‑binding/putative regulatory subunit for GalU
Reads supporting (aligned to +/- strand):  ref base A (64/66);  new base G (2/5);  total (66/71)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

Only 100 of 142 total aligned reads displayed.

CGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGCCATACTTCAC  >  NC_000913/2113211‑2113508
                                                                                                                                                     |                                                                                                                                                     
cgcgCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGcaac                                                                                                                                                      >  1:509965/1‑150 (MQ=255)
cgcgCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGcaac                                                                                                                                                      >  1:1911542/1‑150 (MQ=255)
      aaaCCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTa                                                                                                                                                >  1:398907/1‑150 (MQ=255)
         cccGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACtt                                                                                                                                             >  2:1652316/1‑150 (MQ=255)
          ccGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACttt                                                                                                                                            >  1:856145/1‑150 (MQ=255)
          ccGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACttt                                                                                                                                            >  1:1532204/1‑150 (MQ=255)
           cGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACtttt                                                                                                                                           >  1:764794/1‑150 (MQ=255)
           cGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACtttt                                                                                                                                           <  1:2052353/150‑1 (MQ=255)
             aaaaaaTCAAACGCTAATTGTCTTAACAATCCGCGCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAATACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTtgc                                                                                                                                         <  2:339584/150‑1 (MQ=255)
             aaaaaTTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTtgc                                                                                                                                         <  2:2097047/150‑1 (MQ=255)
              aaaaTTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTtgct                                                                                                                                        >  2:1826277/1‑150 (MQ=255)
                 aTTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCg                                                                                                                                     >  2:1790502/1‑150 (MQ=255)
                     aaaCGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGTAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAAtt                                                                                                                                 <  1:1412500/150‑1 (MQ=255)
                        cGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTc                                                                                                                              <  1:1532156/150‑1 (MQ=255)
                          cTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCAc                                                                                                                            <  2:2042852/150‑1 (MQ=255)
                             aTTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGTGAAAATTCGCAGCAAAAGTAATTTGTTGCGAATCTTcctgc                                                                                                                          <  1:2025949/150‑5 (MQ=25)
                               tGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGTGAAAATTCGCAGCAAAAGTAATTTGTTGCGAATCTTC‑CTGCCg                                                                                                                       <  2:878492/150‑1 (MQ=25)
                               tGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCg                                                                                                                       >  2:1033497/1‑150 (MQ=255)
                                  cTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTACGAATTTTCACTGCCGTta                                                                                                                    >  2:529531/1‑150 (MQ=255)
                                      ccAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAAt                                                                                                                >  1:1670446/1‑150 (MQ=255)
                                      ccAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAAt                                                                                                                >  1:456957/1‑150 (MQ=255)
                                            cGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCtt                                                                                                          <  2:719972/150‑1 (MQ=255)
                                               tctGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATc                                                                                                       >  2:111247/1‑150 (MQ=255)
                                                    aagcAAGGAAAATCCTGGAAAACATTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACTGCAG‑GAAGATTCTCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCg                                                                                                  <  1:109450/147‑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                                                                                            >  2:1413268/1‑150 (MQ=255)
                                                          ggAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACAt                                                                                            >  1:447781/1‑150 (MQ=255)
                                                              aaTCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCgg                                                                                        <  1:517368/150‑1 (MQ=255)
                                                               aTCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGt                                                                                       <  2:878206/150‑1 (MQ=255)
                                                                tCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTc                                                                                      >  2:437898/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)
                                                                        aaCTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTc                                                                              <  2:764794/150‑1 (MQ=255)
                                                                         aCTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCa                                                                             >  1:1466060/1‑150 (MQ=255)
                                                                          cTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCat                                                                            >  2:955164/1‑150 (MQ=255)
                                                                           tttGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCatt                                                                           >  1:535069/1‑150 (MQ=255)
                                                                           tttGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCatt                                                                           >  2:1444252/1‑150 (MQ=255)
                                                                               aaTAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTc                                                                       <  1:1760062/150‑1 (MQ=255)
                                                                                 tAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCAGGTCAGATTTTCATTATTCGc                                                                     >  2:253565/1‑150 (MQ=255)
                                                                                  aaaaCCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCt                                                                    <  2:301575/150‑1 (MQ=255)
                                                                                  aaaaCCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCt                                                                    <  2:398907/150‑1 (MQ=255)
                                                                                  aaaaCCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGGTACATCCGGTCAGATTTTCATTATTCGCt                                                                    <  2:15951/150‑1 (MQ=255)
                                                                                    aaCCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTa                                                                  <  2:1189641/150‑1 (MQ=255)
                                                                                    aaCCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTa                                                                  <  1:1826277/150‑1 (MQ=255)
                                                                                       ccTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACa                                                               <  2:1577351/150‑1 (MQ=255)
                                                                                       ccTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACa                                                               <  2:1888974/150‑1 (MQ=255)
                                                                                           cTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCTGTCAGATTTTCATTATTCGCTTAACAGCtt                                                           >  1:100810/1‑150 (MQ=255)
                                                                                             gCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTtct                                                         >  1:848167/1‑150 (MQ=255)
                                                                                             gCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTtct                                                         <  1:658641/150‑1 (MQ=255)
                                                                                             gCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTtct                                                         >  1:1465092/1‑150 (MQ=255)
                                                                                             gCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATAAACCGTTACATACGGTCAGATTTTCATTATTCGCTTAACAGCTtct                                                         >  1:1988222/1‑150 (MQ=255)
                                                                                               tAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCa                                                       <  1:1033497/150‑1 (MQ=255)
                                                                                               tAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCa                                                       <  1:382572/150‑1 (MQ=255)
                                                                                               tAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCa                                                       <  1:1239167/150‑1 (MQ=255)
                                                                                                 aCTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAAt                                                     >  2:616994/1‑150 (MQ=255)
                                                                                                 aCTCGTTGTTAATCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAAt                                                     <  2:1042632/150‑1 (MQ=255)
                                                                                                  cTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATa                                                    >  1:652357/1‑150 (MQ=255)
                                                                                                   tCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATAc                                                   <  1:1678719/150‑1 (MQ=255)
                                                                                                   tCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATAc                                                   <  2:2121996/150‑1 (MQ=255)
                                                                                                        gttATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTa                                              >  1:1012077/1‑150 (MQ=255)
                                                                                                          tATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCCTCTCAATACCTTTACg                                            <  1:1627642/150‑1 (MQ=255)
                                                                                                           aTTCTGATTGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACgg                                           <  1:635250/150‑1 (MQ=255)
                                                                                                           aTTCTGATGGTTTATATATAACAACTGCAG‑GAAGATTAGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTTTAATTTTCTTATCAACCGTTACATCCGGTCAGTTTTTGATTATTCGGTTAACAGGTTCTCAATACCTTTACgg                                           <  2:154518/150‑1 (MQ=255)
                                                                                                           aTTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACgg                                           <  1:2117852/150‑1 (MQ=255)
                                                                                                           aTTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACgg                                           <  1:1719576/150‑1 (MQ=255)
                                                                                                            ttCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAAATACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGAGTTTCAGTATTCGCTTAACAGCTTCTCATTACCTTTACGGa                                          >  2:620497/1‑150 (MQ=255)
                                                                                                               tGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTAc                                             >  2:654481/1‑144 (MQ=255)
                                                                                                                 aTGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTc                                     >  2:2171468/1‑150 (MQ=255)
                                                                                                                 aTGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTc                                     <  1:1249617/150‑1 (MQ=255)
                                                                                                                  tGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCg                                    >  2:352719/1‑150 (MQ=255)
                                                                                                                   ggTTTATATAAAACAACGGCAT‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCAGTGCCGTTATAATTTTCTTATCAACCGTAACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGc                                   <  2:1804861/150‑1 (MQ=255)
                                                                                                                    gTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGcc                                  <  2:1454662/150‑1 (MQ=255)
                                                                                                                     ttttTATAAAACAACGGAAG‑GAAGATTCGAAAGATTTTACTTTTGCTGCGAATTTTAACTGCAGTTATAATTTTCTTATCAACCGTTACATCCGTTCAGATTTTCATTATTCGCTTAACAGCTTGTCAATACCTTTAGGGATCttatccc                                 <  2:1424474/150‑6 (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)
                                                                                                                          aaattataacGGCAGTGAAAATTCGCAGCAAAAGTAATTTGTTGCGAATCTTC‑CTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCcttct                            >  1:2094923/8‑150 (MQ=18)
                                                                                                                           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)
                                                                                                                                      ctg‑gAAGATTCTCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGt                <  1:299270/148‑1 (MQ=255)
                                                                                                                                           aaGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGc            <  1:786222/150‑1 (MQ=255)
                                                                                                                                           aaGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGc            >  2:972621/1‑150 (MQ=255)
                                                                                                                                            aGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGcc           <  2:1206678/150‑1 (MQ=255)
                                                                                                                                               ttCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGCCATa        >  1:478100/1‑150 (MQ=255)
                                                                                                                                               ttCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGCCATa        >  2:1737027/1‑150 (MQ=255)
                                                                                                                                                 cGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGCCATACt      >  1:1900151/1‑150 (MQ=255)
                                                                                                                                                  gCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGCAATACtt     >  1:204074/1‑150 (MQ=255)
                                                                                                                                                     acaaATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATACGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTAAGGTTGCGTAGGCCATACTTcac  >  1:1998848/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAG‑GAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGCCATACTTCAC  >  NC_000913/2113211‑2113508

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: ‑