Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,058,652 G→C 5.8% intergenic (‑242/+14) 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,6520GC 5.8% 399.7 / 9.8 117intergenic (‑242/+14)serA/rpiAD‑3‑phosphoglycerate dehydrogenase/ribose 5‑phosphate isomerase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (51/59);  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 = 8.25e-01

Only 100 of 119 total aligned reads displayed.

CGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGA  >  NC_000913/3058503‑3058799
                                                                                                                                                     |                                                                                                                                                   
cGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACTCATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAgg                                                                                                                                                     <  2:1209983/150‑1 (MQ=255)
cGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAgg                                                                                                                                                     <  2:190244/150‑1 (MQ=255)
    gcgATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTcc                                                                                                                                                 >  1:1122625/1‑150 (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)
            gAGGAAATATCACGCAAAGGTAATTCTTGGAGGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAg                                                                                                                                         <  2:1175392/150‑1 (MQ=255)
             aGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGa                                                                                                                                        >  1:53137/1‑150 (MQ=255)
              ggAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGAt                                                                                                                                       >  1:749236/1‑150 (MQ=255)
              ggAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGAt                                                                                                                                       <  2:586/150‑1 (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                                                                                                                                 <  1:1250540/150‑1 (MQ=255)
                    atCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTc                                                                                                                                 >  2:1260148/1‑150 (MQ=255)
                    atCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTc                                                                                                                                 <  2:544745/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)
                                tAATTCTTGTATGAATATGATGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACTAAATTTAATGAAGAGAATTTTTTTAACGTTGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTg                                                                                                                     <  2:174279/150‑1 (MQ=255)
                                tAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTg                                                                                                                     >  1:1146803/1‑150 (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:232982/1‑150 (MQ=255)
                                           tGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAgg                                                                                                          <  1:420000/150‑1 (MQ=255)
                                              atatGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGgtgt                                                                                                       <  1:947745/150‑1 (MQ=255)
                                                 tGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGcc                                                                                                    <  1:387196/150‑1 (MQ=255)
                                                 tGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGcc                                                                                                    >  1:164510/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)
                                                         tAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAgcg                                                                                            <  1:140203/150‑1 (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:21709/150‑1 (MQ=255)
                                                            agGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCaa                                                                                         >  1:1089636/1‑150 (MQ=255)
                                                             ggATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAAc                                                                                        <  2:1385711/150‑1 (MQ=255)
                                                             ggATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAAc                                                                                        >  1:1238220/1‑150 (MQ=255)
                                                               atttatttAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGt                                                                                      >  2:376013/1‑150 (MQ=255)
                                                                tttatttaGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTc                                                                                     <  1:472673/150‑1 (MQ=255)
                                                                  tatttaGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCcg                                                                                   >  2:244159/1‑150 (MQ=255)
                                                                      taGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGcc                                                                               >  1:408457/1‑150 (MQ=255)
                                                                             aTCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTa                                                                        >  1:346168/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:174409/1‑150 (MQ=255)
                                                                                  cAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCaaa                                                                   >  2:275284/1‑150 (MQ=255)
                                                                                    aTTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACa                                                                 <  1:382148/150‑1 (MQ=255)
                                                                                         cTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCa                                                            >  1:1646815/1‑150 (MQ=255)
                                                                                         cTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCa                                                            >  1:911017/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)
                                                                                                          caTGTCACCAAATTTTATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAgg                                           >  2:1402687/1‑150 (MQ=255)
                                                                                                             gTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAAt                                        >  2:1645460/1‑150 (MQ=255)
                                                                                                              tCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATc                                       >  2:887399/1‑150 (MQ=255)
                                                                                                              tCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATc                                       >  2:1028485/1‑150 (MQ=255)
                                                                                                               cACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCg                                      <  2:1317762/150‑1 (MQ=255)
                                                                                                                 ccAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCa                                    <  2:1089636/150‑1 (MQ=255)
                                                                                                                  cAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCAt                                   >  1:1579899/1‑150 (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                           <  2:1558860/150‑1 (MQ=255)
                                                                                                                          aTGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCg                           <  1:433001/150‑1 (MQ=255)
                                                                                                                            gAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGtt                         >  1:1702380/1‑150 (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)
                                                                                                                                          tAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCg           >  2:693352/1‑150 (MQ=255)
                                                                                                                                           gaCGGGGGAACCTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCACCGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCgg          >  1:1114620/2‑150 (MQ=255)
                                                                                                                                           aaCGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCgg          <  2:882590/150‑1 (MQ=255)
                                                                                                                                           aaCGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCgg          <  2:232982/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)
                                                                                                                                                     |                                                                                                                                                   
CGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGA  >  NC_000913/3058503‑3058799

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: