Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913535,9180AG43.5% 82.4 / 13.3 63E70G (GAG→GGG) hyihydroxypyruvate isomerase
Reads supporting (aligned to +/- strand):  ref base A (6/29);  new base G (27/0);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 4.39e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCAC  >  NC_000913/535769‑536066
                                                                                                                                                     |                                                                                                                                                        
ttttGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTAGGAGGGGGGGGa                                                                                                                                                          >  2:913297/1‑150 (MQ=255)
  ttGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAgc                                                                                                                                                        <  1:98405/150‑1 (MQ=255)
  ttGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAgc                                                                                                                                                        <  1:526090/150‑1 (MQ=255)
              cgcAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTg                                                                                                                                         <  2:269849/151‑1 (MQ=255)
              cgcAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTg                                                                                                                                         <  2:13493/151‑1 (MQ=255)
               gcAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGGCTGGGCGGCCGGGGGGCGCggcagt                                                                                                                                              >  2:756866/1‑141 (MQ=255)
                    gtgGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTcc                                                                                                                                   <  1:588476/151‑1 (MQ=255)
                    gtgGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACCCGGCGGCGGGGGGGCGCGGGAATGCCTGTATTcc                                                                                                                                   >  2:880658/1‑151 (MQ=255)
                      ggTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTg                                                                                                                                 <  1:71659/151‑1 (MQ=255)
                        ttttCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGCGCGGGATTGCCTGTATTTCttg                                                                                                                                >  2:916470/1‑148 (MQ=255)
                         tttCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGc                                                                                                                               <  1:1146965/150‑1 (MQ=255)
                           tCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCg                                                                                                                             <  1:490053/150‑1 (MQ=255)
                              cggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTgaa                                                                                                                         <  1:917555/151‑1 (MQ=255)
                                  gTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGATGGGCGGGGGGGGGGCGCGGGGTTTTCTGTATTTCttg                                                                                                                                >  1:23021/1‑138 (MQ=255)
                                       aTTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGTGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGttt                                                                                                               <  1:823267/151‑1 (MQ=255)
                                        tttATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGCGCGGGATTTGCTGTATTTCTTGGCGTGAAGAAGAGgtt                                                                                                               >  2:853166/1‑150 (MQ=255)
                                              tttCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGAATGGGCGGGGGGGGGGGGGGGGCTTGCCCGTATTCCTGGCCGTGAAGAAGAGCTTTGGGGTg                                                                                                        >  2:763797/1‑151 (MQ=255)
                                               ttCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGCGCGGGATTGCCTGTATTTCTGGGCGTGAAGAAGAGGTTTGGGGTgg                                                                                                       >  2:940179/1‑151 (MQ=255)
                                               ttCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATg                                                                                                        <  1:1001043/150‑1 (MQ=255)
                                                 ccTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGc                                                                                                      <  1:1099891/150‑1 (MQ=255)
                                                 ccTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGGGCGGGGCTGCCTGTATTTCTGGCCGGGGAGGAGAGGTTCgggg                                                                                                          >  2:587034/1‑145 (MQ=255)
                                                  cTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGCGCGGGATTGCCTGTATTTCTGGCCGTGGAGAAGAGGTTTGGGGTGGCg                                                                                                     >  2:413603/1‑150 (MQ=255)
                                                      tGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGGCTGGGGGGCGGGGGAGCGCGGGATTGACTGTATTTCTGGCCGTGAAGAAGAGTTTCggggtggggt                                                                                                    >  1:147881/1‑147 (MQ=255)
                                                       gACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTagcag                                                                                               <  1:910543/151‑1 (MQ=255)
                                                       gACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGCGCGGGATTGGCTGTTTTCCTGGCCGTGAAGAAGAGGTTTGGGGTgg                                                                                                       >  2:670861/1‑143 (MQ=255)
                                                         cTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGCGGCGGGGGGGCGCGGTATTGCCCGGATTTCCGGCCGGGGAGAAAAGTTTCGGGGTGGCGGag                                                                                                  >  2:1194067/1‑146 (MQ=255)
                                                          tACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGCGCGGGATTTCCTGTATTCCTTGGCGGGAAGAAGAGTTTTGGGGTGGGGTagca                                                                                                >  2:375142/1‑147 (MQ=255)
                                                           aCGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGCGCGGGATTGCCTGTATTTCTGGGCGGGAAGAAGGGGTTCGGGGTGGCGTagcagca                                                                                             >  2:27850/1‑149 (MQ=255)
                                                            cGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGTGGAGGGCGGCGGGGGAGCGCGGGAATTACTGTATTCGTGCGCGTGGAGAAAGAGTTCCGGGAggggtagacaca                                                                                              >  1:1128356/1‑139 (MQ=255)
                                                              aCATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGAt                                                                                        <  2:845928/151‑1 (MQ=255)
                                                                aTTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGGGGGGGGGGGCGGGATTGCCTGTATTTCTGGCCGTGAAGAAGAGGTTTGGGGTGGGGTAGCAGCAGCGGtt                                                                                       >  2:705573/1‑150 (MQ=255)
                                                                       aaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGc                                                                               <  2:800924/151‑1 (MQ=255)
                                                                        aTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGCGCGGGATTGCCTCTTTTCCCTGCCGTGAAGAAGAGGTTTGGGAGGGCGGAGCAGCAGCGAGTAGTTATGcc                                                                              >  2:173629/1‑151 (MQ=255)
                                                                        aTTAAAACATGTGCTGGCGAGGAATAAACTCGAACATACGCTGCACAACTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGGGGGTATTCCCTGATTCCCTTGCCGGATGAAGAGGTTCGGGGGGGCGCtgcagcccagatcagttatg                                                                                >  1:706946/1‑135 (MQ=255)
                                                                          tAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGGGGGGGGGGGGGGGATTACTTGTTTTCCCGGCCGTGAAGAAGAGGTTCGGGGGGGGGGAGCAGCGGGGATTTGGTCTGccc                                                                             >  1:614383/1‑150 (MQ=255)
                                                                            agaCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATGGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTAGTTATGCCCg                                                                            <  1:118794/147‑1 (MQ=255)
                                                                              aCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTgcg                                                                       <  2:998897/151‑1 (MQ=255)
                                                                                  gtgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTgg                                                                   <  2:27720/151‑1 (MQ=255)
                                                                                  gtgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTgg                                                                   <  1:1124332/151‑1 (MQ=255)
                                                                                  gtgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTgg                                                                   <  2:600713/151‑1 (MQ=255)
                                                                                     cTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGGGGGGGGGCGGGGGATTGTCTGTATTCCTGGGCGGGGAGAAGAGGTTGGGGGGGGGGTGGGAGGAGGGATTTGGTATGCCCGTGCGCTggggg                                                                >  1:917111/1‑149 (MQ=255)
                                                                                      tGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGCGGGGGGGGGCGGGATTTCCTGTTTTCCTGGGCGGGAAGAAGAGGTTTGGGGTGGGGGAGCAGCAGCGGTTCGTTTTGCCCGGGCGCTggg                                                                  >  2:225924/1‑148 (MQ=255)
                                                                                        gCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGAATGGGCGGGGGGGGGGAGCGGGAATTCCCTTTTTTCCCGCCGTGGAGAAGAGTTTTGGGGGGGGGTAGGAGCAGCGATTCGTTATTCCCCGGGGCTgggggat                                                              >  1:1193177/1‑146 (MQ=255)
                                                                                          gAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGAATGGGCGGCGGGGGGGCGCGGGATTGCCTGTATTCCTGGGCGGGGAGAAGAGGTTCGGGATGGCGGAGCAGGAGCGATTCGTTATGCCCGTGGGCTGGGgaataa                                                            >  1:559100/1‑150 (MQ=255)
                                                                                            gTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGTGCGGCGGGGGAGCGCGGTATCGCCTGTATTCCTGGCCGTGAAGAAGATTTTCGGGATGGCTTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATaaaaa                                                         <  2:600008/151‑1 (MQ=255)
                                                                                             taataaACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGGTTTGGGGTGGCGGAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATaaaaa                                                         >  2:422960/1‑150 (MQ=255)
                                                                                               ataaACTCGAACATACGCTGCACAATTTACCGGCGGGTGGACGCGCGGCGGGGGGGGGCCGGTTTTCCTGCATTTCTCGCCGGGTAGAAGAGGTTCGGGGTGGCGCGGCATCAGCGATTCGCTATGCCCACGAGCCGGGTAATaaaaaa                                                        >  2:299382/1‑149 (MQ=255)
                                                                                               ataaACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTCTGCCCGTGCGCTGGGTAATaaaaa                                                         <  1:1175381/148‑1 (MQ=255)
                                                                                               ataaACTCGAACATACGCTGCACAATTTACCGGCGGGTGAATGGGCGGCGGGGGGGGGCGGGATTGCCCGTATTCCCGGCCGTGAAGAAGAGGTTCTGGATGGCGTAGCAGCCGCGATTTGTTATGCCCGTGCGCTTGGGAATaaaaaaa                                                       >  2:216137/1‑150 (MQ=255)
                                                                                                  aaCTCGAACATACGCTGCACAATTTACCGGCGGGTGACCGGGCGGGGGGGGAGCGCGGGATTTTCTTTATTCCTGGCCGTGAAGAGGGGGTTTGGGGGGGCGTAGCAGCAGCGATTTGTTATGCCCCTGCGCTGGGGAATaaaaaaa                                                       >  1:478563/1‑147 (MQ=255)
                                                                                                    cTCGAACATACGCTGCACAATTTACCGGCGGGTGGCTTGGCGGCGGGGGAGCGCGGGATTGCTTTTTTTCCCGGCCGTGGAGAAGAGGTTTGGGGAGGGGGAGGAGCAGCGATTTGGTTTGGCCGTGGGGTGGGGGAGCAAAAAAATaa                                                   >  1:473146/1‑149 (MQ=255)
                                                                                                    cTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGAATCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAAc                                                  <  1:521278/150‑1 (MQ=255)
                                                                                                    cTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGCGCGGGATTGCCTGTATTTCTGGCCGTGAAGAAGAGGTTTGGGGTGGCGTAGCAGCAGCGATTTGTTATGCCCGTGGGCTGGGGAATAAAAAAAATaa                                                   >  2:881273/1‑149 (MQ=255)
                                                                                                          cATACGCTGCACAATTTACCGGCGGGGGACTAGGCGGGGGGGGGGGGCGGGATTTTCTTTGTTTCTCGCCCGGGTGAAGCGGTTTTGGATGGCGGAGCCCGAGGGATTTGTTTTTACCTTGCGGCGGGGTATAAAAAAAATTaa                                                   >  2:270848/1‑144 (MQ=255)
                                                                                                           aTACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGt                                          <  2:801235/151‑1 (MQ=255)
                                                                                                           aTACGCTGCACAATTTACCGGCGGGGGGCTTGGGGGGGGGGGGGGGCGGGAATGCTTCTATTCCCGGGGCTGGAGGAGAAGTTTGGGGTGGGGGAGCAGGAGCGATTTGGTTTTGCCGGGCGCTGGGGGATAAAAAAAATaaa                                                  >  1:1013812/1‑142 (MQ=255)
                                                                                                           aTACGATGGATAATTTACCGGCGGGTGACTGGGGGGGGTGTGATCGCGGTTTTGCCTGTATTCCTGGCGGGAAAGAAGAGTTTCGGGATGGCGTAGCAGCATCAATTCTTTATGCCAATGCCCGTGGTAATAAAAAAATTTCCTGTCTgtt                                          <  1:816833/151‑3 (MQ=255)
                                                                                                                  gCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTaa                                    <  1:664323/150‑1 (MQ=255)
                                                                                                                  gCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTAGGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTaaa                                   <  1:291812/151‑1 (MQ=255)
                                                                                                                       aTTTACCGGCGGGTGATTGGGAGGCGGGGGGGCGCGGCATAATCTGTATTCCTGGCCGGGAAGAAGGGGTTCGGGTTGGCGTAGCAGCAGCGATCTGTTTTGACCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTTGGGAAAACGcc                              >  1:1092307/1‑151 (MQ=255)
                                                                                                                              ggcggGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGt                       >  2:989029/1‑151 (MQ=255)
                                                                                                                                   gTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTcag                  <  1:875208/151‑1 (MQ=255)
                                                                                                                                           gcggcgGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAAc          <  1:565393/151‑1 (MQ=255)
                                                                                                                                                   ggAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCAc  <  2:440590/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                        
TTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCAC  >  NC_000913/535769‑536066

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