Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913535,9010TG40.0% 179.8 / 15.0 96G64G (GGT→GGGhyihydroxypyruvate isomerase
Reads supporting (aligned to +/- strand):  ref base T (15/42);  new base G (38/0);  total (54/42)
Fisher's exact test for biased strand distribution p-value = 1.41e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGC  >  NC_000913/535755‑536051
                                                                                                                                                  |                                                                                                                                                       
gATTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACt                                                                                                                                                    >  4:373700/1‑151 (MQ=255)
gATTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACt                                                                                                                                                    >  1:1173526/1‑151 (MQ=255)
 aTTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTg                                                                                                                                                   <  4:446171/151‑1 (MQ=255)
 aTTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTg                                                                                                                                                   <  1:1326147/151‑1 (MQ=255)
   tttCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTgg                                                                                                                                                  <  1:660457/150‑1 (MQ=255)
        cGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTgggggggg                                                                                                                                            >  1:488094/1‑149 (MQ=255)
         gCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGgggggggg                                                                                                                                           >  4:197428/1‑151 (MQ=255)
           ccGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCgggg                                                                                                                                         >  2:538027/1‑151 (MQ=255)
           ccGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGATGGGGGGCgggg                                                                                                                                         >  1:840011/1‑151 (MQ=255)
           ccGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGATGGGCGGCgggg                                                                                                                                         >  3:1101439/1‑151 (MQ=255)
           ccGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCgggg                                                                                                                                         >  3:1335290/1‑151 (MQ=255)
            cGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTgggggggggggg                                                                                                                                        >  3:958347/1‑151 (MQ=255)
            cGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCggggg                                                                                                                                        >  4:573687/1‑151 (MQ=255)
              ttttGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGCggggggg                                                                                                                                      >  1:1234288/1‑149 (MQ=255)
               tttGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTgggggggggggggg                                                                                                                                      >  3:1113255/1‑148 (MQ=255)
                ttGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAgc                                                                                                                                     <  2:177892/150‑1 (MQ=255)
                ttGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGATGGGGCGCggggg                                                                                                                                        >  2:443180/1‑147 (MQ=255)
                ttGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGAggg                                                                                                                                    >  3:59478/1‑149 (MQ=255)
                  gagaAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGcgcg                                                                                                                                  <  2:59638/151‑1 (MQ=255)
                   agaAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGCGGGGGGGGGCgg                                                                                                                                 >  1:405903/1‑151 (MQ=255)
                     aaaGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTTGGCGGGGGGGGGGCGCgg                                                                                                                                 >  3:363779/1‑149 (MQ=255)
                        gCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTg                                                                                                                            <  1:662904/151‑1 (MQ=255)
                          tGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGGGGGGGGGGGGCGGGGGATTGcc                                                                                                                          >  1:1344320/1‑151 (MQ=255)
                          tGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGATCTGGGGGCGGGGGGGGGCgggggggac                                                                                                                          >  1:269489/1‑144 (MQ=255)
                            cgcAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCGGGGCGGCGGGGGGGCGCgggagt                                                                                                                             >  1:1370185/1‑142 (MQ=255)
                             gcAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGt                                                                                                                       <  1:148195/151‑1 (MQ=255)
                              cAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTa                                                                                                                      <  3:1230825/151‑1 (MQ=255)
                                    ggTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCt                                                                                                                 <  1:1101307/150‑1 (MQ=255)
                                    ggTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTg                                                                                                                <  2:552706/151‑1 (MQ=255)
                                         tCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTTACCTGGCGGGCGGGGGGGGCGGTATTGCCGGGAGTCCTg                                                                                                                >  3:216015/1‑146 (MQ=255)
                                            cggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTGCCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTgaa                                                                                                        <  3:195910/151‑1 (MQ=255)
                                            cggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGCGGGGGGGGGCGGGGTTTGCTGTTTTcc                                                                                                                  >  1:907389/1‑141 (MQ=255)
                                            cggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGCGCGGGATTGCCTGTATTTCTGGcc                                                                                                             >  2:512207/1‑146 (MQ=255)
                                             ggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGCGGCGGGGGGGGGCGGGATTTGCTGTTTTTCttggcg                                                                                                            >  4:539140/1‑141 (MQ=255)
                                             ggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGCGGGGGATTGGCTGTATTTCTGGCCGTg                                                                                                          >  1:262404/1‑148 (MQ=255)
                                                gTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGaagaag                                                                                                    <  2:70933/151‑1 (MQ=255)
                                                   gAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGt                                                                                                 <  2:288968/151‑1 (MQ=255)
                                                   gAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGGGGGGGAGCGCGGGATTGCCTGTTTTTCTGGGCggggagga                                                                                                     >  1:325832/1‑140 (MQ=255)
                                                     aTTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGttt                                                                                               <  3:625585/151‑1 (MQ=255)
                                                      tttATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTCGGTGGGGGGGGGGCGCGGGGATGCCTCTTTTTCCGGCCGTGGaggaga                                                                                                   >  1:1356706/1‑146 (MQ=255)
                                                       ttATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTc                                                                                              <  3:1009353/150‑1 (MQ=255)
                                                       ttATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGCGCGGGATTGCCTGTATTTCTGGCCGGGAAGAAGAGgtt                                                                                               >  3:481293/1‑149 (MQ=255)
                                                         aTGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCgg                                                                                            <  3:263760/150‑1 (MQ=255)
                                                          tGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGCGGGGGAGGGCGGGATTGGCTGTATTCCTGGCCGGGGAGAAGAGGTTTgggg                                                                                          >  1:1027407/1‑150 (MQ=255)
                                                                cTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGCGGGGGAGGGCGGGATTTCCTGTATTTCTGGCCGGGAAGAAGAGTTTTGGGGTgg                                                                                       >  4:194308/1‑148 (MQ=255)
                                                                 ttATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGCGGGGGGGGGGGGGATTGGCTGTATTTCTGGCCGGGGAGAAGAGGTTTGGGGTgg                                                                                       >  4:300280/1‑147 (MQ=255)
                                                                    tGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGGGGCGGGGGGGCGCGGGATTGCCTGTATTTCTGGCCGTGGAGAAGAGGTTCGGGGTgggg                                                                                     >  2:151583/1‑144 (MQ=255)
                                                                     gACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTagca                                                                                <  3:741576/150‑1 (MQ=255)
                                                                      aCTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGGGGGGGGGGGGGGCGGGATTGCCTGTATTTCTGGCCGTGAAGAAGAGGTTCGGGGTGGGGTagca                                                                                >  3:47380/1‑149 (MQ=255)
                                                                       cTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGGGGCGGGGGGGGGCGGGATTGGCTTTATTTCTTGGCGGGAAGAAGAGGTTTgggggggggtgggag                                                                               >  4:253697/1‑137 (MQ=255)
                                                                         aCGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGGGGGGGGGGGGGGGATTTCCTGTATTCCTGGGCGTGGAGAAAAGGTTTgggggg                                                                                        >  1:1319395/1‑135 (MQ=255)
                                                                          cGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGAGAATGGGGGGCGGGGGGGGGCGGTGTTGCCTGTATTCCTGGCCGTGGAGGAGAGTTTCgggggggtggagca                                                                                >  4:183900/1‑134 (MQ=255)
                                                                            aCATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGGGAAGAAGAGTTTCGGGATGGCGTAGCAGCATCGAt                                                                        <  4:465581/151‑1 (MQ=255)
                                                                              aTTGAAGAATTAAAACATGTGCTGGCGAGTGATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGCGGGGGGGCGCGGGATTGCCTGTATTCCTGGCCCTGGAGAAGAGGTTTGGGGGGGGGGAGCAGCAGCGAttt                                                                      >  1:927846/1‑150 (MQ=255)
                                                                              aTTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGCGCGGGATTGCCTGTATTTCTGGCCGTGGAGAAGAGGTTCGGGGTGGCGGAGCAGCAGCGAtt                                                                       >  3:570470/1‑150 (MQ=255)
                                                                              aTTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTc                                                                      <  4:545165/151‑1 (MQ=255)
                                                                              aTTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTc                                                                      <  3:1254331/151‑1 (MQ=255)
                                                                                 gaagaaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGGCTGGGCGGCGGGGGAGCGCGGGATTGCCTGTATTTCTGGGCGCGGAGAAGAGGTTTGGGGGGGGGGAGGAGCCGCGGTTTGtt                                                                   >  1:429720/1‑151 (MQ=255)
                                                                                   agaaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATTGGCGGCGGGGGGGCGCGGTGTTGCCCGTATTTCTGGCCGTGAAGAAGAGGCTTGGGCTGGGGGagaagcagc                                                                           >  3:400026/1‑141 (MQ=255)
                                                                                       ttAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGAGGGCGGGATTGCCTGTATTTCTGGCCGGGAAGAAGAGTTTCGGGGTGGGGTAGCAGCAGCGATTCGTTATGccc                                                             >  3:1325526/1‑151 (MQ=255)
                                                                                            aCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTgcg                                                        <  3:904214/151‑1 (MQ=255)
                                                                                            aCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTgcg                                                        <  3:575299/151‑1 (MQ=255)
                                                                                             cATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTgcgc                                                       <  3:786848/151‑1 (MQ=255)
                                                                                             cATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGCGGGGGGGGGCGGGATTGCCTGTATTTCTCGGCGGGGAGAAGAGGTTCGGGGTGGCGGAGGCGGAGCGGTTCGTTATTGCCGGgcg                                                        >  1:462526/1‑150 (MQ=255)
                                                                                               tgtgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACCGGGCGGGGGGGGGGGGCGGGATTGCCTGTATTTCTTGCCGTGGAGAAGAGGTTTGGGGGGGGGGAGCAGCAGCGATTCGTTATGCCCGTgcgc                                                       >  1:1132999/1‑149 (MQ=255)
                                                                                                gtgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTg                                                     <  3:1036213/150‑1 (MQ=255)
                                                                                                 tgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGGGGGGGGGGGGGGGGGGATTGCCTGTATTTCTGGCCGTGGAGAAGAGGTTTGGGGTGGCGGAGCAGCAGCGGTTCGGTATGGCCGTGCGCCggg                                                   >  3:182852/1‑151 (MQ=255)
                                                                                                  gCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGCGCGGGATTGCCTGTTTTTCTGGCCGGGAAGAAGGGGTTCGGGGGGGGGGAGGAGCAGCGATTTGTTTTGCCCGTGCGCTgg                                                    >  4:13260/1‑149 (MQ=255)
                                                                                                     ggCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGta                                                 <  3:770050/149‑1 (MQ=255)
                                                                                                      gCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGtaata                                              <  3:1071668/151‑1 (MQ=255)
                                                                                                       cGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGtaataa                                             <  3:567530/151‑1 (MQ=255)
                                                                                                       cGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGtaataa                                             <  3:1347873/151‑1 (MQ=255)
                                                                                                       cGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGtaataa                                             <  1:181023/151‑1 (MQ=255)
                                                                                                         aGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGGGGGGCGGGGGGGGGGGGTATTGCCTGTATTCCTTGGCGGGAAGAAGAGGTTTGGGATGGGGGAGCAGCAGCGATTTGTTATGCCCGTGCGCTGGGGAATaaa                                            >  4:75413/1‑150 (MQ=255)
                                                                                                            aataaACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATaaaaaaa                                        <  1:761165/151‑1 (MQ=255)
                                                                                                             ataaACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAAt                                       <  4:308776/151‑1 (MQ=255)
                                                                                                               aaaCTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTa                                     <  4:127514/151‑1 (MQ=255)
                                                                                                               aaaCTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGGAGGAGCAGCGATTCGTTATGCCCGGGCGCTGGGGAATAAAAAAAtt                                      >  2:193172/1‑150 (MQ=255)
                                                                                                                aaCTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTaa                                    <  3:552228/151‑1 (MQ=255)
                                                                                                                     gAACATACGCTGCACAATTTACCGGCGGGGGACTTGGCGGCGGGGGAGCGCGGGAATTCCTTTATTTCTGGGCGGGGAAAAAAGGTTCGGGGTGGGGGAGGAGCAGCCATTTGTTTTGTCCGTGGGGCGGGGGATaaaaaaaatta                                    >  4:201729/1‑142 (MQ=255)
                                                                                                                      aaCATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTca                              <  1:657803/151‑2 (MQ=255)
                                                                                                                      aaCATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCt                              <  1:836853/151‑1 (MQ=255)
                                                                                                                        cATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTgg                            <  3:1002435/151‑1 (MQ=255)
                                                                                                                         aTACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGGGGGGGGGCGCGGGATTTCCTGTATTCCTGGCCGGGGAGAAGAGGTTTGGGGTGGCGGAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGGAATAAAAAAATTAACTGTCTgg                            >  2:285964/1‑150 (MQ=255)
                                                                                                                         aTACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGGGCGGGATTGCCCGTATTCCTGGCCGTGGAGAAGAGGTTCGGGGTGGCGGAGCAACAGCGAATTGTTATGCCCGGGCGCTGGGGAATAAAAAAAATAAATGt                                >  3:725828/1‑146 (MQ=255)
                                                                                                                         aTACGCTGCACAATTTACCGGCGGGGGAATGGGCGGCGGGGGGGGGCGGGCTTGCCTTTATTCCTGGCCGGGAAGAAAAGTTTCGGGGTGGGGGAGCAGCAACGGTTCGGTATGCCCGGGCGCTGGGGAATAAAAAAAATAACTGTCTgg                            >  4:159440/1‑150 (MQ=255)
                                                                                                                                gCACAATTTACCGGCGGGTGGCTGGGCGGGGGGGGGGGGGGGGCATGCCTTTATTTCTGGCCGGGGAGAAGAGGTTTGGGGGGGGGGAGGAGGAGCGGTTTGGTTTGCCCGGGCGCTGGGGGATAAAAAAAATAAATGTTTGGTCGGTaa                     >  1:64021/1‑150 (MQ=255)
                                                                                                                                 cacaATTTACCGGCGGGTGACCGGGCGGCGGGGGGGGGCGGTATTGACTTTTTTTCTTGCCGCGGAGAAGAGTTTCGGGGGGGGGGAGCAGCAGCGGTTTGTTTTGCCCCTGCGCTGGGGAATAAAAAAATTAACTGTCTGGTGGGTaaa                    >  4:301990/1‑150 (MQ=255)
                                                                                                                                 cacaATTTACCGGCGGGGGACTGGGGGGCGGGGGAGCGCGGTATTTCCTGTATTTCTGGGCGGGAAGAAGAGATTCGGGATGGCGTAGCAGCAGCGATTCGTTTTGCCCGTGCGCCGGGGAATAAAAAAAATAAATGTCTGGTCgggcaa                    >  2:148294/1‑146 (MQ=255)
                                                                                                                                        tACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCgg             <  3:405965/150‑1 (MQ=255)
                                                                                                                                         aCCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGc            <  4:547163/150‑1 (MQ=255)
                                                                                                                                         aCCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCt           <  2:351068/151‑1 (MQ=255)
                                                                                                                                            ggcggGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGGTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTgg         >  1:539519/1‑150 (MQ=255)
                                                                                                                                             gcggGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGtt       <  1:342543/151‑1 (MQ=255)
                                                                                                                                                ggTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTca    <  3:743913/151‑1 (MQ=255)
                                                                                                                                                 gTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTcag   <  3:322208/151‑1 (MQ=255)
                                                                                                                                                  tGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTcagc  >  2:331915/1‑151 (MQ=255)
                                                                                                                                                  tGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTcagc  >  3:493606/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                       
GATTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGC  >  NC_000913/535755‑536051

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