Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,477,147 T→C 5.9% S161S (TCA→TCG) ‡ dsdC ← dsd operon activator; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,477,1470TC 5.9% 249.3 / 5.1 85S161S (TCA→TCG) ‡dsdCdsd operon activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (41/39);  new base C (3/2);  total (44/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 = 8.88e-01

AGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCACCAACCAGCATTGGGCG  >  NC_000913/2476999‑2477291
                                                                                                                                                    |                                                                                                                                                
aGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTTATTACCCTGTTGGTTAAAGCATGTCATTGAGCGTAATACGGGCTACATACTGGCAGGATTTATTCATACATAAGAAAGTGATGAGTGAGTTGCGCtt                                                                                                                                                 >  2:1093507/1‑149 (MQ=255)
  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:657423/150‑1 (MQ=255)
              gCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCAtcatca                                                                                                                                   <  1:651443/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                                                                                            <  1:1322189/150‑1 (MQ=255)
                                                     aTTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTg                                                                                            <  2:553722/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:1257273/1‑150 (MQ=255)
                                                                                                     aCTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCgg                                            >  1:818189/1‑150 (MQ=255)
                                                                                                     aCTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCgg                                            >  2:32268/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)
                                                                                                                                                    |                                                                                                                                                
AGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTAGATCGCCAAATCGATTCCGGCACGTTGCAGATTGACGTTGTCATTACCAGTGAGCACGGTGAGCGAAATAGACGGATACCGGCGTGTAAAGTCACCTAATGCGGGCACCAACCAGCATTGGGCG  >  NC_000913/2476999‑2477291

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: