Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,018,256 G→C 5.5% T15S (ACC→AGC)  ycbZ ← putative peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,018,2560GC 5.5% 213.7 / 3.3 73T15S (ACC→AGC) ycbZputative peptidase
Reads supporting (aligned to +/- strand):  ref base G (31/38);  new base C (2/2);  total (33/40)
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.53e-01

GGCCTTCGCCAGCATAAAAGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTC  >  NC_000913/1018111‑1018405
                                                                                                                                                 |                                                                                                                                                     
ggCCTTCGCCAGCATAAAAGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATc                                                                                                                                                   <  1:1134059/150‑1 (MQ=255)
   cTTCGCCAGCATAAAAGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCATCGCAAATTGCAGGCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATACCGATAGCTATCAgg                                                                                                                                                <  1:258312/150‑1 (MQ=255)
     tCGCCAGCATAAAAGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGaa                                                                                                                                              >  2:427005/1‑150 (MQ=255)
       gCCAGCATAAAAGAGGCGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGCAcc                                                                                                                                            >  2:779554/1‑150 (MQ=255)
       gCCAGCATAAAAGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAAcc                                                                                                                                            >  2:608535/1‑150 (MQ=255)
        ccAGCATAAAAGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCTAATATTTCCTGATAGCTATCGGTATCAGGAACCa                                                                                                                                           >  2:1080285/1‑150 (MQ=255)
        ccAGCATAAAAGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCa                                                                                                                                           <  2:809913/150‑1 (MQ=255)
        ccAGCATAAAAGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATACCGATAGCTATCAGGaaata                                                                                                                                           <  2:568276/150‑4 (MQ=255)
                 aaGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTCTCGGTATCAGGAACCAGGTCACGcc                                                                                                                                  <  1:1245507/150‑1 (MQ=255)
                  agcggaggaGTCTCGCGTATGCAGCAACTGCTGCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCa                                                                                                                                 <  2:657938/147‑1 (MQ=255)
                       aggaTGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCaa                                                                                                                            >  2:855433/1‑150 (MQ=255)
                       aggaTGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCaa                                                                                                                            <  1:325934/150‑1 (MQ=255)
                          aTGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGtt                                                                                                                         <  1:1299041/150‑1 (MQ=255)
                             cTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCg                                                                                                                      <  1:855433/150‑1 (MQ=255)
                               cgcgTATGCAGCAACTGCTCCAGCGCAAATTGCCGCCGCGGTTGAGGATCACTGAATAAAGGATCGGTTTCGTCAAGCAAATGTGGCTGCGCAAATATTTCCTGATAGCTTTCGGTATCAGGAACCCGGTCACGCCATGCAAGTTTCGTa                                                                                                                    >  2:164967/1‑150 (MQ=255)
                                    aTGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGt                                                                                                               <  1:239185/150‑1 (MQ=255)
                                     tGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTc                                                                                                              <  2:416535/150‑1 (MQ=255)
                                        agcaACTGCTCCAGCGCAAATTGCAGCCGCGGGTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCTAGTTTCGTAATGGTCaaa                                                                                                           <  2:949090/150‑1 (MQ=255)
                                          caACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGt                                                                                                         >  2:1194195/1‑150 (MQ=255)
                                             cTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGa                                                                                                      <  1:1299627/150‑1 (MQ=255)
                                              tGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGAt                                                                                                     <  1:1038247/150‑1 (MQ=255)
                                                 tCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGtt                                                                                                  >  1:515272/1‑150 (MQ=255)
                                                  ccAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGttt                                                                                                 >  2:955990/1‑150 (MQ=255)
                                                   catcgcTAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCATATATTTCCTTATAGCTATCGGTATCAGGATCCAGGTCACGCCATGCAAGTTTCGTAATTGTCAAAGTTGATGtttt                                                                                                <  1:96532/147‑1 (MQ=255)
                                                     gcgcAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGtttttt                                                                                              <  2:283349/150‑1 (MQ=255)
                                                     gcgcAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGtttttt                                                                                              >  1:1051093/1‑150 (MQ=255)
                                                     gcgcAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGtttttt                                                                                              <  2:788199/150‑1 (MQ=255)
                                                      cgcAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTa                                                                                             >  2:698454/1‑150 (MQ=255)
                                                       gcAATTTGCAGACGAGGTTGAGTATCACTGAAAAAAGGATCGTTTTCGTCAATCAAATTTGGATGAGGAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCTAAGTTGATGTTTTTTAg                                                                                            <  2:231450/150‑1 (MQ=255)
                                                       gcAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAg                                                                                            <  2:128703/150‑1 (MQ=255)
                                                           aTTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCtg                                                                                        <  1:655548/150‑1 (MQ=255)
                                                                  ccGCGGATGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCaa                                                                                 <  1:761903/150‑1 (MQ=255)
                                                                         tGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAcgccgcc                                                                          >  1:1155534/1‑149 (MQ=255)
                                                                         tGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGAAAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCcgcg                                                                          <  1:815834/150‑1 (MQ=255)
                                                                           aGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGAt                                                                        >  2:557661/1‑150 (MQ=255)
                                                                            gTATCACTGATTAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGCTGTTTTTTAGTCTGTAGTCAAAGCCGCGAtt                                                                       >  2:203995/1‑150 (MQ=255)
                                                                             tATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGAGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCACCAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTCTTGTCAAAGCCGCGATta                                                                      >  2:318338/1‑150 (MQ=255)
                                                                               tCACTTAATAAAGGATCTTTTTAGTCAATCAAATGTGGCTGAGCATATATTTCCTGATAGCTATCGGTATCATGAACCAGGTCACGCCATGCAAGTTTCGTAATGGACAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATtata                                                                    <  1:217274/150‑1 (MQ=255)
                                                                                cACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATAc                                                                   <  1:31127/150‑1 (MQ=255)
                                                                                cACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATAc                                                                   <  1:226437/150‑1 (MQ=255)
                                                                                   tGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGt                                                                <  2:778576/150‑1 (MQ=255)
                                                                                   tGAATAAAGGATCGTTTTCGTCAATCAAATGTGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTa                                                               <  1:198055/150‑1 (MQ=255)
                                                                                    gAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACAGTa                                                               >  2:1113654/1‑150 (MQ=255)
                                                                                      aTAAAGGATCGTGATCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGGTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAAc                                                             <  2:204593/150‑1 (MQ=255)
                                                                                        aaaGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCg                                                           >  2:1178849/1‑150 (MQ=255)
                                                                                         aaGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATATCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCgg                                                          <  1:874702/150‑1 (MQ=255)
                                                                                                gTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGGTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACa                                                   <  2:737832/150‑1 (MQ=255)
                                                                                                  tttCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGc                                                 <  2:1177703/150‑1 (MQ=255)
                                                                                                      gTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCacac                                             <  2:515272/150‑1 (MQ=255)
                                                                                                      gTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCacac                                             >  1:1132878/1‑150 (MQ=255)
                                                                                                        caatcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGt                                           >  1:871000/1‑150 (MQ=255)
                                                                                                          atcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAg                                         >  1:1076145/1‑150 (MQ=255)
                                                                                                                  gggCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCAcc                                 <  1:291768/149‑1 (MQ=255)
                                                                                                                   ggCTGAGCCAATATTTCCTGATAGCTAGCGGTATCAGGAACCATGTCACGCAATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTATAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCCCTACAGCACACGTAGAAAGCACCg                                >  2:169293/1‑150 (MQ=255)
                                                                                                                     cTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTTATGGTCAAACTTGATGTTTTTTAGTCTGTTGTCAAAGCAGCGATTATACCGAAACCGGCACTACAGCACACGTAGAAAGCACCGAc                              >  2:49757/1‑150 (MQ=255)
                                                                                                                      tGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACa                             >  1:1148942/1‑150 (MQ=255)
                                                                                                                        aGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAAt                           >  2:275504/1‑150 (MQ=255)
                                                                                                                        aGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTATTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGACTGCACCGACAAt                           >  1:801481/1‑150 (MQ=255)
                                                                                                                          cAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATAc                         <  1:608535/150‑1 (MQ=255)
                                                                                                                            aaTATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTc                       >  2:1007958/1‑150 (MQ=255)
                                                                                                                              tatTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCt                     >  2:330727/1‑150 (MQ=255)
                                                                                                                                tttaatGATAGCTTTCGGTATCATGAACCAGTTCACGACATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTTTTGTCAAAGCCGCGTTTATACGGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTgg                   <  2:110573/145‑1 (MQ=255)
                                                                                                                                gttCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTgg                   >  2:1265830/2‑150 (MQ=255)
                                                                                                                                 ttCCTGATACCGATAGCTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAACGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGc                  >  2:690724/1‑150 (MQ=255)
                                                                                                                                  tCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCa                 <  2:520202/150‑1 (MQ=255)
                                                                                                                                   ccTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCAACAATACTCCTGGCAt                >  1:1200774/1‑150 (MQ=255)
                                                                                                                                        tAGCTATCGGTAGCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATTGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCg           <  2:688252/150‑1 (MQ=255)
                                                                                                                                          gCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGCCAATACTCCTGGCATGGGCGtt         >  1:673989/1‑150 (MQ=255)
                                                                                                                                          gCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGtt         >  1:1150207/1‑150 (MQ=255)
                                                                                                                                           cTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTGCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTa        <  1:1123828/150‑1 (MQ=255)
                                                                                                                                            tATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCATCGACAATACTCCTGGCATGGGCGTTaa       >  1:482939/1‑150 (MQ=255)
                                                                                                                                             aTCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGTGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTaaa      <  1:1093741/150‑1 (MQ=255)
                                                                                                                                              ccGGTATCAGGAACCAGGTCACTCCATGCAAGTTTCGTAATGGTCAAAGAAGATGTTTATTAGTCGGGTGTCAAAGCCGCGATTATACCGTAACCCGCACTACAGCACTCGTAGTCAGCACCGACAATACTCCTGGCATGGGCGTTAAAg     <  2:185081/149‑1 (MQ=255)
                                                                                                                                                ggTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAACGCCCCGACAATACTCCTGGCATGGGCGTTAAAGCt   >  2:941800/1‑150 (MQ=255)
                                                                                                                                                 gTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTc  >  1:583588/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GGCCTTCGCCAGCATAAAAGAGGAGGATGCTCGCGTATGCAGCAACTGCTCCAGCGCAAATTGCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTC  >  NC_000913/1018111‑1018405

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