Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,070,413 T→G 6.2% intergenic (‑12/+102) yihR ← / ← yihS putative sulphoquinovose mutarotase/sulphoquinovose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,070,4130TG 6.2% 356.3 / 10.0 107intergenic (‑12/+102)yihR/yihSputative sulphoquinovose mutarotase/sulphoquinovose isomerase
Reads supporting (aligned to +/- strand):  ref base T (46/54);  new base G (4/3);  total (50/57)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

Only 100 of 108 total aligned reads displayed.

CCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGC  >  NC_000913/4070266‑4070559
                                                                                                                                                   |                                                                                                                                                   
ccTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGt                                                                                                                                                     >  2:807974/1‑148 (MQ=255)
 cTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTg                                                                                                                                                    <  2:90353/148‑1 (MQ=255)
 cTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGa                                                                                                                                                   <  3:476471/149‑1 (MQ=255)
   gAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAg                                                                                                                                                  >  2:707175/1‑148 (MQ=255)
        gTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑a                                                                                                                                            <  1:778449/148‑1 (MQ=255)
        gTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑a                                                                                                                                            <  2:579516/148‑1 (MQ=255)
        gTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑a                                                                                                                                            <  1:694574/148‑1 (MQ=255)
        gTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑a                                                                                                                                            <  2:691468/148‑1 (MQ=255)
             ttCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCaa                                                                                                                                      >  3:862942/1‑149 (MQ=255)
              tCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCaa                                                                                                                                      >  2:88575/1‑148 (MQ=255)
                     gACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCaaa                                                                                                                               >  2:61615/1‑148 (MQ=255)
                        cAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACt                                                                                                                           >  3:851699/1‑149 (MQ=255)
                         aGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACt                                                                                                                           >  1:78269/1‑148 (MQ=255)
                           cACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTTCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCa                                                                                                                        <  3:321406/149‑1 (MQ=255)
                              cGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACAT‑CAGCCTCACGGGAGGCTgta                                                                                                                                            >  4:339470/1‑124 (MQ=255)
                              cGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACAT‑CAGCCTCACGGGAGGCTgta                                                                                                                                            <  3:339470/126‑3 (MQ=255)
                                aCGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGa                                                                                                                   >  4:293281/1‑149 (MQ=255)
                                    tCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCaaaaa                                                                                                                             >  2:547057/1‑135 (MQ=255)
                                    tCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCaaaaa                                                                                                                             <  1:547057/135‑1 (MQ=255)
                                     cacaATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCaa                                                                                                                                      >  2:450635/1‑125 (MQ=255)
                                     cacaATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCaa                                                                                                                                      <  1:450635/125‑1 (MQ=255)
                                     cacaATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTc                                                                                                              <  3:105743/149‑1 (MQ=255)
                                       caATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTccc                                                                                                            >  3:23091/1‑149 (MQ=255)
                                         aTGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTcccc                                                                                                           >  2:14978/1‑148 (MQ=255)
                                         aTGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTcccc                                                                                                           >  2:664190/1‑148 (MQ=255)
                                         aTGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTcccc                                                                                                           >  2:5637/1‑148 (MQ=255)
                                         aTGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTcccc                                                                                                           >  1:667840/1‑148 (MQ=255)
                                         aTGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTcccc                                                                                                           >  1:314549/1‑148 (MQ=255)
                                           ggTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCaaaaa                                                                                                                             >  2:397632/1‑128 (MQ=255)
                                           ggTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCaaaaa                                                                                                                             <  1:397632/128‑1 (MQ=255)
                                            gTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCg                                                                                                        <  2:78269/148‑1 (MQ=255)
                                                  tGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCttt                                                                                                 <  3:170931/149‑1 (MQ=255)
                                                  tGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCtt                                                                                                  <  1:664190/148‑1 (MQ=255)
                                                   gATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                <  3:2833/149‑1 (MQ=255)
                                                    aTAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                <  1:880690/148‑1 (MQ=255)
                                                    aTAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                <  2:829769/148‑1 (MQ=255)
                                                    aTAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                <  2:171282/148‑1 (MQ=255)
                                                    aTAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                <  2:472463/148‑1 (MQ=255)
                                                     tAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                                                                                              <  4:514997/149‑1 (MQ=255)
                                                     tAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                                                                                              <  4:23091/149‑1 (MQ=255)
                                                      aaTCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                                                                                              <  1:807974/148‑1 (MQ=255)
                                                           cccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGagag                                                                                         <  1:517821/148‑1 (MQ=255)
                                                           cccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAgg                                                                                        <  4:444531/149‑1 (MQ=255)
                                                            ccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCaaaa                                                                                                                              >  2:830989/1‑110 (MQ=255)
                                                            ccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCaaaa                                                                                                                              <  1:830989/110‑1 (MQ=255)
                                                            ccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAgg                                                                                        >  2:521104/1‑148 (MQ=255)
                                                                 gCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                                                                                  <  4:774799/149‑1 (MQ=255)
                                                                 gCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAg                                                                                   <  1:14978/148‑1 (MQ=255)
                                                                 gCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAg                                                                                   <  1:521647/148‑1 (MQ=255)
                                                                   aaGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg                                                                                 >  2:531093/1‑148 (MQ=255)
                                                                      tGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGgaga                                                                              <  2:710476/148‑1 (MQ=255)
                                                                      tGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGagag                                                                             >  4:743159/1‑149 (MQ=255)
                                                                       gAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGagag                                                                             <  2:667840/148‑1 (MQ=255)
                                                                       gAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgg                                                                            <  4:21308/149‑1 (MQ=255)
                                                                          tgggTCTGTCCTCTGCAATGCATTTTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACGGCCTCCCGGGAGGCTG‑ATGTCAAACCCAAAAACGCCCCCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGa                                                                         <  4:851699/147‑1 (MQ=38)
                                                                                 gTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACTACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgg                                                                            >  1:486259/1‑139 (MQ=255)
                                                                                 gTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAGCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgg                                                                            <  2:486259/139‑1 (MQ=255)
                                                                                  tCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑At                                                                                                                                           >  4:52140/1‑75 (MQ=255)
                                                                                  tCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑At                                                                                                                                           <  3:52140/75‑1 (MQ=255)
                                                                                  tCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGc                                                                 >  4:456074/1‑149 (MQ=255)
                                                                                    cACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAc                                                               >  3:619562/1‑149 (MQ=255)
                                                                                     aCTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAcc                                                              >  4:776563/1‑149 (MQ=255)
                                                                                     aCTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAc                                                               >  2:257913/1‑148 (MQ=255)
                                                                                      cTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaaaa                                                                      <  1:420135/140‑1 (MQ=255)
                                                                                            tGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGtt                                                                                     >  4:116351/1‑119 (MQ=255)
                                                                                            tGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGtt                                                                                     <  3:116351/119‑1 (MQ=255)
                                                                                             gCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCCCCGGACGGTCCCCTCGCCCCTTGGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt                                                       <  1:201806/148‑1 (MQ=255)
                                                                                             gCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt                                                       <  1:257913/148‑1 (MQ=255)
                                                                                             gCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt                                                       >  2:691085/1‑148 (MQ=255)
                                                                                              cATATTGGTAATTTGCATACAAGGCCCTTTAATAAACGACATACAGCCCCCCGTGAGGCTG‑TGGTCAACCCCAAAAACTGCACCGGCCGGCCCCCCCGCCCCTTGGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTc                                                      <  1:655298/148‑1 (MQ=255)
                                                                                                  ttGGTAATTTGCAAACAGGGCACCTTAATTAACAACATACCGCCTCCCGTGTGGCTG‑ATGTCAAACCCAAAACCGCCACCGGCCGCGCCCCTCGCCCCTTTGGGGGGAGGGTAAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGAt                                                  <  1:707175/148‑1 (MQ=255)
                                                                                                   tGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAgg                                                                                                                                                 >  2:54138/1‑53 (MQ=255)
                                                                                                   tGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAgg                                                                                                                                                 <  1:54138/53‑1 (MQ=255)
                                                                                                   tGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                <  2:717730/101‑1 (MQ=255)
                                                                                                   tGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                                                                                >  1:717730/1‑101 (MQ=255)
                                                                                                   tGGGGATTTGCATACAAGGCACCTTAATAACCAACATACCCCCCCCGTTGAGGCTG‑TGGACAACCCCAAAAGGCGCACCGGGCGATCCCCTCCCCTTTTTGGAAAGGGGGTAAGGGAGAGGGAAAAACCCGCACCCTATTTTCAGATc                                                 <  1:848598/148‑1 (MQ=255)
                                                                                                     gTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGt                                              <  3:220150/149‑1 (MQ=255)
                                                                                                              aTACAGGGCTCCTTAATTAACGACATACGGCCCCCCGGGGGGCGG‑ATGTAAACCCCAAAAACTGCGCCGGACGGCCCCCTCGCCCTTTGGGGGAGAGGGTTAGGGAGGGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCa                                     <  3:456074/149‑1 (MQ=255)
                                                                                                                aCAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGATAGGGAGAGGGGaaa                                                                       >  1:100200/1‑113 (MQ=255)
                                                                                                                aCAAGGCACCTAAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTGGGGGAGAGGGTTAGGGAGAGGGGaaa                                                                       <  2:100200/113‑1 (MQ=255)
                                                                                                                  aaGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTaa                                 >  3:584451/1‑149 (MQ=255)
                                                                                                                     gCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAata                               >  2:754645/1‑148 (MQ=255)
                                                                                                                        ctttaattaaCGACATCCCGCCTCCCGTGAGGCTG‑ATGAAAAACCCAAAACCTGCACCGGACGGCCCCCTCGCCCTTTGGGGGAGGGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATcc                            <  1:521104/146‑1 (MQ=255)
                                                                                                                        ccTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaaaa                                                                      <  2:455185/106‑1 (MQ=255)
                                                                                                                        ccTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaaaa                                                                      >  1:455185/1‑106 (MQ=255)
                                                                                                                        ccTTAATTAACGACATACAGCCCCCCGTGGGGCTG‑ATGTCAAACCAAAAAACGGCCCCGGCCGGCCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATcc                            <  1:5637/148‑1 (MQ=38)
                                                                                                                          ttaattaaCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCcag                          >  1:292600/1‑148 (MQ=255)
                                                                                                                              ttaaCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGt                                              >  4:227995/1‑124 (MQ=255)
                                                                                                                              ttaaCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGt                                              <  3:227995/124‑1 (MQ=255)
                                                                                                                                 aCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAGCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaaaa                                                                      <  1:16641/97‑1 (MQ=255)
                                                                                                                                 aCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaaaa                                                                      >  2:16641/1‑97 (MQ=255)
                                                                                                                                  ttgacatcAGCCTCACGGGAGGCTGTATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGAccc                   >  1:196418/8‑148 (MQ=38)
                                                                                                                                     cATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCtcc              <  4:619562/149‑3 (MQ=255)
                                                                                                                                        tcAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCaa            >  2:323043/2‑148 (MQ=255)
                                                                                                                                        tcAGCCTCACGGGAGGCTGTATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCa             >  2:75373/2‑148 (MQ=255)
                                                                                                                                        tcAGCCTCACGGGAGGCTGTATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGAGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCa             >  2:477547/2‑148 (MQ=39)
                                                                                                                                            ccTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCt       >  3:846472/1‑149 (MQ=255)
                                                                                                                                              tCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTgg     >  3:150334/1‑149 (MQ=255)
                                                                                                                                              tCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTgg     >  3:712687/1‑149 (MQ=255)
                                                                                                                                                 cGTGAGGCTG‑ATGTAAAACCCAAAAACGGCACCGGACGGCCCCCTCGCCCCTTTGGGGAGAGGTTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGc  <  4:712687/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGC  >  NC_000913/4070266‑4070559

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