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.

CGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGCCATACTTCAC  >  NC_000913/2113211‑2113508
                                                                                                                                                     |                                                                                                                                                     
cgcgCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGcaac                                                                                                                                                      >  1:509965/1‑150 (MQ=255)
cgcgCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGcaac                                                                                                                                                      >  1:1911542/1‑150 (MQ=255)
      aaaCCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTa                                                                                                                                                >  1:398907/1‑150 (MQ=255)
         cccGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACtt                                                                                                                                             >  2:1652316/1‑150 (MQ=255)
          ccGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACttt                                                                                                                                            >  1:856145/1‑150 (MQ=255)
          ccGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACttt                                                                                                                                            >  1:1532204/1‑150 (MQ=255)
           cGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACtttt                                                                                                                                           >  1:764794/1‑150 (MQ=255)
           cGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACtttt                                                                                                                                           <  1:2052353/150‑1 (MQ=255)
             aaaaaaTCAAACGCTAATTGTCTTAACAATCCGCGCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAATACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTtgc                                                                                                                                         <  2:339584/150‑1 (MQ=255)
             aaaaaTTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTtgc                                                                                                                                         <  2:2097047/150‑1 (MQ=255)
              aaaaTTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTtgct                                                                                                                                        >  2:1826277/1‑150 (MQ=255)
                 aTTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCg                                                                                                                                     >  2:1790502/1‑150 (MQ=255)
                     aaaCGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGTAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAAtt                                                                                                                                 <  1:1412500/150‑1 (MQ=255)
                        cGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTc                                                                                                                              <  1:1532156/150‑1 (MQ=255)
                          cTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCAc                                                                                                                            <  2:2042852/150‑1 (MQ=255)
                             aTTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGTGAAAATTCGCAGCAAAAGTAATTTGTTGCGAATCTTcctgc                                                                                                                          <  1:2025949/150‑5 (MQ=25)
                               tGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGTGAAAATTCGCAGCAAAAGTAATTTGTTGCGAATCTTCCTGCCg                                                                                                                       <  2:878492/150‑1 (MQ=25)
                               tGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCg                                                                                                                       >  2:1033497/1‑150 (MQ=255)
                                  cTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTACGAATTTTCACTGCCGTta                                                                                                                    >  2:529531/1‑150 (MQ=255)
                                      ccAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAAt                                                                                                                >  1:1670446/1‑150 (MQ=255)
                                      ccAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAAt                                                                                                                >  1:456957/1‑150 (MQ=255)
                                            cGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCtt                                                                                                          <  2:719972/150‑1 (MQ=255)
                                               tctGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATc                                                                                                       >  2:111247/1‑150 (MQ=255)
                                                    aagcAAGGAAAATCCTGGAAAACATTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACTGCAGGAAGATTCTCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCg                                                                                                  <  1:109450/147‑1 (MQ=255)
                                                    aaaCAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCg                                                                                                  <  1:1857242/150‑1 (MQ=255)
                                                    aaaCAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACAGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCg                                                                                                  <  2:223191/150‑1 (MQ=255)
                                                      aCAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGtt                                                                                                <  2:497386/150‑1 (MQ=255)
                                                       cAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTa                                                                                               <  1:906498/150‑1 (MQ=255)
                                                       cAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTa                                                                                               >  1:1056720/1‑150 (MQ=255)
                                                       aaaGGAAAATCCTTGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTa                                                                                               <  2:491441/149‑1 (MQ=255)
                                                        aaGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTAc                                                                                              >  1:561840/1‑150 (MQ=255)
                                                          ggAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACAt                                                                                            >  2:1413268/1‑150 (MQ=255)
                                                          ggAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACAt                                                                                            >  1:447781/1‑150 (MQ=255)
                                                              aaTCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCgg                                                                                        <  1:517368/150‑1 (MQ=255)
                                                               aTCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGt                                                                                       <  2:878206/150‑1 (MQ=255)
                                                                tCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTc                                                                                      >  2:437898/1‑150 (MQ=255)
                                                                      aaaaCTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGAttt                                                                                >  1:1779475/1‑150 (MQ=255)
                                                                      aaaaCTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATAAACCGTTACATCCGGTCAGAttt                                                                                >  2:204083/1‑150 (MQ=255)
                                                                        aaCTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTc                                                                              <  2:1601714/150‑1 (MQ=255)
                                                                        aaCTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTc                                                                              <  2:764794/150‑1 (MQ=255)
                                                                         aCTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCa                                                                             >  1:1466060/1‑150 (MQ=255)
                                                                          cTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCat                                                                            >  2:955164/1‑150 (MQ=255)
                                                                           tttGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCatt                                                                           >  1:535069/1‑150 (MQ=255)
                                                                           tttGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCatt                                                                           >  2:1444252/1‑150 (MQ=255)
                                                                               aaTAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTc                                                                       <  1:1760062/150‑1 (MQ=255)
                                                                                 tAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCAGGTCAGATTTTCATTATTCGc                                                                     >  2:253565/1‑150 (MQ=255)
                                                                                  aaaaCCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCt                                                                    <  2:301575/150‑1 (MQ=255)
                                                                                  aaaaCCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCt                                                                    <  2:398907/150‑1 (MQ=255)
                                                                                  aaaaCCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGGTACATCCGGTCAGATTTTCATTATTCGCt                                                                    <  2:15951/150‑1 (MQ=255)
                                                                                    aaCCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTa                                                                  <  2:1189641/150‑1 (MQ=255)
                                                                                    aaCCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTa                                                                  <  1:1826277/150‑1 (MQ=255)
                                                                                       ccTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACa                                                               <  2:1577351/150‑1 (MQ=255)
                                                                                       ccTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACa                                                               <  2:1888974/150‑1 (MQ=255)
                                                                                           cTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCTGTCAGATTTTCATTATTCGCTTAACAGCtt                                                           >  1:100810/1‑150 (MQ=255)
                                                                                             gCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTtct                                                         >  1:848167/1‑150 (MQ=255)
                                                                                             gCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTtct                                                         <  1:658641/150‑1 (MQ=255)
                                                                                             gCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTtct                                                         >  1:1465092/1‑150 (MQ=255)
                                                                                             gCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATAAACCGTTACATACGGTCAGATTTTCATTATTCGCTTAACAGCTtct                                                         >  1:1988222/1‑150 (MQ=255)
                                                                                               tAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCa                                                       <  1:1033497/150‑1 (MQ=255)
                                                                                               tAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCa                                                       <  1:382572/150‑1 (MQ=255)
                                                                                               tAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCa                                                       <  1:1239167/150‑1 (MQ=255)
                                                                                                 aCTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAAt                                                     >  2:616994/1‑150 (MQ=255)
                                                                                                 aCTCGTTGTTAATCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAAt                                                     <  2:1042632/150‑1 (MQ=255)
                                                                                                  cTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATa                                                    >  1:652357/1‑150 (MQ=255)
                                                                                                   tCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATAc                                                   <  1:1678719/150‑1 (MQ=255)
                                                                                                   tCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATAc                                                   <  2:2121996/150‑1 (MQ=255)
                                                                                                        gttATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTa                                              >  1:1012077/1‑150 (MQ=255)
                                                                                                          tATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCCTCTCAATACCTTTACg                                            <  1:1627642/150‑1 (MQ=255)
                                                                                                           aTTCTGATTGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACgg                                           <  1:635250/150‑1 (MQ=255)
                                                                                                           aTTCTGATGGTTTATATATAACAACTGCAGGAAGATTAGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTTTAATTTTCTTATCAACCGTTACATCCGGTCAGTTTTTGATTATTCGGTTAACAGGTTCTCAATACCTTTACgg                                           <  2:154518/150‑1 (MQ=255)
                                                                                                           aTTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACgg                                           <  1:2117852/150‑1 (MQ=255)
                                                                                                           aTTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACgg                                           <  1:1719576/150‑1 (MQ=255)
                                                                                                            ttCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAAATACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGAGTTTCAGTATTCGCTTAACAGCTTCTCATTACCTTTACGGa                                          >  2:620497/1‑150 (MQ=255)
                                                                                                               tGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTAc                                             >  2:654481/1‑144 (MQ=255)
                                                                                                                 aTGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTc                                     >  2:2171468/1‑150 (MQ=255)
                                                                                                                 aTGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTc                                     <  1:1249617/150‑1 (MQ=255)
                                                                                                                  tGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCg                                    >  2:352719/1‑150 (MQ=255)
                                                                                                                   ggTTTATATAAAACAACGGCATGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCAGTGCCGTTATAATTTTCTTATCAACCGTAACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGc                                   <  2:1804861/150‑1 (MQ=255)
                                                                                                                    gTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGcc                                  <  2:1454662/150‑1 (MQ=255)
                                                                                                                     ttttTATAAAACAACGGAAGGAAGATTCGAAAGATTTTACTTTTGCTGCGAATTTTAACTGCAGTTATAATTTTCTTATCAACCGTTACATCCGTTCAGATTTTCATTATTCGCTTAACAGCTTGTCAATACCTTTAGGGATCttatccc                                 <  2:1424474/150‑6 (MQ=255)
                                                                                                                         tataAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCcttg                             >  1:1707418/1‑149 (MQ=255)
                                                                                                                         tataAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCcttc                             >  1:695871/1‑150 (MQ=255)
                                                                                                                          ataAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCcttct                            >  2:1967223/1‑150 (MQ=255)
                                                                                                                          ataAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCcttct                            <  1:1790502/150‑1 (MQ=255)
                                                                                                                          aaattataacGGCAGTGAAAATTCGCAGCAAAAGTAATTTGTTGCGAATCTTCCTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCcttct                            >  1:2094923/8‑150 (MQ=18)
                                                                                                                           taAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGATTTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGACCcttctt                           >  1:211660/1‑150 (MQ=255)
                                                                                                                           taAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAAGAGCTTCTCAATACCTTTACGGAACTTCGCCCcttctt                           >  2:624426/1‑150 (MQ=255)
                                                                                                                           taAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCcttctt                           >  2:680259/1‑150 (MQ=255)
                                                                                                                            aaaaCAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCttt                          <  1:2071060/150‑1 (MQ=255)
                                                                                                                             aaaCAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTc                         >  2:1756516/1‑150 (MQ=255)
                                                                                                                             aaaCAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTc                         >  1:629909/1‑150 (MQ=255)
                                                                                                                                      ctggAAGATTCTCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGt                <  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)
                                                                                                                                                     |                                                                                                                                                     
CGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGCCATACTTCAC  >  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: