Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,222,725 T→C 17.7% R31R (CGT→CGC ebgA → evolved beta‑D‑galactosidase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,222,7250TC17.7% 1942.6 / inf 811R31R (CGT→CGCebgAevolved beta‑D‑galactosidase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (339/328);  new base C (69/75);  total (408/403)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 819 total aligned reads displayed.

TTACACAAACTACCATCACCATGAATGGTTCCGATTTCTCTCTACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGG  >  NC_000913/3222577‑3222872
                                                                                                                                                    |                                                                                                                                                   
ttACACAAACTACCATCACCATGAATGGTTCCGATTTCTCTCTACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGt                                                                                                                                                     >  1:6337633/1‑149 (MQ=255)
   cacaAACTACCATCACCATGAATGGTTCCGATTTCTCTCTACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTAcc                                                                                                                                                  >  4:3267722/1‑149 (MQ=255)
   cacaAACTACCATCACCATGAATGGTTCCGATTTCTCTCTACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTAcc                                                                                                                                                  >  1:6115564/1‑149 (MQ=255)
   cacaAACTACCATCACCATGAATGGTTCCGATTTCTCTCTACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGCAcc                                                                                                                                                  <  4:2777551/149‑1 (MQ=255)
      aaaCTACCATCACCATGAATGGTTCCGATTTCTCTCTACCGGGCGGCCCTATGAATCGCTGGGCAAACATTCTGCTCACCCACGAAACCCGTCTTGCGCCGCGCGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGACAGTACCttt                                                                                                                                               >  4:4857068/1‑149 (MQ=255)
      aaaCTACCATCACCATGAATGGTTCCGATTTCTCTCTACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCttt                                                                                                                                               <  2:5534064/149‑1 (MQ=255)
       aaCTACCATCACCATGAATGGTTCCGATTTCTCTCTACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTg                                                                                                                                              >  4:4617937/1‑149 (MQ=255)
        aCTACCATCACCATGAATGGTTCCGATTTCTCTCTACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGc                                                                                                                                             <  1:4181994/149‑1 (MQ=255)
             cATCACCATGAATGGTTCCGATTTCTCTCTACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCcgcg                                                                                                                                        >  3:2692878/1‑149 (MQ=255)
                cACCATGAATGGTTCCGATTTCTCTCTACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGaaa                                                                                                                                     >  1:1712133/1‑149 (MQ=255)
                  ccATGAATGGTTCCGATTTCTCTCTACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAAcc                                                                                                                                   >  4:1051338/1‑149 (MQ=255)
                  ccATGAATGGTTCCGATTTCTCTCTACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGCACCTTTGCCCGCGAAAcc                                                                                                                                   >  3:6456040/1‑149 (MQ=255)
                   cATGAATGGTTCCGATTTCTCTCTACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACca                                                                                                                                  <  4:1697148/149‑1 (MQ=255)
                    aTGAATGGTTCCGATTTCTCTCTACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGGAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATACTGTTGCCCAAGCGCGTACCTTTGCCCGCGAAACcag                                                                                                                                 >  4:6892914/1‑149 (MQ=255)
                    aTGAATGGTTCCGATTTCTCTCTACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACcag                                                                                                                                 >  4:572737/1‑149 (MQ=255)
                     tGAATGGTTCCGATTTCTCTCTACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACcagc                                                                                                                                >  4:4278326/1‑149 (MQ=255)
                      gAATGGTTCCGATTTCTCTCTACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAAc                                                                                                                                    >  1:2391098/1‑144 (MQ=255)
                       aaTGGTTCCGATTTCTCTCTACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACcagcag                                                                                                                              >  1:7145886/1‑149 (MQ=255)
                            ttCCGATTTCTCTCTACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGt                                                                                                                         >  4:6252517/1‑149 (MQ=255)
                            ttCCGATTTCTCTCTACCGAGAGGCCCAATGAATGGCTGGGAAATCATTCAGCTCACGCACGACAACCGACTTGCGTCGCGTGCGCGCTTTTTTTCATATGATTCTGTTGCGCACGCGCCTACCTTTGCCCGCGAAGCCAGCCGCgtgt                                                                                                                         >  2:3316616/1‑149 (MQ=255)
                             tCCGATTTCTCTCTACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGtt                                                                                                                        <  1:853246/149‑1 (MQ=255)
                              ccGATTTCTCTCTACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACcagca                                                                                                                               >  4:2939163/1‑141 (MQ=255)
                              ccGATTTCTCTCTACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGttt                                                                                                                       >  1:365282/1‑149 (MQ=255)
                                    tctctcTACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGccc                                                                                                                 <  1:5110251/149‑1 (MQ=255)
                                     ctctctACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCt                                                                                                                <  1:1858523/149‑1 (MQ=255)
                                     ctctctACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCt                                                                                                                <  1:264185/149‑1 (MQ=255)
                                       ctctACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTa                                                                                                              >  3:3837522/1‑149 (MQ=255)
                                        tctACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTaa                                                                                                             >  1:2985837/1‑149 (MQ=255)
                                         ctACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAg                                                                                                            >  3:5142884/1‑149 (MQ=255)
                                          tACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGCACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGc                                                                                                           <  4:71346/149‑1 (MQ=255)
                                           aCCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCg                                                                                                          <  4:5076156/149‑1 (MQ=255)
                                              gggAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGATAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTc                                                                                                       <  3:6139757/149‑1 (MQ=255)
                                                gAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGt                                                                                                        <  2:3653568/146‑1 (MQ=255)
                                                 aGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGt                                                                                                    <  2:2292638/149‑1 (MQ=255)
                                                   gCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTgg                                                                                                  <  1:385367/149‑1 (MQ=255)
                                                     ccTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGccc                                                                                                                                           <  2:1294358/106‑1 (MQ=255)
                                                     ccTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGaa                                                                                                <  3:1087316/149‑1 (MQ=255)
                                                     ccTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGaa                                                                                                <  4:3657117/149‑1 (MQ=255)
                                                      cTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGCACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAAt                                                                                               >  4:1910507/1‑149 (MQ=255)
                                                       tATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAAtt                                                                                              <  1:4777860/149‑1 (MQ=255)
                                                       tATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAAtt                                                                                              >  2:5842087/1‑149 (MQ=255)
                                                           aaTCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCa                                                                                          <  3:4974857/149‑1 (MQ=255)
                                                                cTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACtttt                                                                                     >  4:846608/1‑149 (MQ=255)
                                                                cTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGCACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACtttt                                                                                     <  3:2432273/149‑1 (MQ=255)
                                                                 tGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACttttt                                                                                    >  3:2141059/1‑149 (MQ=255)
                                                                  gggAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACtttttt                                                                                   <  2:3942655/149‑1 (MQ=255)
                                                                   ggAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTg                                                                                  >  2:3987895/1‑149 (MQ=255)
                                                                    gAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGa                                                                                 >  2:7202312/1‑149 (MQ=255)
                                                                    gAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGa                                                                                 <  2:1869033/149‑1 (MQ=255)
                                                                    gAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGACTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGa                                                                                 >  4:1203329/1‑149 (MQ=255)
                                                                     aaaaCATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGAc                                                                                <  2:6298137/149‑1 (MQ=255)
                                                                     aaaaCATTCAGCACACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGAc                                                                                >  3:5461363/1‑149 (MQ=255)
                                                                               gCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGCACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGc                                                                      >  3:2979001/1‑149 (MQ=255)
                                                                                     ccACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTAc                                                                >  1:1053089/1‑149 (MQ=255)
                                                                                     ccACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTAc                                                                >  3:3634805/1‑149 (MQ=255)
                                                                                     ccACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGCACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTAc                                                                <  3:6619646/149‑1 (MQ=255)
                                                                                     ccACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGCACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTAc                                                                <  2:2518225/149‑1 (MQ=255)
                                                                                      cACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTAcc                                                               >  2:6153502/1‑149 (MQ=255)
                                                                                       aCGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCa                                                              >  2:893666/1‑149 (MQ=255)
                                                                                          aaaaCCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGaa                                                           <  2:2735655/149‑1 (MQ=255)
                                                                                           aaaCCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAACCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAg                                                          >  2:5991569/1‑149 (MQ=255)
                                                                                            aaCCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGCTTCTGTTGCGCAAGCGCGCACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATGCGCTGCAAGTACCAGAAGc                                                         <  2:6684260/149‑1 (MQ=255)
                                                                                             aCCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGCACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGcc                                                        <  1:6280341/149‑1 (MQ=255)
                                                                                               cGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCa                                                                                          <  4:3826252/113‑1 (MQ=255)
                                                                                                gACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTc                                                     <  4:5421215/149‑1 (MQ=255)
                                                                                                 aCTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCa                                                    <  4:6580007/149‑1 (MQ=255)
                                                                                                        ccgcGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGa                                             <  4:1337974/149‑1 (MQ=255)
                                                                                                        ccgcGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGa                                             <  2:4580634/149‑1 (MQ=255)
                                                                                                        ccgcGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGCACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAGGTACCAGAAGCCTTCACCTCTGa                                             <  1:6940248/149‑1 (MQ=255)
                                                                                                          gcgtgcgtACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGt                                           <  3:468086/149‑1 (MQ=255)
                                                                                                          gcgtgcgtACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGt                                           <  1:5899563/149‑1 (MQ=255)
                                                                                                           cgtgcgtACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGtt                                          <  3:828675/149‑1 (MQ=255)
                                                                                                               cgtACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTGTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATg                                      <  1:447462/149‑1 (MQ=255)
                                                                                                                gtACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTGCCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATgg                                     <  2:1070212/149‑1 (MQ=255)
                                                                                                                gtACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATgg                                     >  2:4976878/1‑149 (MQ=255)
                                                                                                                  aCTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCt                                   <  4:2110191/149‑1 (MQ=255)
                                                                                                                   cTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTg                                  >  3:3307236/1‑149 (MQ=255)
                                                                                                                    tttttttCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGa                                 >  1:5440716/1‑149 (MQ=255)
                                                                                                                    tttttttCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGa                                 >  2:430177/1‑149 (MQ=255)
                                                                                                                    tttttttCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGa                                 >  3:1970428/1‑149 (MQ=255)
                                                                                                                    tttttttCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGa                                 >  4:2485177/1‑149 (MQ=255)
                                                                                                                    tttttttCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTCACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGa                                 >  4:5959508/1‑149 (MQ=255)
                                                                                                                    tttttttCATATGATTCTGTTGCGCAAGCGCGCACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGaa                                                           >  1:2208408/1‑123 (MQ=255)
                                                                                                                     ttttttCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGAc                                <  1:1233152/149‑1 (MQ=255)
                                                                                                                     ttttttCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGAc                                <  2:2964646/149‑1 (MQ=255)
                                                                                                                          tCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTgggg                           >  1:1471027/1‑149 (MQ=255)
                                                                                                                             tatGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCa                        >  2:1501086/1‑149 (MQ=255)
                                                                                                                             tatGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCa                        >  3:153371/1‑149 (MQ=255)
                                                                                                                             tatGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCa                        >  2:6873174/1‑149 (MQ=255)
                                                                                                                               tGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCata                      <  4:401883/149‑1 (MQ=255)
                                                                                                                               tGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCata                      >  1:3442834/1‑149 (MQ=255)
                                                                                                                                gATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCatat                     >  3:1842596/1‑149 (MQ=255)
                                                                                                                                    cTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTAcc                 <  3:4306823/149‑1 (MQ=255)
                                                                                                                                     tGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCg                >  2:2102041/1‑149 (MQ=255)
                                                                                                                                             aaGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTc                                                     >  1:2786524/1‑104 (MQ=255)
                                                                                                                                             aaGCGCGCACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGcc        >  3:6913492/1‑149 (MQ=255)
                                                                                                                                                cgcgTACCGTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCAtg     <  2:1681998/149‑1 (MQ=255)
                                                                                                                                                 gcgTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCAtgt    <  1:6698803/149‑1 (MQ=255)
                                                                                                                                                 gcgTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCAtgt    <  3:1373285/149‑1 (MQ=255)
                                                                                                                                                   gTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTgg  >  3:5770206/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTACACAAACTACCATCACCATGAATGGTTCCGATTTCTCTCTACCGGGAGGCCCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGG  >  NC_000913/3222577‑3222872

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: