Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,477,149 A→T 5.6% S161T (TCA→ACA) ‡ dsdC ← dsd operon activator; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,477,1490AT 5.6% 251.1 / 2.2 87S161T (TCA→ACA) ‡dsdCdsd operon activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base A (43/39);  new base T (3/2);  total (46/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-01

TCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCACCAACCAGCATTGGGCGATAGAG  >  NC_000913/2477001‑2477297
                                                                                                                                                    |                                                                                                                                                    
tCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGAc                                                                                                                                                     <  1:464070/150‑1 (MQ=255)
 cGTTGCTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACg                                                                                                                                                    >  1:785664/1‑150 (MQ=255)
  gTTGCTCCATGCCTGTCTGTCATGGTGCAACGTACTGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGATTTGCGCTGACgg                                                                                                                                                   <  1:901327/150‑1 (MQ=255)
  gTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTTCCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATTCTGGCAGGATTTCTTAATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACgg                                                                                                                                                   >  1:1539625/1‑150 (MQ=255)
  gTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTACTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACgg                                                                                                                                                   <  1:29800/150‑1 (MQ=255)
     gCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc                                                                                                                                                <  1:730369/150‑1 (MQ=255)
     gCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc                                                                                                                                                <  2:1679311/150‑1 (MQ=255)
     gCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc                                                                                                                                                <  2:1672838/150‑1 (MQ=255)
     gCTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc                                                                                                                                                >  2:1684452/1‑150 (MQ=255)
      cTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGca                                                                                                                                               >  1:1233846/1‑150 (MQ=255)
      cTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGca                                                                                                                                               >  1:1496829/1‑150 (MQ=255)
      cTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGca                                                                                                                                               >  2:727707/1‑150 (MQ=255)
      cTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATAAGTGAGTTGCGCTGACGGCGca                                                                                                                                               >  1:110749/1‑150 (MQ=255)
        ccATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATATCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCCGTCAGCGcaac                                                                                                                                             <  1:1409262/150‑3 (MQ=255)
          aTGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGATAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCCGTCAGCGcaactc                                                                                                                                           <  1:1480501/150‑5 (MQ=255)
           tGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCatcatc                                                                                                                                          >  2:1601101/1‑150 (MQ=255)
           tGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCatcatc                                                                                                                                          >  2:189951/1‑150 (MQ=255)
            gCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCAtcatca                                                                                                                                         <  1:651443/150‑1 (MQ=255)
            gCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCAtcatca                                                                                                                                         <  1:657423/150‑1 (MQ=255)
             ccTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCaa                                                                                                                                        >  2:1522295/1‑150 (MQ=255)
                 tctgtcATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTa                                                                                                                                    >  2:1578715/1‑150 (MQ=255)
                      cATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCg                                                                                                                               <  2:667922/150‑1 (MQ=255)
                      cATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCAGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCg                                                                                                                               >  2:1082936/1‑150 (MQ=255)
                            gCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAAt                                                                                                                         >  1:1547482/1‑150 (MQ=255)
                                  tACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTTACGGCGCATCATCAAAGTAGATCGCCAAATCGATTc                                                                                                                   <  1:1516628/150‑1 (MQ=255)
                                  tACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATAAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTc                                                                                                                   <  2:110749/150‑1 (MQ=255)
                                   aCAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTcc                                                                                                                  <  1:937816/150‑1 (MQ=255)
                                   aCAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTcc                                                                                                                  >  2:994375/1‑150 (MQ=255)
                                    cAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCg                                                                                                                 >  1:956757/1‑150 (MQ=255)
                                    cAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCg                                                                                                                 >  2:712471/1‑150 (MQ=255)
                                    cAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCg                                                                                                                 >  1:1506722/1‑150 (MQ=255)
                                     aGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATTTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCgg                                                                                                                <  2:1345231/150‑1 (MQ=255)
                                          cacaGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGt                                                                                                           >  1:1288249/1‑150 (MQ=255)
                                            caGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTg                                                                                                         >  1:230099/1‑150 (MQ=255)
                                                   aTTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTg                                                                                                  <  2:553722/150‑1 (MQ=255)
                                                   aTTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTg                                                                                                  <  1:1322189/150‑1 (MQ=255)
                                                   aTTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATAATTGCAGGATTTCTTCATCCATCAGAAATTGATGAGTGAGTTGCGCTGAAGGCGCAACATCAAAGTAGATCGCAAAATCGATTCCGGCACGTTGCAGAtta                                                                                                  >  1:1001642/1‑149 (MQ=255)
                                                       ccGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGt                                                                                              <  1:1418940/150‑1 (MQ=255)
                                                            ttGGTTAAAGCATGTCTTTGAGCGTATTCCTGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCa                                                                                         <  1:727707/150‑1 (MQ=255)
                                                             tGGTTAAAGCATGTCTTTTAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCAt                                                                                        <  2:209660/150‑1 (MQ=255)
                                                                 tAAAGCAGGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTGCTTCATCCCTCAGAAAGTGATGAGTGAGTGGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTAcc                                                                                    <  2:449554/150‑1 (MQ=255)
                                                                   aaGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAg                                                                                  <  2:1496829/150‑1 (MQ=255)
                                                                       aTGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAg                                                                              <  1:1093507/150‑1 (MQ=255)
                                                                          tCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCAc                                                                           >  2:1641615/1‑150 (MQ=255)
                                                                           cTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACg                                                                          >  2:1210650/1‑150 (MQ=255)
                                                                              tGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTg                                                                       >  2:1719620/1‑150 (MQ=255)
                                                                               gAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGa                                                                      >  1:670977/1‑150 (MQ=255)
                                                                                  cGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCg                                                                   <  1:646687/150‑1 (MQ=255)
                                                                                   gTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGa                                                                  >  2:1703588/1‑150 (MQ=255)
                                                                                     aTTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGaaa                                                                <  2:670977/150‑1 (MQ=255)
                                                                                      ttCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAAt                                                               >  1:1455197/1‑150 (MQ=255)
                                                                                      ttCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCAAGTTGCAGATTGACGTTGTCATTACCAGTGACCACGGTGAGCGAAAt                                                               >  1:965791/1‑150 (MQ=255)
                                                                                       tCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATa                                                              <  2:230099/150‑1 (MQ=255)
                                                                                        ccGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAg                                                             >  2:1632553/1‑150 (MQ=255)
                                                                                          gggCTACTTACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGAc                                                           <  1:1036705/150‑1 (MQ=255)
                                                                                           ggCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACg                                                          <  1:994375/150‑1 (MQ=255)
                                                                                           ggCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACg                                                          <  2:780574/150‑1 (MQ=255)
                                                                                           ggCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACg                                                          <  1:1641615/150‑1 (MQ=255)
                                                                                            gCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACgg                                                         <  1:1662872/150‑1 (MQ=255)
                                                                                                cATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATAc                                                     >  2:24852/1‑150 (MQ=255)
                                                                                                  tACTGGCATGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGAATCATCAAAGTAGATCGCCAAATCGATTCCGGCTCGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGGGAGCGAAATAGACGGATACCg                                                   <  1:816146/150‑1 (MQ=255)
                                                                                                   aCTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCgg                                                  >  2:32269/1‑150 (MQ=255)
                                                                                                   aCTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCgg                                                  >  2:32268/1‑150 (MQ=255)
                                                                                                   aCTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCgg                                                  >  1:818189/1‑150 (MQ=255)
                                                                                                   aCTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCgg                                                  >  2:1257273/1‑150 (MQ=255)
                                                                                                     tGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCg                                                <  1:631790/150‑1 (MQ=255)
                                                                                                          ggATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTaa                                           >  1:759603/1‑150 (MQ=255)
                                                                                                          ggATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGATGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTaa                                           <  1:413281/150‑1 (MQ=255)
                                                                                                               tcttcATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCa                                      >  2:20927/1‑150 (MQ=255)
                                                                                                                   catccatcAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTa                                  <  1:1601101/150‑1 (MQ=255)
                                                                                                                   caaccatcAGAAAGTGATGAGTGAGATCCGATGACGGCCGCTCATCAAAGTAGATCACCAAATCGATTCCGGCACGTTGAAGATTGACGTAGTCATTACCAGTTAGAACGGTGAGCGAAATAGACGGATACCGGCGTGAAAAGTCAccaa                                  >  1:320536/4‑148 (MQ=255)
                                                                                                                     tccatcAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAAt                                >  2:1402464/1‑150 (MQ=255)
                                                                                                                           aGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGc                          >  1:1098038/1‑150 (MQ=255)
                                                                                                                           aGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGc                          <  1:730399/150‑1 (MQ=255)
                                                                                                                            gAATGTGAGGAGTGAGGTGCGCTGACGGCGCATCATCAACGTAGATCGCCAAATCGATTCCGGAACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAATTAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCa                         <  1:334193/150‑1 (MQ=255)
                                                                                                                              aaGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCacc                       <  2:785664/150‑1 (MQ=255)
                                                                                                                                gTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCaccaa                     >  1:956294/1‑150 (MQ=255)
                                                                                                                                    tgagtgagTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCACCAACCAg                 >  2:713775/1‑150 (MQ=255)
                                                                                                                                      ttgatgatGCGCCGTCAGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCCCCTTATGCGGGCCCCAACCAGCa               >  2:711162/8‑150 (MQ=255)
                                                                                                                                      ttgatgatGCGCCGTCAGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCACCAACCAGCa               >  2:538200/8‑150 (MQ=255)
                                                                                                                                         atgatGCGCCGTCAGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCACCAACCAGCATTg            >  2:579705/5‑150 (MQ=255)
                                                                                                                                           gtTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCACCAACCAGCATTggg          <  1:513864/150‑1 (MQ=255)
                                                                                                                                           gtTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCACCAACCAGCATTggg          >  2:1349027/1‑150 (MQ=255)
                                                                                                                                            ttGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCACCAACCAGCATTGGGc         <  1:20927/150‑1 (MQ=255)
                                                                                                                                             tGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCACCAACCAGCATTGGGCg        <  1:1210650/150‑1 (MQ=255)
                                                                                                                                                   gACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTTAATTCACCTAATGCGGGCACCAACCAGCATTGGGCGATagag  >  2:529807/1‑150 (MQ=255)
                                                                                                                                                   gACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCACCAACCAGCATTGGGCGATagag  >  2:1661379/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCACCAACCAGCATTGGGCGATAGAG  >  NC_000913/2477001‑2477297

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