Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,058,653 T→C 6.0% intergenic (‑243/+13) serA ← / ← rpiA D‑3‑phosphoglycerate dehydrogenase/ribose 5‑phosphate isomerase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,058,6530TC 6.0% 327.3 / 9.5 117intergenic (‑243/+13)serA/rpiAD‑3‑phosphoglycerate dehydrogenase/ribose 5‑phosphate isomerase, constitutive
Reads supporting (aligned to +/- strand):  ref base T (51/58);  new base C (3/4);  total (54/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.02e-01

Only 100 of 118 total aligned reads displayed.

CGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGATTT  >  NC_000913/3058506‑3058802
                                                                                                                                                   |                                                                                                                                                     
cgcgATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTc                                                                                                                                                     <  1:243309/150‑1 (MQ=255)
 gcgATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTcc                                                                                                                                                    >  1:1122625/1‑150 (MQ=255)
     tGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCg                                                                                                                                                <  1:1666342/150‑1 (MQ=255)
     tGTTGAGGAAATACGACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGTGGATTTATTTAGAGCAATCGACATTTGCCTGGTAAAAGCGTGACACATGTCACGAAATTTAATGAAGAGAATTTTTTTAACGGGGGTGGTTCCCCCg                                                                                                                                                <  2:802860/150‑1 (MQ=255)
     tGTTGAGGAAATACCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCg                                                                                                                                                <  2:802855/150‑1 (MQ=255)
       ttGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTc                                                                                                                                              <  2:478815/150‑1 (MQ=255)
       ttGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAACCTCCCCCGtt                                                                                                                                              <  1:12881/150‑2 (MQ=255)
        tGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCa                                                                                                                                             >  1:1015900/1‑150 (MQ=255)
        tGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCa                                                                                                                                             <  1:783472/150‑1 (MQ=255)
        tGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAACCTCCCCCGtta                                                                                                                                             <  1:1237090/150‑3 (MQ=255)
           ggAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAACCTCCCCCGttaaaa                                                                                                                                          <  2:334175/150‑6 (MQ=255)
            gAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAACCTCCCCCGttaaaaa                                                                                                                                         <  1:1495277/150‑7 (MQ=255)
             aaaTATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCa                                                                                                                                        >  2:1130826/1‑150 (MQ=255)
                tatCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCAttt                                                                                                                                     <  1:621804/150‑1 (MQ=255)
                 atCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTc                                                                                                                                    <  2:544745/150‑1 (MQ=255)
                 atCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTc                                                                                                                                    >  2:1260148/1‑150 (MQ=255)
                 atCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTc                                                                                                                                    <  1:1250540/150‑1 (MQ=255)
                     cGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACaa                                                                                                                                >  2:36409/1‑150 (MQ=255)
                       cAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATg                                                                                                                              >  2:830177/1‑150 (MQ=255)
                        aaaGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGGATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATgg                                                                                                                             <  1:952868/150‑1 (MQ=255)
                           ggTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGttt                                                                                                                          <  1:791380/150‑1 (MQ=255)
                            gTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGtttt                                                                                                                         <  1:741665/150‑1 (MQ=255)
                             tAATTCTTGTATGAATATGATGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACTAAATTTAATGAAGAGAATTTTTTTAACGTTGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTg                                                                                                                        <  2:174279/150‑1 (MQ=255)
                             tAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTg                                                                                                                        >  1:1146803/1‑150 (MQ=255)
                             tAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTg                                                                                                                        <  2:1677045/150‑1 (MQ=255)
                             tAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTg                                                                                                                        >  1:1317762/1‑150 (MQ=255)
                               aTTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGac                                                                                                                      <  2:1130135/150‑1 (MQ=255)
                               aTTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGac                                                                                                                      <  2:1416476/150‑1 (MQ=255)
                                    tGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGt                                                                                                                 <  1:374965/150‑1 (MQ=255)
                                      gATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCTTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTATTGACACCGTCa                                                                                                               >  2:1696389/1‑150 (MQ=255)
                                      gATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCa                                                                                                               >  1:882590/1‑150 (MQ=255)
                                      gATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCa                                                                                                               >  2:592958/1‑150 (MQ=255)
                                      gATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCa                                                                                                               >  1:232982/1‑150 (MQ=255)
                                        tGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAgg                                                                                                             <  1:420000/150‑1 (MQ=255)
                                              tGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGcc                                                                                                       <  1:387196/150‑1 (MQ=255)
                                              tGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGcc                                                                                                       >  1:164510/1‑150 (MQ=255)
                                              tGCTGGTTTAGAGGATNTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGcc                                                                                                       >  1:1558860/1‑150 (MQ=255)
                                               gCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGGTCCCCCGTCAGATCATTTCACAATGGTTATGACACCGTCAGGTGTGCCa                                                                                                      <  2:74078/150‑1 (MQ=255)
                                                cTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAAt                                                                                                    >  1:407357/1‑150 (MQ=255)
                                                  ggTTTAGATGATTTATTTATAGCAATCGACAATTTCCTTGTAAAAGCGTGACACATGTCACCAAATTTTATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCTTTTCACAATGGTTTTGACACCGTCAGTTGTGCCAATc                                                                                                   <  1:655535/150‑1 (MQ=255)
                                                  ggTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGa                                                                                                                       >  2:1614771/1‑130 (MQ=255)
                                                  ggTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGa                                                                                                                       <  1:1614771/130‑1 (MQ=255)
                                                  ggTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATc                                                                                                   <  1:1419811/150‑1 (MQ=255)
                                                   gTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCa                                                                                                  <  1:746222/150‑1 (MQ=255)
                                                       agagGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCGATCAgcgc                                                                                              >  1:397900/1‑150 (MQ=255)
                                                       agagGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAgcgc                                                                                              >  1:979268/1‑150 (MQ=255)
                                                       agagGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAgcgc                                                                                              >  2:1575834/1‑150 (MQ=255)
                                                        gagGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACAGTCAGGTGTGCCAATCAgcgaa                                                                                             >  1:1040672/1‑148 (MQ=255)
                                                         agGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCaa                                                                                            >  1:1089636/1‑150 (MQ=255)
                                                         agGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCaa                                                                                            <  1:21709/150‑1 (MQ=255)
                                                          ggATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAAc                                                                                           <  2:1385711/150‑1 (MQ=255)
                                                          ggATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAAc                                                                                           >  1:1238220/1‑150 (MQ=255)
                                                          ggATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGACAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAAc                                                                                           >  2:215474/1‑150 (MQ=255)
                                                           gATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACg                                                                                          <  2:1286262/150‑1 (MQ=255)
                                                             tttatttaGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTc                                                                                        <  1:472673/150‑1 (MQ=255)
                                                               tatttaGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCcg                                                                                      >  2:244159/1‑150 (MQ=255)
                                                                          aTCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTa                                                                           <  2:1122625/150‑1 (MQ=255)
                                                                          aTCGACAATTGCCTGGTAAAAGCGTGACACATGTAACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTa                                                                           >  1:1426398/1‑150 (MQ=255)
                                                                           tCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAg                                                                          <  2:1238220/150‑1 (MQ=255)
                                                                           tCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAg                                                                          <  2:899191/150‑1 (MQ=255)
                                                                            cGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGCGGTTCACCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAga                                                                         >  1:710813/1‑149 (MQ=255)
                                                                            cGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCATGTGTGCCAATCAGAGCAACGTCCGCGCCACGGTTAGc                                                                         >  1:241632/1‑150 (MQ=255)
                                                                             gACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCa                                                                        <  2:1523812/150‑1 (MQ=255)
                                                                               cAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCaaa                                                                      >  2:275284/1‑150 (MQ=255)
                                                                                 aTTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACa                                                                    <  1:382148/150‑1 (MQ=255)
                                                                                      cTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCa                                                               >  1:911017/1‑150 (MQ=255)
                                                                                      cTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCa                                                               >  1:1646815/1‑150 (MQ=255)
                                                                                          tAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAACCAAGCCAACAg                                                           <  2:1695158/150‑1 (MQ=255)
                                                                                           aatttggTGACATGTGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGt                                                          >  1:799493/7‑150 (MQ=255)
                                                                                               gCGGGACACATGTCACCAAATTTAATGAAGAGAATTTTTATAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACTAGCCAACAGTcacc                                                      <  1:72878/150‑1 (MQ=255)
                                                                                                  tGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACg                                                   >  2:433001/1‑150 (MQ=255)
                                                                                                         tGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAAt                                           <  2:407357/150‑1 (MQ=255)
                                                                                                           tCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATc                                          >  2:887399/1‑150 (MQ=255)
                                                                                                            cACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCg                                         <  2:1317762/150‑1 (MQ=255)
                                                                                                              ccAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCa                                       <  2:1089636/150‑1 (MQ=255)
                                                                                                               cAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCAt                                      >  2:465831/1‑150 (MQ=255)
                                                                                                                aaaTTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCAtt                                     <  1:376013/150‑1 (MQ=255)
                                                                                                                aaaTTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCAtt                                     <  1:1641969/150‑1 (MQ=255)
                                                                                                                    ttAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATg                                 <  2:1015900/150‑1 (MQ=255)
                                                                                                                     tAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACTCCAGGAATCGCATTTATgg                                <  1:41580/150‑1 (MQ=255)
                                                                                                                       aTGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCg                              <  1:433001/150‑1 (MQ=255)
                                                                                                                       aTGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCg                              <  2:1558860/150‑1 (MQ=255)
                                                                                                                             cgaATTTTTTTAACGGGTGAGGTTCCCCCGGCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCCACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGTATCGCATTTATGGCGTTTTcc                        <  1:9981/149‑1 (MQ=255)
                                                                                                                             agaATTTTTTTAACGGGGTAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCTACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCAAGCCAGGAATCGCATTTATGGCGTTTTcc                        <  1:15939/150‑1 (MQ=255)
                                                                                                                             agaATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTcc                        >  2:14335/1‑150 (MQ=255)
                                                                                                                                  ttttttAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGc                   <  1:1402687/150‑1 (MQ=255)
                                                                                                                                      ttAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATc               <  1:1130826/150‑1 (MQ=255)
                                                                                                                                      ctgaCGGGGGAACCTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATc               >  1:956644/4‑150 (MQ=255)
                                                                                                                                      ctgaCGGGGGAACCTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGAAAACAAGCCAACAGTCACCCCGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATc               >  1:433955/4‑150 (MQ=255)
                                                                                                                                       tAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCg              >  2:693352/1‑150 (MQ=255)
                                                                                                                                        gaCGGGGGAACCTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCACCGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCgg             >  1:1114620/2‑150 (MQ=255)
                                                                                                                                        aaCGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCgg             <  2:232982/150‑1 (MQ=255)
                                                                                                                                        aaCGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCgg             <  2:882590/150‑1 (MQ=255)
                                                                                                                                         aCGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCggg            <  2:1136556/150‑1 (MQ=255)
                                                                                                                                           gggggAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGATTCGCATTTATGGCGTTTTCCAGAGCTATCGGGTc          >  2:877727/1‑150 (MQ=255)
                                                                                                                                           gggggAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTc          >  2:900675/1‑150 (MQ=255)
                                                                                                                                            ggggAGGTGCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGGTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCa         <  1:174409/150‑1 (MQ=255)
                                                                                                                                               gAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAgg      <  2:1646815/150‑1 (MQ=255)
                                                                                                                                                aGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGa     >  1:984173/1‑150 (MQ=255)
                                                                                                                                                   ttCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGAttt  <  2:1146803/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGATTT  >  NC_000913/3058506‑3058802

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: