Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,282,202 C→T 6.0% A76V (GCC→GTC)  agaS → tagatose‑6‑phosphate ketose/aldose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,282,2020CT 6.0% 260.4 / 2.9 82A76V (GCC→GTC) agaStagatose‑6‑phosphate ketose/aldose isomerase
Reads supporting (aligned to +/- strand):  ref base C (35/42);  new base T (2/3);  total (37/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATCCGTTCACTCACCAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGTCGATCCGGCAACAGCCCGGAAAGCGTCGCAGCCGTGGAACTG  >  NC_000913/3282055‑3282350
                                                                                                                                                   |                                                                                                                                                    
ggATCCGTTCACTCACCAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCgc                                                                                                                                                      >  2:165761/1‑148 (MQ=255)
 gATCCGTTCACTCACCAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACGGGAGCCGGAACGTCGGCATTTATCGGTGTCATCATCGCGCCGTGGCTCgcca                                                                                                                                                    <  4:163984/149‑1 (MQ=255)
 gATCCGTTCACTCACCAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCgc                                                                                                                                                      >  7:115570/1‑147 (MQ=255)
 gATCCGTTCACTCACCAACATCGACGCGCTACGTTCCGCGCGCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCgcca                                                                                                                                                    <  1:102979/149‑1 (MQ=255)
  aTCCGTTCACTCACCAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGTACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCgccag                                                                                                                                                   <  4:18458/149‑1 (MQ=255)
   tCCGTTCACTCACCAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCgccagc                                                                                                                                                  >  7:244560/1‑149 (MQ=255)
      gTTCACTCACCAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCgccagcca                                                                                                                                                >  1:278740/1‑148 (MQ=255)
            tcacCAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCgg                                                                                                                                          >  6:63343/1‑148 (MQ=255)
              acaaaCATCGACGCGCTACGTTCCGCGCTCAATAACTTCCGTGAACCGGTACAGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCAGGGACATCATCGCGCCGTGGCTCCCCAGCCAGAACGGTa                                                                                                                                        >  5:121562/4‑148 (MQ=255)
              acCAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCATCCATACCGGt                                                                                                                                         <  7:422052/147‑1 (MQ=255)
              acCAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTa                                                                                                                                        >  7:95579/1‑148 (MQ=255)
              acCAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACAGTTACTGCGCAAAGAGAATCGGCGGATCATCCTGCCCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGTCAGCCATACCGGTa                                                                                                                                        >  7:293625/1‑148 (MQ=255)
                cAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTaaa                                                                                                                                      >  5:347427/1‑148 (MQ=255)
                cAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTaaa                                                                                                                                      >  2:136695/1‑148 (MQ=255)
                cAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTaaa                                                                                                                                      >  3:349334/1‑148 (MQ=255)
                cAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTGGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTaaa                                                                                                                                      >  5:8109/1‑148 (MQ=255)
                 aaCATCGACGCGCGACGGTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGACGCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTaaaaa                                                                                                                                    <  3:315594/149‑1 (MQ=255)
                    aTCGACGCGCTACGTTCCACGATCAATAACGTCCTTGAACAGTTACTTCGCAAAGAGAATCTGCGGATCATCCTCACCGGAGCCGGACCGTCGGCATTTATCGGTGCCATCATCGCGCCGTGGCTCGCCAGCCATCCCGGGAAAAACtt                                                                                                                                 >  6:281634/1‑149 (MQ=255)
                     tCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTc                                                                                                                                <  6:88328/149‑1 (MQ=255)
                     tCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTc                                                                                                                                <  8:173044/149‑1 (MQ=255)
                      cGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACtt                                                                                                                                 >  7:315071/1‑147 (MQ=255)
                       gACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCa                                                                                                                               <  4:247422/148‑1 (MQ=255)
                       gACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCa                                                                                                                               <  2:352265/148‑1 (MQ=255)
                       gACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAg                                                                                                                              >  4:324376/1‑149 (MQ=255)
                         cgcgATACGTTCCGCGCTCTATAACTTCCTTTAACCGTTACTGCGCAAAGAGAATCTGCGGAGCATCCTGACCTGAGCCTGAACGTCGGCAGTTATCGGTGACATCATCGCGCCGTGGATCGCCAGCAATACCGGTAAAAACTTc                                                                                                                                <  5:281634/145‑1 (MQ=255)
                           cgcTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGcc                                                                                                                          <  6:167607/149‑1 (MQ=255)
                             cTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCg                                                                                                                         >  1:189911/1‑148 (MQ=255)
                             cTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGt                                                                                                                        <  3:230313/149‑1 (MQ=255)
                              tACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGAGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTa                                                                                                                       >  1:313110/1‑149 (MQ=255)
                                cGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTa                                                                                                                       <  2:278740/147‑1 (MQ=255)
                                cGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTAcc                                                                                                                     <  5:204073/149‑1 (MQ=255)
                                        gcTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCg                                                                                                             <  8:318577/149‑1 (MQ=255)
                                         cTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGa                                                                                                            >  1:51527/1‑149 (MQ=255)
                                          tCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGAt                                                                                                           >  3:399471/1‑149 (MQ=255)
                                          tCAATAACTTCCTTGAACCGTTACTGCGCAAAAAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCAGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGAt                                                                                                           >  6:45412/1‑149 (MQ=255)
                                            aaTAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGAt                                                                                                           <  4:349334/147‑1 (MQ=255)
                                             aTAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTg                                                                                                        <  3:56283/149‑1 (MQ=255)
                                              tAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGac                                                                                                                  <  4:250116/138‑1 (MQ=255)
                                              tAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGac                                                                                                                  >  3:250116/1‑138 (MQ=255)
                                                 cTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCa                                                                                                    <  7:166862/149‑1 (MQ=255)
                                                        gAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGCCCGGAGCCGGAACGACGGCATTTATCGGTGACATGATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATcc                                                                                             >  2:8755/1‑149 (MQ=255)
                                                        gAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATc                                                                                              <  1:177059/148‑1 (MQ=255)
                                                        gAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATc                                                                                              <  8:244560/148‑1 (MQ=255)
                                                           ccGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGAt                                                                                          <  4:310787/149‑1 (MQ=255)
                                                                    cgcAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTAcc                                                                                  <  6:347427/148‑1 (MQ=255)
                                                                    cgcAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTAc                                                                                   <  4:175113/147‑1 (MQ=255)
                                                                       aaaGAGTATGTGCGGATAATACTGACCGGAGGAGTAACGTGGGCATTTATCGGTGACATCATTGCGCAGTGGTTCTTCTGACATACCGGTAAAAACTTAATCTTCGTACTGACCACCGATCTTGTCAACTTTCCTATGTACTACCt                                                                                 <  6:121562/146‑1 (MQ=255)
                                                                       aaaGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGaa                                                                              <  3:324376/149‑1 (MQ=255)
                                                                         agagAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAAcc                                                                            <  2:189911/149‑1 (MQ=255)
                                                                             aaTCTGCGGATCATCCTGCCCGGAGCCGGAACGTCGGCATTTATCGGTGACATCCTCGCGCCGTGGCTCGCCAGCCATACCGGTAACAACTTCAGCGCCGTCCCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAg                                                                         >  2:244333/1‑148 (MQ=255)
                                                                             aaTCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGc                                                                        <  1:40326/149‑1 (MQ=255)
                                                                                   cGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGACACCAATCCGATGGACTACCTGAACCCAGCTCATcc                                                                  >  2:298163/1‑149 (MQ=255)
                                                                                      atcatcCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCgct                                                               >  1:238873/1‑149 (MQ=255)
                                                                                            cTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGtt                                                         <  8:303672/149‑1 (MQ=255)
                                                                                             tGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGtt                                                         >  6:371362/1‑148 (MQ=255)
                                                                                              gACCGGTGCCGGAACGGCGGCAGTTAGCGGTGACATCATCGCGCCGTGGCTCGTCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGa                                                       <  8:293625/149‑1 (MQ=255)
                                                                                                   gAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCGAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTcc                                                  <  2:97359/149‑1 (MQ=255)
                                                                                                   gAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGaa                                                                              <  7:60405/121‑1 (MQ=255)
                                                                                                   gAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGaa                                                                              >  8:60405/1‑121 (MQ=255)
                                                                                                   gAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGAtct                                                    >  4:136468/1‑147 (MQ=255)
                                                                                                   gAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGAtct                                                    <  3:136468/147‑1 (MQ=255)
                                                                                                     gCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCtt                                                >  3:358591/1‑149 (MQ=255)
                                                                                                      ccGGAACGTCGTCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCTGTATAAACTTCAGCGCCGTACCGACCTCCGATCTGGTCACCTATCCGATGTACTAGCTGAAGCCAGCTCATCCGCTGCTGTTGATCTCCTTc                                               <  1:110147/149‑1 (MQ=255)
                                                                                                        ggAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCgg                                             <  2:209219/149‑1 (MQ=255)
                                                                                                         gAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGt                                            <  8:95579/149‑1 (MQ=255)
                                                                                                         gAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGt                                            >  3:78435/1‑149 (MQ=255)
                                                                                                             gTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGAAGGACTACCTGAACCCAGCTCATACGCTGCTCTTGATCTCCTTCGGTCGa                                         >  1:321884/1‑148 (MQ=255)
                                                                                                                ggCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGTCGATcc                                      <  1:165761/148‑1 (MQ=255)
                                                                                                                ggCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGTCGATCCg                                     <  5:371362/149‑1 (MQ=255)
                                                                                                                   aTTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGTCGATCCGGc                                   >  5:172336/1‑148 (MQ=255)
                                                                                                                    tttATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGGCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGTCGATCCGGCaa                                 <  1:136695/149‑1 (MQ=255)
                                                                                                                       aTCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGTCGATCCGGCAACAg                              >  4:53251/1‑149 (MQ=255)
                                                                                                                               catcatCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGTCGATCCGGCAACAGCCCGGaaa                      <  8:165209/149‑1 (MQ=255)
                                                                                                                                atcatcGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGTCGATCCGGCAACAGCCCGGAAAg                     >  7:231665/1‑149 (MQ=255)
                                                                                                                                 tcatcGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGAAGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGTCGATCCGGCAACAGCCCGGAAAGc                    <  5:311556/149‑1 (MQ=255)
                                                                                                                                   atcGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGTCGATCCGGCAACAGCCCGGAAAGCGt                  >  5:293170/1‑149 (MQ=255)
                                                                                                                                    tcGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGTCGATCCGGCAACAGCCCGGAAAGCGTc                 <  2:51527/149‑1 (MQ=255)
                                                                                                                                    tcGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGTCGATCCGGCAACAGCCCGGAAAGCGTc                 <  4:78435/149‑1 (MQ=255)
                                                                                                                                    tcGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGTCGATCCGGCAACAGCCCGGAAAGCGTc                 <  8:231665/149‑1 (MQ=255)
                                                                                                                                      gcgcCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGTCGATCCGGCAACAGCCCGGAAAGCGTCg                <  3:53251/148‑1 (MQ=255)
                                                                                                                                        gcCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGTCGATCCGGCAACAGCCCGGAAAGCGTCGCAg             >  4:264050/1‑149 (MQ=255)
                                                                                                                                            tGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGTCGATCCGGCATCAGCCCGGAAAGCGTCGCAGCCGt         <  8:115570/149‑1 (MQ=255)
                                                                                                                                            tGGCTCGCCAGCCATACAGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGTCGATCCGGCAACAGCCCGGAAAGCGTCGCAGCCGt         <  4:399471/149‑1 (MQ=255)
                                                                                                                                              gCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGACGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGTCGATCCGGCAACAGCCCGGAAAGCGTCGCAGCCGTg        >  4:275096/1‑148 (MQ=255)
                                                                                                                                                tCGTCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGTCGATCCGGCAACAGCCCGGAAAGCGTCGCAGCCGt         >  4:388250/1‑145 (MQ=255)
                                                                                                                                                tCGTCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGTCGATCCGGCAACAGCCCGGAAAGCGTCGCAGCCGt         <  3:388250/145‑1 (MQ=255)
                                                                                                                                                tCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGCACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGTCGATCCGGCAACAGCCCGGAAAGCGTCGCAGCCGTGGaa     >  7:389079/1‑149 (MQ=255)
                                                                                                                                                tCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGTCGATCCGGCCACAGCCCGGAAAGCGTCGCAGCCGTAGaa     >  7:372040/1‑149 (MQ=255)
                                                                                                                                                tCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGTCGATCCGGCAACAGCCCGGAAAGCGTCGCAGCCGTGGaa     >  1:164934/1‑149 (MQ=255)
                                                                                                                                                  gccatccaAACCGGTAAAAACTTCAGCGCCGTACCGACCCCCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGTCGATCCGGCAACAGCCCGGAAAGCGTCGCAGCCGTGGaa     <  6:212604/147‑1 (MQ=255)
                                                                                                                                                  gccagccaTACCGGTATAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACGAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGTCGATCCGGCAACAGCCCGGAAAGCGGCGCAGCCGTGGAAc    <  6:418685/148‑1 (MQ=255)
                                                                                                                                                  gccagccaTACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGTCGATCCGGCAACAGCCCGGAAAGCGTCGCAGCCGTGGAACt   <  5:45412/149‑1 (MQ=255)
                                                                                                                                                   ccagccaTACCGGTAAAAACTTCAGCGCCGTACAGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGTCGATCCGGCAACAGCCCGGAAAGCGTCGCAGCCGTGGAACTg  >  2:39761/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGATCCGTTCACTCACCAACATCGACGCGCTACGTTCCGCGCTCAATAACTTCCTTGAACCGTTACTGCGCAAAGAGAATCTGCGGATCATCCTGACCGGAGCCGGAACGTCGGCATTTATCGGTGACATCATCGCGCCGTGGCTCGCCAGCCATACCGGTAAAAACTTCAGCGCCGTACCGACCACCGATCTGGTCACCAATCCGATGGACTACCTGAACCCAGCTCATCCGCTGCTGTTGATCTCCTTCGGTCGATCCGGCAACAGCCCGGAAAGCGTCGCAGCCGTGGAACTG  >  NC_000913/3282055‑3282350

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: ‑