Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,072,411 C→T 5.3% A81T (GCC→ACC)  yihT ← 6‑deoxy‑6‑sulphofructose‑1‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,072,4110CT 5.3% 244.1 / 3.0 76A81T (GCC→ACC) yihT6‑deoxy‑6‑sulphofructose‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (39/32);  new base T (2/2);  total (42/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.28e-01

ACCAGCAGTTTTAAGGCTTTGCCACCGTCTTGTTTGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGTTAACTTTGAAATCAGTTAAAACGCTATCGGCTACCGGAGCGGGTGCCCCAGCCGCGG  >  NC_000913/4072272‑4072557
                                                                                                                                           |                                                                                                                                                  
aCCAGCAGTTTTAAGGCTTTGCCACCGTCTTGTTTGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGtt                                                                                                                                           <  1:246216/149‑1 (MQ=255)
  cagcagTTTTAAGGCTTTGCCACCGTCTTGTTTGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGtttt                                                                                                                                         >  3:359851/1‑149 (MQ=255)
    gcagTTTTAAGGCTTTGCCACCGTCTTGTTTGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGc                                                                                                                                       <  4:280886/149‑1 (MQ=255)
    gcagTTTTAAGGCTTTGCCACCGTCTTGTTTGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGc                                                                                                                                       >  4:189442/1‑149 (MQ=255)
    gcagTTTTAAGGCTTTGCCACCGTCTTGTTTGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGc                                                                                                                                       <  1:313166/149‑1 (MQ=255)
    gcagTTTTAAGGCTTTGCCACCGTCTTGGTTGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACCTTCGGCAATCGCGTTTTGc                                                                                                                                       >  2:33696/1‑149 (MQ=255)
        ttttAAGGCTTTGCCACCGTCTTGTTTGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCaa                                                                                                                                   >  2:235233/1‑149 (MQ=255)
         tttAAGGCTTTGCCACCGTCTTGTTTGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAac                                                                                                                                  <  1:280222/149‑1 (MQ=255)
         tttAAGGCTTTGCCACCGTCTTGTTTGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAac                                                                                                                                  <  3:235781/149‑1 (MQ=255)
              ggCTTTGCCACCGTCTTGTTTGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCt                                                                                                                             >  3:200053/1‑149 (MQ=255)
               gCTTTGCCACCGTCTTGTTTGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTg                                                                                                                            >  1:219017/1‑149 (MQ=255)
                cTTTGCCACCGTCTTGTTTGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTgg                                                                                                                           <  1:446533/149‑1 (MQ=255)
                cTTTGCCACCGTCTTGTTTGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTgg                                                                                                                           >  4:537526/1‑149 (MQ=255)
                cTTTGCCACCGTCTTGTTTGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTgg                                                                                                                           <  3:16449/149‑1 (MQ=255)
                 tttGCCACCGTCTTGTTTGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGc                                                                                                                          >  3:551825/1‑149 (MQ=255)
                 tttGCCACCGTCTTGTTTGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGc                                                                                                                          >  3:98832/1‑149 (MQ=255)
                 tttGCCACCGTCTTGTTTGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGc                                                                                                                          >  1:132796/1‑149 (MQ=255)
                     ccACCGTCTTGTTTGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTa                                                                                                                      <  2:40669/149‑1 (MQ=255)
                        ccGTCTTGTTTGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCa                                                                                                                   <  4:165771/149‑1 (MQ=255)
                             ttGTTTGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAAtt                                                                                                              >  3:106239/1‑149 (MQ=255)
                             ttGTTTGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAAtt                                                                                                              >  1:243157/1‑149 (MQ=255)
                                tttGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAAttgtt                                                                                                           <  2:219017/149‑1 (MQ=255)
                                tttGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAAttgtt                                                                                                           >  2:321543/1‑149 (MQ=255)
                                tttGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAAttgtt                                                                                                           >  4:19523/1‑149 (MQ=255)
                                 ttGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGTAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAAttgttg                                                                                                          >  4:133051/1‑149 (MQ=255)
                                 ttGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGTAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAAttgttg                                                                                                          >  4:96342/1‑149 (MQ=255)
                                  tGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGa                                                                                                         <  2:308949/149‑1 (MQ=255)
                                  tGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGa                                                                                                         <  2:243288/149‑1 (MQ=255)
                                   gATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGAt                                                                                                        <  2:447072/149‑1 (MQ=255)
                                    aTCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATc                                                                                                       >  1:576604/1‑149 (MQ=255)
                                      cTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTa                                                                                                     <  2:6611/149‑1 (MQ=255)
                                        gTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTAcc                                                                                                   <  4:106239/149‑1 (MQ=255)
                                                 tgattTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGcc                                                                                          <  2:279111/149‑1 (MQ=255)
                                                     tttGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAgg                                                                                      <  2:67576/149‑1 (MQ=255)
                                                          ggTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATaa                                                                                 <  3:463638/149‑1 (MQ=255)
                                                          ggTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATaa                                                                                 >  2:454575/1‑149 (MQ=255)
                                                                  aCCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGagag                                                                         <  1:358711/149‑1 (MQ=255)
                                                                  aCCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGagag                                                                         >  1:471867/1‑149 (MQ=255)
                                                                         tATCGACCGGAATACCGTTGCCAGGAATGAACTCCTCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCttt                                                                  >  2:323363/1‑149 (MQ=255)
                                                                         tATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCttt                                                                  <  4:343322/149‑1 (MQ=255)
                                                                         tATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAACCGCCGAGCCATAACGCTAGAGGGCCttt                                                                  <  2:269600/149‑1 (MQ=255)
                                                                              aCCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAAttgt                                                                                                            <  3:123151/102‑1 (MQ=255)
                                                                              aCCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAAttgt                                                                                                            >  4:123151/1‑102 (MQ=255)
                                                                                cGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGt                                                           >  2:459071/1‑149 (MQ=255)
                                                                                      aCCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGTTAACtt                                                     >  1:283529/1‑149 (MQ=255)
                                                                                         gTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGTTAACTTTGa                                                  >  4:84138/1‑149 (MQ=255)
                                                                                          ttGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGTTAACTTTGaa                                                 <  2:471867/149‑1 (MQ=255)
                                                                                             ccAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGTTAACTTTGAAATc                                              >  2:410135/1‑149 (MQ=255)
                                                                                              cAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGTTAACTTTGAAATCa                                             >  2:316074/1‑149 (MQ=255)
                                                                                                ggAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAac                                                                                                                                  <  4:171133/62‑1 (MQ=255)
                                                                                                ggAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAac                                                                                                                                  >  3:171133/1‑62 (MQ=255)
                                                                                                      aaCTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGTTAACTTTGAAATCAGTTAAAAc                                     <  3:96858/149‑1 (MQ=255)
                                                                                                        cTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGTTAACTTTGAAATCAGTTAAAACGc                                   >  4:71327/1‑149 (MQ=255)
                                                                                                         tCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGTTAACTTTGAAATCAGTTAAAACGCt                                  >  1:585817/1‑149 (MQ=255)
                                                                                                         tCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGTTAACTTTGAAATCAGTTAAAACGCt                                  >  2:561888/1‑149 (MQ=255)
                                                                                                         tCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGTTAACTTTGAAATCAGTTAAAACGCt                                  >  4:134166/1‑149 (MQ=255)
                                                                                                           aTCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCgaga                                                                          <  3:21330/107‑1 (MQ=255)
                                                                                                           aTCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCgaga                                                                          >  4:21330/1‑107 (MQ=255)
                                                                                                               gcggcgACAAGCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGGTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGTTAACTTTGAAATCAGTTAAAACGCTATCGGc                            <  2:479644/149‑1 (MQ=255)
                                                                                                                      cAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGTTAACTTTGAAATCAGTTAAAACGCTATCGGCTACCGGa                     <  4:359851/149‑1 (MQ=255)
                                                                                                                       aaTCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGTTAACttt                                                    >  3:72705/1‑117 (MQ=255)
                                                                                                                       aaTCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGTTAACttt                                                    <  4:72705/117‑1 (MQ=255)
                                                                                                                       aaTCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGTTAACTTTGAAATCAGTTAAAACGCTATCGGCTACCGGAg                    >  4:125383/1‑149 (MQ=255)
                                                                                                                       aaTCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGTTAACTTTGAAATCAGTTAAAACGCTATCGGCTACCGGAg                    >  2:354113/1‑149 (MQ=255)
                                                                                                                       aaTCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGTTAACTTTGAAATCAGTTAAAACGCTATCGGCTACCGGAg                    >  4:397544/1‑149 (MQ=255)
                                                                                                                       aaTCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGTTAACTTTGAAATCAGTTAAAACGCTATCGGCTACCGGAg                    >  3:19525/1‑149 (MQ=255)
                                                                                                                        aTCATGGCGCAACTTTTGGTAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGTTAACTTTGAAATCAGTTAAAACGCTATCGGCTACCGGAGc                   <  3:133051/149‑1 (MQ=255)
                                                                                                                        aTCATGGCGCAACTTTTGGTAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGTTAACTTTGAAATCAGTTAAAACGCTATCGGCTACCGGAGc                   <  3:96342/149‑1 (MQ=255)
                                                                                                                         tCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGTTAACTTTGAAATCAGTTAAAACGCTATCGGCTACCGGAGCg                  >  1:250017/1‑149 (MQ=255)
                                                                                                                          cATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCTCCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGTTAACTTTGAAATCAGTTAAAACGCTATCGGCTACCGGAGCgg                 <  4:57547/149‑1 (MQ=255)
                                                                                                                           aTGGCGCAACTTTTGGAAATCGCGTCTTGCTCAACCACCTGGCGGTAGGAGAATTGTTGATCTACCAGAATCCCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGTTAACGTTGAAATCAGTTAAAACGCTATCGGCTACCGGAGCggg                >  2:185892/1‑149 (MQ=255)
                                                                                                                                gcAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGTTAACTTTGAAATCAGTTAAAACGCTATCGGCTACCGGAGCGGGTGccc           >  3:44363/1‑149 (MQ=255)
                                                                                                                                  aaCTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGTTAACTTTGAAATCAGTTAAAACGCTATCGGCTACCGGAGCGGGTGCCCCa         <  2:403960/149‑1 (MQ=255)
                                                                                                                                     ttttGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGTTAACTTTGAAATCAGTTAAAACGCTATCGGCTACCGGAGCGGGTGCCCCAGcc      >  4:212252/1‑149 (MQ=255)
                                                                                                                                     tttgGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATGTACCAGAATCGCCGAGTCATAAGGGGAGAGGGCCTTTGCAGCGTTAACTTTGAAATCAGTTAAAACGCTATCGGCTACCGGAGCGGGTGCCCCAGcc      >  2:557390/1‑149 (MQ=255)
                                                                                                                                         ggCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGTTAACTTTGAAATCAGTTAAAACGCTATCGGCTACCGGAGCGGGTGCCCCAGCCGCgg  <  2:283529/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                  
ACCAGCAGTTTTAAGGCTTTGCCACCGTCTTGTTTGATCTGTAGCGGATTGATTTTGCGGTCAATAACCACGCTATCGACCGGAATACCGTTGCCAGGAATGAACTCATCGGCGGCGACAATCATGGCGCAACTTTTGGCAATCGCGTTTTGCTCAACCACCTGGCGGTAGCAGAATTGTTGATCTACCAGAATCGCCGAGGCATAAGGCGAGAGGGCCTTTGCAGCGTTAACTTTGAAATCAGTTAAAACGCTATCGGCTACCGGAGCGGGTGCCCCAGCCGCGG  >  NC_000913/4072272‑4072557

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: