Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913535,9010TG45.0% 124.0 / 13.6 80G64G (GGT→GGGhyihydroxypyruvate isomerase
Reads supporting (aligned to +/- strand):  ref base T (13/31);  new base G (36/0);  total (49/31)
Fisher's exact test for biased strand distribution p-value = 4.19e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATGATTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTT  >  NC_000913/535751‑536047
                                                                                                                                                      |                                                                                                                                                  
atatGATTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGt                                                                                                                                                    <  2:287011/151‑1 (MQ=255)
atatGATTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGt                                                                                                                                                    <  2:673943/151‑1 (MQ=255)
 tatGATTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTg                                                                                                                                                   >  2:419823/1‑151 (MQ=255)
  atGATTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGa                                                                                                                                                  <  2:288694/151‑1 (MQ=255)
     aTTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTg                                                                                                                                                   >  2:144030/1‑147 (MQ=255)
       tttCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTggg                                                                                                                                             >  1:277699/1‑151 (MQ=255)
          ctcGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTgggggg                                                                                                                                          >  1:58075/1‑151 (MQ=255)
            cGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGggcggc                                                                                                                                         <  1:411696/150‑1 (MQ=255)
            cGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGgcggcg                                                                                                                                        <  2:27668/151‑1 (MQ=255)
              cccGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGGTGGGGGGCggg                                                                                                                                      >  1:354014/1‑151 (MQ=255)
              cccGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAGGGGGGGGCggg                                                                                                                                      >  1:598595/1‑151 (MQ=255)
               ccGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCgggg                                                                                                                                     <  1:652304/151‑1 (MQ=255)
                    ttGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGGGGGGCggggggg                                                                                                                                  >  1:703835/1‑147 (MQ=255)
                       agaAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGGTGGGGGGGGGGGGGGGGCgg                                                                                                                             >  1:307237/1‑151 (MQ=255)
                         aaaGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGggcg                                                                                                                              >  2:930496/1‑148 (MQ=255)
                         aaaGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTa                                                                                                                           <  1:129151/151‑1 (MQ=255)
                          aaGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGCGGGGGGGGGCgg                                                                                                                             >  1:168880/1‑148 (MQ=255)
                           aGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTAtt                                                                                                                         <  2:785448/151‑1 (MQ=255)
                           aGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTTGGCGGCGGGGGGGGGCGGTGtt                                                                                                                         >  2:551839/1‑151 (MQ=255)
                              tGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGGGGGGGGGCGGGGGTTTGc                                                                                                                       >  2:353300/1‑150 (MQ=255)
                              tGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAGTGGGGGGCGGGGGAGGGCGGTCTTGcc                                                                                                                      >  1:876883/1‑151 (MQ=255)
                               gcgcAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCt                                                                                                                     <  1:657685/151‑1 (MQ=255)
                                        ggTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCt                                                                                                             <  1:848834/150‑1 (MQ=255)
                                           tttCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGGGGCGGGGGGGGGCGGGATTGGCTGTTTTTCTGGcc                                                                                                         >  2:383446/1‑151 (MQ=255)
                                            ttCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGAGCGCGGGATTGCCTGTATTCCttg                                                                                                           >  2:699948/1‑146 (MQ=255)
                                                cggcgTTGAATTTATGTGTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTCCCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTgaa                                                                                                    <  1:551290/151‑1 (MQ=255)
                                                 ggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGCGCGGGATTGCCTGTATTTCTGGGCggg                                                                                                      >  2:673301/1‑146 (MQ=255)
                                                  gcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTgaaga                                                                                                  <  2:396149/151‑1 (MQ=255)
                                                  gcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGGTGGGGGGGGGGGGGGGGGGGGATTGCCTGTCtttcctgc                                                                                                          >  1:622240/1‑137 (MQ=255)
                                                   cgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTgaagaa                                                                                                 <  2:480257/151‑1 (MQ=255)
                                                   cgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGATTGGGGGCGGGGGGGGGCGGGAATGCCTGTATTCCTGGGCg                                                                                                        >  1:64891/1‑144 (MQ=255)
                                                       gAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGt                                                                                             <  1:161066/151‑1 (MQ=255)
                                                            tATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCg                                                                                         <  2:1010239/150‑1 (MQ=255)
                                                            tATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGGGCGGCGGGGGGGCGCGGGATTGCCTTTATTTCTGGCCGTGGAGAAGAGTTTTgg                                                                                        >  2:460658/1‑151 (MQ=255)
                                                               gTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATTCGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGAt                                                                                     <  1:881321/151‑1 (MQ=255)
                                                               gTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGAt                                                                                     <  1:227850/151‑1 (MQ=255)
                                                               gTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGAt                                                                                     <  2:813217/151‑1 (MQ=255)
                                                                tttCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCGCAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGAt                                                                                     <  1:806326/150‑1 (MQ=255)
                                                                tttCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGGGCGGCGGGGGAGCGCGGGATTGCCTGGATTTCTCGCCGTGAAGGAAAGGTTTGGGGTg                                                                                    >  2:10356/1‑151 (MQ=255)
                                                                      tATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGCGGCGGGGGGGCGCGGGATTGCCTGTATTCCTGGCCGGGAAGAAGAGGTTCGGGGTGGCg                                                                                 >  2:317938/1‑148 (MQ=255)
                                                                        tGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGCGGGGGAGCGCGGGATTGCCCGTATTCCTGGCCGTGGAGAAGAGTTTTGGGGTgggggag                                                                              >  2:342691/1‑144 (MQ=255)
                                                                         gACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTagca                                                                            <  2:834771/150‑1 (MQ=255)
                                                                         gACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAGTGGGCGGCGGGGGGGGGCGGGATTTCCTGTATTTCCGGGCGGGGAGAAGAGGTTTggggggg                                                                                   >  1:610368/1‑139 (MQ=255)
                                                                          aCTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGGGCGGGATTTCCTGTATTCCTGGCCGGGGAGAAGAGGTTCGGGGTgggggag                                                                              >  2:93156/1‑142 (MQ=255)
                                                                           cTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGGGGGGGAGCGCGGGATTGCCTGTATTTCTGGCCGTGGAGAAGAGGTTTGGGGTGGCGTagca                                                                            >  1:990398/1‑148 (MQ=255)
                                                                            tACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGGGGGGGAGCGCGGGATTTGCTGGATTTCTGGCCGGGAAGAAAAGGTTCGGGGTgggg                                                                                 >  2:259052/1‑140 (MQ=255)
                                                                             aCGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTagcagca                                                                         <  1:25015/149‑1 (MQ=255)
                                                                             aCGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGGGCGGCGGGGGGGCGCGGGATTGCCTGTATTTCTGGCCGTGAAGAAGAGGTTCGGGGTGGCg                                                                                 >  2:477836/1‑141 (MQ=255)
                                                                              cGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGCGGGGGGGCGCGGGATTGACTGTATTCCTGGCCGGGGAGGAGAGGTTTggggtggggtgggcac                                                                          >  1:143736/1‑141 (MQ=255)
                                                                               gACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGGGGGGGGGCGGGGGATTTCCTGTTTTTCTGGGCGGGGAGAAAAGGTTTGGGGGGGGGTTGGAGCAGCg                                                                      >  1:566830/1‑150 (MQ=255)
                                                                                    tGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGCGGCGGGGGGGGGCGGGATTGGCTGTTTTCCTTGCCGTGGAGAGGAGGTTTGGGGTGGCGTAGGAGCAGCGAtttg                                                                 >  1:279395/1‑148 (MQ=255)
                                                                                       agaaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTAt                                                             <  1:532234/151‑1 (MQ=255)
                                                                                        gaaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGGGGCGGGGGGGCGCGGGATTTCCTGTATTTCTTGCCGTGAAGAAGAGGTTCGGGGTGGGGGAGCAGCAGCGATTCGTTAt                                                             >  2:14742/1‑150 (MQ=255)
                                                                                         aaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATg                                                            <  1:995383/150‑1 (MQ=255)
                                                                                               aaCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGTCTGGGCGGGGGGGGGGCGCGGGATTGCCTGTTTTCCTGGCCGGGAAGAAGAGGTTCGGGGTGGCGGAGCAGCAGCGATTCGTTATGccc                                                         >  2:143967/1‑147 (MQ=255)
                                                                                               aaCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGCGGGGGGGGGCGGGATTGCCTGTATTTCTGGCCGTGGAGAAGAGGTTTGGGGTGGGGGAGGAGCAGCGATTCGTTATGccc                                                         >  2:603105/1‑147 (MQ=255)
                                                                                                aCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGGGGGGGGGGGCGGGATTGCCTTTATTTCCGGCCGGGGAGGAGAGTTTCGGGGTGGGGGGGCAGCAGGGGTTTGtttttgcccgg                                                      >  1:491570/1‑140 (MQ=255)
                                                                                                   tgtgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGCGGGGGGGCGCGGGATTGCCTGTATTTCCGGGCGTGGAGAAGAGGTTTGGGGGGGGGGGGGAGCAGCGATTCGTTATGCCCGTGGGCt                                                  >  1:640198/1‑150 (MQ=255)
                                                                                                      tcTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGGGCGGGAATGCCTGTATTTCTGGCCGCGGAGAAGAGGTTTGGGGTGGGGGAGGCGCAGCGAGTTGTTATGGCCCCGCGGTggg                                               >  1:340520/2‑150 (MQ=255)
                                                                                                      gCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTggg                                               <  2:395377/150‑1 (MQ=255)
                                                                                                       cTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGGGGGGGGGGGCGGGGTTGGCTTTTTTTCTGGCCGGGGAGGAGGGTTTCGGGGGGGCGCTGCAGGAGGGGATTGTTTTGCCCCTGCGCCggg                                               >  1:87230/1‑149 (MQ=255)
                                                                                                            gAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATaaa                                        <  2:206436/151‑1 (MQ=255)
                                                                                                             aGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGCGGGGGGGCGCGGGATTGCCTGTATTTCTGGCCGGGGAGAAGAGGTTTGGGGTGGGGGAGGAGCAGCGGTTTGTTTTGCCCGTGCGCTGGGGAATaaa                                        >  1:381594/1‑150 (MQ=255)
                                                                                                              gTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGATGGGGGGCGGGGGGGGGCGGGGTTTCCTGTATTTCTGGCCGGGGAGGAGAGGTTTGGGGGGGGGGAGGCGGAGCGGTTCGTTTTGTCCGGGCGCTGGGGGATaaaa                                       >  1:598525/1‑150 (MQ=255)
                                                                                                                 ataaACTCGAACATACGCTGCACAATTTACCGGCGGGTGACAGGGCGGCGGGGGGGGGGCGGGTTGGCTGTTCTCCTGGCCGGGGTGAAGAGGGTAGGGGTGGCGTAGGCGGGGCGATTCGTTTTGGCCGGGGGCTGGGGATTaaaataa                                    >  1:55557/1‑150 (MQ=255)
                                                                                                                   aaaCTCGAACATACGCTGCACAATTTACCGGCGGGTGGCTGGGGGGCGGGGGAGAGCGGGAATTCCCGTATTCCTGGCCGCGAAGGAGAGGTTTGGGATGGCGGAGCAGCAGCGATTCGTTATGCCCGCGCGCTGGGGAATaaaaaaa                                    >  2:654630/1‑148 (MQ=255)
                                                                                                                    aaCTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTaa                                <  2:168895/151‑1 (MQ=255)
                                                                                                                     aCTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGCGGGGGGGCGCGGGGTTTCCTTTATTCCTTGCCGGGGAGAAGAGGTTTGGGGGGGCGTAGCAGCAACGATTTGGTTTGCCCGCGCGCTGGGGAATaaaaaaa                                    >  1:561353/1‑146 (MQ=255)
                                                                                                                      cTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGAGGGCGGGATTGACTTTATTCCTGGGCGTGGAGGAGAGGTTTGGGGTGGGGGAGCAGCAGGGATTCTGTATTGCCGTGCGCTGGGGGATAAAAAAAATAACt                              >  1:632036/1‑151 (MQ=255)
                                                                                                                        cGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGt                            <  2:640251/151‑1 (MQ=255)
                                                                                                                        cGAACATACGCTGCACAATTTACCGGCGGGTGACTAGGCGGCGGGGGGGGGCGGGTTTAGCTGTATTCCTGGCCGTGCAGAAGGAGTTAGGGGTAGCGGAGCAGCAGCGATTCGTTATGACCGTGTGCTGGGGAATAAAAAAAATAACt                              >  2:586608/1‑149 (MQ=255)
                                                                                                                             aTACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGAGCGCGGGATTGGCTGTATTTCTGGCCGTGAAGAAGAGGTTTGGGGTGGCGTAGCAGCAGCGATTCGTTATGCCCGGGCGCTGGGGAATAAAAAAAATAAATGTCTgg                        >  2:945629/1‑150 (MQ=255)
                                                                                                                                cGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCgg                    <  2:16098/151‑1 (MQ=255)
                                                                                                                                 gCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGt                   <  2:703891/151‑1 (MQ=255)
                                                                                                                                  cTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGGGCGGGATTGCCTTTATTTCTGGCCGTGGAGAAGAGGTTCGGGGTGGGGGAGCAGGAGCGATTCGTTATGCCCGGGCGCGGGGGGATAAAAAAAATaaaatgctgggtg                     >  1:911016/1‑137 (MQ=255)
                                                                                                                                         aTTTACCGGCGGGGGACTGGGGGGGGGGGGGGGGCGGGATTTCCTGTATTTCTGGCCGGGAAGAAGAGGTTCGGGGTGGGGGAGCAGGAGCGGTTTGTTTTGCCCGTGCGCTGGGTAATAAAAAAAATAAATGTCTGGTCGGGAAAAc              >  1:1006974/1‑148 (MQ=255)
                                                                                                                                           ttACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCg          <  2:43591/150‑1 (MQ=255)
                                                                                                                                             aCCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCt       <  2:598645/151‑1 (MQ=255)
                                                                                                                                              ccGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTg      >  2:735927/1‑151 (MQ=255)
                                                                                                                                              ccGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTg      <  2:881059/151‑1 (MQ=255)
                                                                                                                                                ggcggGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGt    <  2:370512/151‑1 (MQ=255)
                                                                                                                                                  cggGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGttt  >  1:673166/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
ATATGATTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTT  >  NC_000913/535751‑536047

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: