Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,113,373 C→T 5.0% intergenic (‑312/+61) 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,3730CT 5.0% 446.7 / 6.9 138intergenic (‑312/+61)rfbB/wcaNdTDP‑glucose 4,6 dehydratase, NAD(P)‑binding/putative regulatory subunit for GalU
Reads supporting (aligned to +/- strand):  ref base C (61/70);  new base T (2/5);  total (63/75)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 143 total aligned reads displayed.

AAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGCCATACTTCACAAACGCCTGCATA  >  NC_000913/2113224‑2113521
                                                                                                                                                      |                                                                                                                                                    
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)
          aCGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAAtttt                                                                                                                                            >  1:1797359/1‑150 (MQ=255)
           cGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTc                                                                                                                                           <  1:1532156/150‑1 (MQ=255)
             cTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCAc                                                                                                                                         <  2:2042852/150‑1 (MQ=255)
               aaTTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTg                                                                                                                                       >  1:1888974/1‑150 (MQ=255)
                aTTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGTGAAAATTCGCAGCAAAAGTAATTTGTTGCGAATCTTcctgc                                                                                                                                       <  1:2025949/150‑5 (MQ=25)
                aTTGTCTTAACAATCCGCTCTGGAAACAAGGATAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGTGAAACTTCGCAGCAAAAGTAATTTGTTGCGAATCTTcctgc                                                                                                                                       <  1:317746/150‑5 (MQ=25)
                 ttGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGTGAAAATTCGCAGCAAAAGTAATTTGTTGCGAATCTTCCTGcc                                                                                                                                     <  2:2208256/150‑1 (MQ=25)
                 ttGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGcc                                                                                                                                     <  1:1913735/150‑1 (MQ=255)
                  tGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGTGAAAATTCGCAGCAAAAGTAATTTGTTGCGAATCTTCCTGCCg                                                                                                                                    <  2:878492/150‑1 (MQ=25)
                  tGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCg                                                                                                                                    >  2:480438/1‑150 (MQ=255)
                     cTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGttt                                                                                                                                 >  2:382572/1‑149 (MQ=255)
                     cTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTACGAATTTTCACTGCCGTta                                                                                                                                 >  2:529531/1‑150 (MQ=255)
                         ccAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAAt                                                                                                                             >  1:1670446/1‑150 (MQ=255)
                           aaTCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAAttt                                                                                                                           >  1:2113166/1‑150 (MQ=255)
                               cGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCtt                                                                                                                       <  2:719972/150‑1 (MQ=255)
                                gCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTa                                                                                                                      <  1:805525/150‑1 (MQ=255)
                                gCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTa                                                                                                                      <  2:1003629/150‑1 (MQ=255)
                                       aaaCAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCg                                                                                                               <  1:488670/150‑1 (MQ=255)
                                       aaaCAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACAGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCg                                                                                                               <  2:223191/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                                                                                                         >  1:447781/1‑150 (MQ=255)
                                             ggAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACAt                                                                                                         >  2:1413268/1‑150 (MQ=255)
                                                 aaTCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCgg                                                                                                     <  1:517368/150‑1 (MQ=255)
                                                  aTCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGt                                                                                                    <  2:878206/150‑1 (MQ=255)
                                                     cTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTctg                                                                                                 >  2:2021710/1‑148 (MQ=255)
                                                       ggAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGAt                                                                                               >  1:2063449/1‑150 (MQ=255)
                                                         aaaaCTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGAttt                                                                                             >  1:1779475/1‑150 (MQ=255)
                                                         aaaaCTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATAAACCGTTACATCCGGTCAGAttt                                                                                             >  2:204083/1‑150 (MQ=255)
                                                           aaCTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTc                                                                                           <  2:764794/150‑1 (MQ=255)
                                                           aaCTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTc                                                                                           <  2:1601714/150‑1 (MQ=255)
                                                             cTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCat                                                                                         >  2:955164/1‑150 (MQ=255)
                                                              tttGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCatt                                                                                        >  1:535069/1‑150 (MQ=255)
                                                              tttGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCatt                                                                                        >  2:1444252/1‑150 (MQ=255)
                                                               ttGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCatta                                                                                       >  2:1082056/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)
                                                                     aaaaCCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCt                                                                                 <  2:1425017/150‑1 (MQ=255)
                                                                       aaCCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTa                                                                               <  1:1826277/150‑1 (MQ=255)
                                                                       aaCCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTa                                                                               <  2:1189641/150‑1 (MQ=255)
                                                                          ccTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACa                                                                            <  2:1888974/150‑1 (MQ=255)
                                                                          ccTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACa                                                                            <  2:1577351/150‑1 (MQ=255)
                                                                             aCTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCt                                                                         <  1:1686583/150‑1 (MQ=255)
                                                                              cTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCTGTCAGATTTTCATTATTCGCTTAACAGCtt                                                                        >  1:100810/1‑150 (MQ=255)
                                                                                gCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTtct                                                                      <  1:658641/150‑1 (MQ=255)
                                                                                gCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTtct                                                                      >  1:1465092/1‑150 (MQ=255)
                                                                                  tAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCa                                                                    <  1:1033497/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)
                                                                                         ttgttATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTTCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCTGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCtt                                                             <  2:2097417/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:1719576/150‑1 (MQ=255)
                                                                                               ttCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAAATACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGAGTTTCAGTATTCGCTTAACAGCTTCTCATTACCTTTACGGa                                                       >  2:620497/1‑150 (MQ=255)
                                                                                                tCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGaa                                                      >  2:679390/1‑150 (MQ=255)
                                                                                                  tGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTAc                                                          <  1:654481/144‑1 (MQ=255)
                                                                                                  tGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTAc                                                          >  2:654481/1‑144 (MQ=255)
                                                                                                    aTGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTc                                                  >  2:2171468/1‑150 (MQ=255)
                                                                                                      ggTTTATATAAAACAACGGCATGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCAGTGCCGTTATAATTTTCTTATCAACCGTAACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGc                                                <  2:1804861/150‑1 (MQ=255)
                                                                                                       gTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGcc                                               <  2:1454662/150‑1 (MQ=255)
                                                                                                            tataAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCcttg                                          >  1:1707418/1‑149 (MQ=255)
                                                                                                            tataAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCcttc                                          >  1:695871/1‑150 (MQ=255)
                                                                                                             ataAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCcttct                                         >  2:1967223/1‑150 (MQ=255)
                                                                                                             aaattataacGGCAGTGAAAATTCGCAGCAAAAGTAATTTGTTGCGAATCTTCCTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCcttct                                         >  1:2094923/8‑150 (MQ=18)
                                                                                                              taAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGATTTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGACCcttctt                                        >  1:211660/1‑150 (MQ=255)
                                                                                                              taAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAAGAGCTTCTCAATACCTTTACGGAACTTCGCCCcttctt                                        >  2:624426/1‑150 (MQ=255)
                                                                                                                aaaCAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTTACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTc                                      <  1:1413268/150‑1 (MQ=255)
                                                                                                                aaaCAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTc                                      >  1:629909/1‑150 (MQ=255)
                                                                                                                aaaCAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTc                                      >  2:1756516/1‑150 (MQ=255)
                                                                                                                 tataacGGCAGTGAAAATTCGCAGCAAAAGTAATTTGTTGCGAATCTTCCTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCa                                     >  2:1230330/4‑150 (MQ=21)
                                                                                                                         ctggAAGATTCTCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGt                             <  1:299270/148‑1 (MQ=255)
                                                                                                                         cAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCTTTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGt                             <  2:1599319/150‑1 (MQ=255)
                                                                                                                              aaGATTCGCTACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTGCGCCCCTTCTTTCAGGTTGCGTAGGc                         <  1:1153008/150‑1 (MQ=255)
                                                                                                                              aaGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGc                         >  2:972621/1‑150 (MQ=255)
                                                                                                                              aaGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGc                         <  1:786222/150‑1 (MQ=255)
                                                                                                                                  ttCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGCCATa                     >  2:1737027/1‑150 (MQ=255)
                                                                                                                                  ttCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGCCATa                     >  1:478100/1‑150 (MQ=255)
                                                                                                                                    cGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGCCATACt                   >  1:1900151/1‑150 (MQ=255)
                                                                                                                                     gCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGCAATACtt                  >  1:204074/1‑150 (MQ=255)
                                                                                                                                       aacaaATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGCCATACTTca                <  1:437898/150‑1 (MQ=255)
                                                                                                                                        acaaATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATACGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTAAGGTTGCGTAGGCCATACTTcac               >  1:1998848/1‑150 (MQ=255)
                                                                                                                                              tACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGCCATACTTCACAAACGc         <  2:1056720/150‑1 (MQ=255)
                                                                                                                                                 ttttGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGCCATACTTCACAAACGCCTg      <  2:48388/150‑1 (MQ=255)
                                                                                                                                                  tttGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGCCATACTTCACAAACGCCTGc     >  2:2129516/1‑150 (MQ=255)
                                                                                                                                                  tttGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGCCATACTTCACAAACGCCTGc     <  2:1189791/150‑1 (MQ=255)
                                                                                                                                                   ttGCTGCGAATTTTCACTGCCGTTATATTTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGCCATACTTCACAAACGCCTGCa    <  2:509965/150‑1 (MQ=255)
                                                                                                                                                   ttGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCTGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGCCATACTTCACAAACGCCTGCa    <  1:1489820/150‑1 (MQ=255)
                                                                                                                                                     gctgcGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGGCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGCCATACTTCACAAACGCCTGCata  <  1:1077867/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
AAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATCCGGTCAGATTTTCATTATTCGCTTAACAGCTTCTCAATACCTTTACGGAACTTCGCCCCTTCTTTCAGGTTGCGTAGGCCATACTTCACAAACGCCTGCATA  >  NC_000913/2113224‑2113521

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: