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→GGG) hyihydroxypyruvate 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

GATTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGC  >  NC_000913/535755‑536051
                                                                                                                                                  |                                                                                                                                                       
gATTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑Ct                                                                                                                                                    >  4:373700/1‑151 (MQ=255)
gATTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGA‑Ct                                                                                                                                                    >  1:1173526/1‑151 (MQ=255)
 aTTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTg                                                                                                                                                   <  4:446171/151‑1 (MQ=255)
 aTTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTg                                                                                                                                                   <  1:1326147/151‑1 (MQ=255)
   tttCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTgg                                                                                                                                                  <  1:660457/150‑1 (MQ=255)
        cGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGG‑CTgggggggg                                                                                                                                            >  1:488094/1‑149 (MQ=255)
         gCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGA‑ATGgggggggg                                                                                                                                           >  4:197428/1‑151 (MQ=255)
           ccGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCgggg                                                                                                                                         >  2:538027/1‑151 (MQ=255)
           ccGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGG‑ATGGGGGGCgggg                                                                                                                                         >  1:840011/1‑151 (MQ=255)
           ccGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGG‑ATGGGCGGCgggg                                                                                                                                         >  3:1101439/1‑151 (MQ=255)
           ccGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGA‑CTGGGGGGCgggg                                                                                                                                         >  3:1335290/1‑151 (MQ=255)
            cGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGG‑CTgggggggggggg                                                                                                                                        >  3:958347/1‑151 (MQ=255)
            cGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGA‑CTGGGGGGCggggg                                                                                                                                        >  4:573687/1‑151 (MQ=255)
              ttttGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGG‑CTGGGGGGCggggggg                                                                                                                                      >  1:1234288/1‑149 (MQ=255)
               tttGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGA‑CTgggggggggggggg                                                                                                                                      >  3:1113255/1‑148 (MQ=255)
                ttGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAgc                                                                                                                                     <  2:177892/150‑1 (MQ=255)
                ttGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGG‑ATGGGGCGCggggg                                                                                                                                        >  2:443180/1‑147 (MQ=255)
                ttGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGA‑CTGGGGGGCGGGGGAggg                                                                                                                                    >  3:59478/1‑149 (MQ=255)
                  gagaAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGcgcg                                                                                                                                  <  2:59638/151‑1 (MQ=255)
                   agaAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGG‑CTGGGGGGCGGGGGGGGGCgg                                                                                                                                 >  1:405903/1‑151 (MQ=255)
                     aaaGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTTGGCGGGGGGGGGGCGCgg                                                                                                                                 >  3:363779/1‑149 (MQ=255)
                        gCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTg                                                                                                                            <  1:662904/151‑1 (MQ=255)
                          tGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGGGGGGGGGGGGCGGGGGATTGcc                                                                                                                          >  1:1344320/1‑151 (MQ=255)
                          tGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGATCTGGG‑GGCGGGGGGGGGCgggggggac                                                                                                                          >  1:269489/1‑144 (MQ=255)
                            cgcAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGG‑CGGGGCGGCGGGGGGGCGCgggagt                                                                                                                             >  1:1370185/1‑142 (MQ=255)
                             gcAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGt                                                                                                                       <  1:148195/151‑1 (MQ=255)
                              cAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTa                                                                                                                      <  3:1230825/151‑1 (MQ=255)
                                    ggTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCt                                                                                                                 <  1:1101307/150‑1 (MQ=255)
                                    ggTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTg                                                                                                                <  2:552706/151‑1 (MQ=255)
                                         tCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTTA‑CCTGGCGGGCGGGGGGGGCGGTATTGCCGGGAGTCCTg                                                                                                                >  3:216015/1‑146 (MQ=255)
                                            cggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTGCCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTgaa                                                                                                        <  3:195910/151‑1 (MQ=255)
                                            cggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGG‑CTGGGGGGCGGGGGGGGGCGGGGTTTGCTGTTTTcc                                                                                                                  >  1:907389/1‑141 (MQ=255)
                                            cggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGA‑CTGGGCGGCGGGGGGGCGCGGGATTGCCTGTATTTCTGGcc                                                                                                             >  2:512207/1‑146 (MQ=255)
                                             ggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGG‑CTGGGCGGCGGGGGGGGGCGGGATTTGCTGTTTTTCttggcg                                                                                                            >  4:539140/1‑141 (MQ=255)
                                             ggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGA‑CTGGGGGGCGGGGGGGCGGGGGATTGGCTGTATTTCTGGCCGTg                                                                                                          >  1:262404/1‑148 (MQ=255)
                                                gTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGaagaag                                                                                                    <  2:70933/151‑1 (MQ=255)
                                                   gAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGt                                                                                                 <  2:288968/151‑1 (MQ=255)
                                                   gAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGG‑CTGGGGGGGGGGGGAGCGCGGGATTGCCTGTTTTTCTGGGCggggagga                                                                                                     >  1:325832/1‑140 (MQ=255)
                                                     aTTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGttt                                                                                               <  3:625585/151‑1 (MQ=255)
                                                      tttATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGG‑CTCGGTGGGGGGGGGGCGCGGGGATGCCTCTTTTTCCGGCCGTGGaggaga                                                                                                   >  1:1356706/1‑146 (MQ=255)
                                                       ttATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTc                                                                                              <  3:1009353/150‑1 (MQ=255)
                                                       ttATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGA‑CTGGGGGGCGGGGGGGCGCGGGATTGCCTGTATTTCTGGCCGGGAAGAAGAGgtt                                                                                               >  3:481293/1‑149 (MQ=255)
                                                         aTGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCgg                                                                                            <  3:263760/150‑1 (MQ=255)
                                                          tGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGG‑CTGGGGGGCGGGGGAGGGCGGGATTGGCTGTATTCCTGGCCGGGGAGAAGAGGTTTgggg                                                                                          >  1:1027407/1‑150 (MQ=255)
                                                                cTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGG‑CTGGGGGGCGGGGGAGGGCGGGATTTCCTGTATTTCTGGCCGGGAAGAAGAGTTTTGGGGTgg                                                                                       >  4:194308/1‑148 (MQ=255)
                                                                 ttATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGG‑CTGGGGGGCGGGGGGGGGGGGGATTGGCTGTATTTCTGGCCGGGGAGAAGAGGTTTGGGGTgg                                                                                       >  4:300280/1‑147 (MQ=255)
                                                                    tGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGGGGCGGGGGGGCGCGGGATTGCCTGTATTTCTGGCCGTGGAGAAGAGGTTCGGGGTgggg                                                                                     >  2:151583/1‑144 (MQ=255)
                                                                     gACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTagca                                                                                <  3:741576/150‑1 (MQ=255)
                                                                      aCTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGGGGGGGGGGGGGGCGGGATTGCCTGTATTTCTGGCCGTGAAGAAGAGGTTCGGGGTGGGGTagca                                                                                >  3:47380/1‑149 (MQ=255)
                                                                       cTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGGGGCGGGGGGGGGCGGGATTGGCTTTATTTCTTGGCGGGAAGAAGAGGTTTgggggggggtgggag                                                                               >  4:253697/1‑137 (MQ=255)
                                                                         aCGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGG‑CTGGGGGGGGGGGGGGGGGGGGATTTCCTGTATTCCTGGGCGTGGAGAAAAGGTTTgggggg                                                                                        >  1:1319395/1‑135 (MQ=255)
                                                                          cGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGAGA‑ATGGGGGGCGGGGGGGGGCGGTGTTGCCTGTATTCCTGGCCGTGGAGGAGAGTTTCgggggggtggagca                                                                                >  4:183900/1‑134 (MQ=255)
                                                                            aCATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGGGAAGAAGAGTTTCGGGATGGCGTAGCAGCATCGAt                                                                        <  4:465581/151‑1 (MQ=255)
                                                                              aTTGAAGAATTAAAACATGTGCTGGCGAGTGATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGG‑CTGGGGGGCGGGGGGGCGCGGGATTGCCTGTATTCCTGGCCCTGGAGAAGAGGTTTGGGGGGGGGGAGCAGCAGCGAttt                                                                      >  1:927846/1‑150 (MQ=255)
                                                                              aTTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGGGCGCGGGATTGCCTGTATTTCTGGCCGTGGAGAAGAGGTTCGGGGTGGCGGAGCAGCAGCGAtt                                                                       >  3:570470/1‑150 (MQ=255)
                                                                              aTTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTc                                                                      <  4:545165/151‑1 (MQ=255)
                                                                              aTTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTc                                                                      <  3:1254331/151‑1 (MQ=255)
                                                                                 gaagaaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGG‑CTGGGCGGCGGGGGAGCGCGGGATTGCCTGTATTTCTGGGCGCGGAGAAGAGGTTTGGGGGGGGGGAGGAGCCGCGGTTTGtt                                                                   >  1:429720/1‑151 (MQ=255)
                                                                                   agaaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGA‑ATTGGCGGCGGGGGGGCGCGGTGTTGCCCGTATTTCTGGCCGTGAAGAAGAGGCTTGGGCTGGGGGagaagcagc                                                                           >  3:400026/1‑141 (MQ=255)
                                                                                       ttAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGA‑CTGGGCGGCGGGGGAGGGCGGGATTGCCTGTATTTCTGGCCGGGAAGAAGAGTTTCGGGGTGGGGTAGCAGCAGCGATTCGTTATGccc                                                             >  3:1325526/1‑151 (MQ=255)
                                                                                            aCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTgcg                                                        <  3:904214/151‑1 (MQ=255)
                                                                                            aCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTgcg                                                        <  3:575299/151‑1 (MQ=255)
                                                                                             cATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTgcgc                                                       <  3:786848/151‑1 (MQ=255)
                                                                                             cATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGG‑CTGGGGGGCGGGGGGGGGCGGGATTGCCTGTATTTCTCGGCGGGGAGAAGAGGTTCGGGGTGGCGGAGGCGGAGCGGTTCGTTATTGCCGGgcg                                                        >  1:462526/1‑150 (MQ=255)
                                                                                               tgtgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CCGGGCGGGGGGGGGGGGCGGGATTGCCTGTATTTCTTGCCGTGGAGAAGAGGTTTGGGGGGGGGGAGCAGCAGCGATTCGTTATGCCCGTgcgc                                                       >  1:1132999/1‑149 (MQ=255)
                                                                                                gtgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTg                                                     <  3:1036213/150‑1 (MQ=255)
                                                                                                 tgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGGGGGGGGGGGGGGGGGGATTGCCTGTATTTCTGGCCGTGGAGAAGAGGTTTGGGGTGGCGGAGCAGCAGCGGTTCGGTATGGCCGTGCGCCggg                                                   >  3:182852/1‑151 (MQ=255)
                                                                                                  gCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGA‑CTGGGGGGCGGGGGGGCGCGGGATTGCCTGTTTTTCTGGCCGGGAAGAAGGGGTTCGGGGGGGGGGAGGAGCAGCGATTTGTTTTGCCCGTGCGCTgg                                                    >  4:13260/1‑149 (MQ=255)
                                                                                                     ggCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGta                                                 <  3:770050/149‑1 (MQ=255)
                                                                                                      gCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGtaata                                              <  3:1071668/151‑1 (MQ=255)
                                                                                                       cGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGtaataa                                             <  3:567530/151‑1 (MQ=255)
                                                                                                       cGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGtaataa                                             <  3:1347873/151‑1 (MQ=255)
                                                                                                       cGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGtaataa                                             <  1:181023/151‑1 (MQ=255)
                                                                                                         aGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGA‑ATGGGGGGCGGGGGGGGGGGGTATTGCCTGTATTCCTTGGCGGGAAGAAGAGGTTTGGGATGGGGGAGCAGCAGCGATTTGTTATGCCCGTGCGCTGGGGAATaaa                                            >  4:75413/1‑150 (MQ=255)
                                                                                                            aataaACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATaaaaaaa                                        <  1:761165/151‑1 (MQ=255)
                                                                                                             ataaACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAAt                                       <  4:308776/151‑1 (MQ=255)
                                                                                                               aaaCTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTa                                     <  4:127514/151‑1 (MQ=255)
                                                                                                               aaaCTCGAACATACGCTGCACAATTTACCGGCGGGGGA‑CTGGGCGGCGGGGGGGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGGAGGAGCAGCGATTCGTTATGCCCGGGCGCTGGGGAATAAAAAAAtt                                      >  2:193172/1‑150 (MQ=255)
                                                                                                                aaCTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTaa                                    <  3:552228/151‑1 (MQ=255)
                                                                                                                     gAACATACGCTGCACAATTTACCGGCGGGGGA‑CTTGGCGGCGGGGGAGCGCGGGAATTCCTTTATTTCTGGGCGGGGAAAAAAGGTTCGGGGTGGGGGAGGAGCAGCCATTTGTTTTGTCCGTGGGGCGGGGGATaaaaaaaatta                                    >  4:201729/1‑142 (MQ=255)
                                                                                                                      aaCATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTca                              <  1:657803/151‑2 (MQ=255)
                                                                                                                      aaCATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCt                              <  1:836853/151‑1 (MQ=255)
                                                                                                                        cATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTgg                            <  3:1002435/151‑1 (MQ=255)
                                                                                                                         aTACGCTGCACAATTTACCGGCGGGGGG‑CTGGGGGGGGGGGGGGCGCGGGATTTCCTGTATTCCTGGCCGGGGAGAAGAGGTTTGGGGTGGCGGAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGGAATAAAAAAATTAACTGTCTgg                            >  2:285964/1‑150 (MQ=255)
                                                                                                                         aTACGCTGCACAATTTACCGGCGGGGGA‑CTGGGCGGCGGGGGGGGGCGGGATTGCCCGTATTCCTGGCCGTGGAGAAGAGGTTCGGGGTGGCGGAGCAACAGCGAATTGTTATGCCCGGGCGCTGGGGAATAAAAAAAATAAATGt                                >  3:725828/1‑146 (MQ=255)
                                                                                                                         aTACGCTGCACAATTTACCGGCGGGGGA‑ATGGGCGGCGGGGGGGGGCGGGCTTGCCTTTATTCCTGGCCGGGAAGAAAAGTTTCGGGGTGGGGGAGCAGCAACGGTTCGGTATGCCCGGGCGCTGGGGAATAAAAAAAATAACTGTCTgg                            >  4:159440/1‑150 (MQ=255)
                                                                                                                                gCACAATTTACCGGCGGGTGG‑CTGGGCGGGGGGGGGGGGGGGGCATGCCTTTATTTCTGGCCGGGGAGAAGAGGTTTGGGGGGGGGGAGGAGGAGCGGTTTGGTTTGCCCGGGCGCTGGGGGATAAAAAAAATAAATGTTTGGTCGGTaa                     >  1:64021/1‑150 (MQ=255)
                                                                                                                                 cacaATTTACCGGCGGGTGA‑CCGGGCGGCGGGGGGGGGCGGTATTGACTTTTTTTCTTGCCGCGGAGAAGAGTTTCGGGGGGGGGGAGCAGCAGCGGTTTGTTTTGCCCCTGCGCTGGGGAATAAAAAAATTAACTGTCTGGTGGGTaaa                    >  4:301990/1‑150 (MQ=255)
                                                                                                                                 cacaATTTACCGGCGGGGGA‑CTGGGGGGCGGGGGAGCGCGGTATTTCCTGTATTTCTGGGCGGGAAGAAGAGATTCGGGATGGCGTAGCAGCAGCGATTCGTTTTGCCCGTGCGCCGGGGAATAAAAAAAATAAATGTCTGGTCgggcaa                    >  2:148294/1‑146 (MQ=255)
                                                                                                                                        tACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCgg             <  3:405965/150‑1 (MQ=255)
                                                                                                                                         aCCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGc            <  4:547163/150‑1 (MQ=255)
                                                                                                                                         aCCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCt           <  2:351068/151‑1 (MQ=255)
                                                                                                                                            ggcggGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGGTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTgg         >  1:539519/1‑150 (MQ=255)
                                                                                                                                             gcggGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGtt       <  1:342543/151‑1 (MQ=255)
                                                                                                                                                ggTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTca    <  3:743913/151‑1 (MQ=255)
                                                                                                                                                 gTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTcag   <  3:322208/151‑1 (MQ=255)
                                                                                                                                                  tGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTcagc  >  2:331915/1‑151 (MQ=255)
                                                                                                                                                  tGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTcagc  >  3:493606/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                       
GATTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGA‑CTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGC  >  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: ‑