Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913535,9010TG42.1% 119.5 / 13.7 78G64G (GGT→GGG) hyihydroxypyruvate isomerase
Reads supporting (aligned to +/- strand):  ref base T (14/30);  new base G (32/0);  total (48/30)
Fisher's exact test for biased strand distribution p-value = 9.11e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCC‑TGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTA‑TGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCA  >  NC_000913/535755‑536049
                                                                                                                                                  |                                                                                                                                                      
gATTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACt                                                                                                                                                    >  2:1281642/1‑151 (MQ=255)
 aTTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGgagaat                                                                                                                                                    >  2:912044/1‑148 (MQ=255)
   tttCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGagaggg                                                                                                                                                 >  2:727279/1‑146 (MQ=255)
     tctcGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGggcg                                                                                                                                               >  2:14030/1‑151 (MQ=255)
        cGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGggcg                                                                                                                                               <  2:35286/148‑1 (MQ=255)
         gCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGcggcgg                                                                                                                                           <  1:1042616/151‑1 (MQ=255)
         gCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGgggcgg                                                                                                                                           >  2:1224375/1‑151 (MQ=255)
           ccGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCgggg                                                                                                                                         <  1:859072/151‑1 (MQ=255)
           ccGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACCGGGGGGCgggg                                                                                                                                         >  1:287680/1‑151 (MQ=255)
           ccGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGGGCGGCgggg                                                                                                                                         >  2:1090355/1‑151 (MQ=255)
             gTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGa                                                                                                                                       <  2:1551815/151‑1 (MQ=255)
             gTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCgggggg                                                                                                                                       >  2:700020/1‑150 (MQ=255)
             gTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGCGGATgggggggggggg                                                                                                                                        >  1:561051/1‑150 (MQ=255)
                  gagaAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGGGGGGGGGGGGggcg                                                                                                                                  >  2:931672/1‑151 (MQ=255)
                  gagaAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGGGGGGGGGggcg                                                                                                                                  >  2:965181/1‑151 (MQ=255)
                   agaAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACCGGGGGGCGGGGGAGGGCgg                                                                                                                                 >  2:910267/1‑151 (MQ=255)
                    gaAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGCGCgg                                                                                                                                 >  2:298054/1‑150 (MQ=255)
                      aaGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAGTGGGGGGCGGGGGAGAGCggg                                                                                                                                >  1:396222/1‑148 (MQ=255)
                       aGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGGTGGGCGGGGGGGGGGCggggg                                                                                                                                >  1:1325797/1‑147 (MQ=255)
                          tGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGCGGGGGGGCGCGGGGTTGcc                                                                                                                          >  2:1135497/1‑151 (MQ=255)
                          tGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATCTGTGGGGGGGGGGCGTCGTATTgtc                                                                                                                          >  2:185442/1‑149 (MQ=255)
                          tGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGAGGGATGGGGGGGGGGGGGGGGGGGGATTGcc                                                                                                                          >  1:530932/1‑151 (MQ=255)
                          tGCGCAGGGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGTGGACTtgggggtgggggggggtggcaatgt                                                                                                                           >  1:647919/1‑137 (MQ=255)
                            cgcAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGCGGGCCGGGCGGGGGGGGAGCgggg                                                                                                                                 >  1:435235/1‑142 (MQ=255)
                              cAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGGGGCGGGGGGGGGCGGGAtt                                                                                                                             >  2:1242761/1‑144 (MQ=255)
                              cAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGATGGGGGGGGGGGGGGGGGGGGATTGCCTGt                                                                                                                       >  2:1062423/1‑150 (MQ=255)
                                    ggTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCC‑t                                                                                                                <  2:28476/150‑1 (MQ=255)
                                    ggTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCC‑Tg                                                                                                               <  2:590025/151‑1 (MQ=255)
                                     gTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGCGCGGGATTGCCTGTATTTC‑ttg                                                                                                              >  1:1045149/1‑149 (MQ=255)
                                         tCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCC‑TGGCCg                                                                                                           <  2:472558/150‑1 (MQ=255)
                                             ggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGGGGCGGGGGAGCGCGGGGTTGCCTGTATTTC‑TGGCCGt                                                                                                          >  2:1503382/1‑147 (MQ=255)
                                             ggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCC‑TGGCCGTgaag                                                                                                      <  1:655081/151‑1 (MQ=255)
                                               cgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGTGGGGGGGGGGGGGATTTCCTGTATTTC‑TTGGCGTgaag                                                                                                      >  1:1288419/1‑149 (MQ=255)
                                                gTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGGGCGGGATTGCCTGTATTCC‑TGGCCg                                                                                                           >  2:702759/1‑143 (MQ=255)
                                                  tGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCC‑TGGCCGTGaagaag                                                                                                   <  1:46381/149‑1 (MQ=255)
                                                   gAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCC‑TGGCCGTGAAGAAGAGt                                                                                                <  1:1489295/151‑1 (MQ=255)
                                                     aTTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGGATGGGGGGGGGGGGGGGGCGGGATTGCCTGTTTTTC‑TTGCCGGGGAGAAAAGgtt                                                                                              >  1:618319/1‑151 (MQ=255)
                                                            tttCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGCGGGGGGGCGCGGGATTGCCTGTTTTTC‑TGGCCGGGGAGAAGAGGGTCggg                                                                                          >  1:898295/1‑148 (MQ=255)
                                                            tttCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGGCGGGGAGCGCGGGTTTGGCTGTATTTC‑CGGCCGCGGAGAAAAGtttttg                                                                                           >  2:149346/1‑144 (MQ=255)
                                                               ccTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGGGCGGGATTGGCTGTATTTC‑TGGGCGGGGAGGAGAGGTTTGGGGTgggg                                                                                    >  1:874276/1‑149 (MQ=255)
                                                                 ttATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGTATGGGGGGGCGGGGAAGGCCGTATTTTCTGTATTCC‑TTGACGGGGTGAAAAGTTtggggttggggt                                                                                   >  1:306084/1‑150 (MQ=255)
                                                                   aTGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCC‑TGGCCGTGAAGAAGAGTTTCGGGATGGCGTag                                                                                 <  2:38331/150‑1 (MQ=255)
                                                                   aTGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCC‑TGGCCGTGAAGAAGAGTTTCGGGATGGCGTag                                                                                 <  1:204520/150‑1 (MQ=255)
                                                                     gACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCC‑TGGCCGTGAAGAAGAGTTTCGGGATGGCGTagcag                                                                              <  1:665100/151‑1 (MQ=255)
                                                                     gACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCC‑TGGCCGTGAAGAAGAGTTTCGGGATGGCGTagca                                                                               <  1:190372/150‑1 (MQ=255)
                                                                         aCGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGCGCGGGATTGCCTGTATTCC‑TGGCCGTGGAGAAGAGTTTTGGGGTGGGGTagcagc                                                                             >  2:1253568/1‑148 (MQ=255)
                                                                         aCGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCC‑TGGCCGTGAAGAAGAGTTTCGGGATGGCGTagcagcag                                                                           <  2:1493955/150‑1 (MQ=255)
                                                                         aCGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCC‑TGGCCGTGAAGAAGAGTTTCGGGATGGCGTagcagca                                                                            <  2:601960/149‑1 (MQ=255)
                                                                            aCATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCC‑TGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGAt                                                                       <  2:94710/151‑1 (MQ=255)
                                                                            aCATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGCGCGGGGGGGCGCGGGATTGGCCGTTTTCC‑TGGGCGTGTAGAAGAGTTTTGTGGGGGCGTaggagcag                                                                           >  2:1416093/1‑147 (MQ=255)
                                                                              aTTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGGGGCGGGGGAGCGCGGGATTGCCTGTATTCC‑CTGCCGTGGAGAAGAGTTTTGGGATGGCGTAGCAGCAGCGCtt                                                                      >  2:1165800/1‑150 (MQ=255)
                                                                                 gaagaaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCC‑TGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGt                                                                   <  2:504634/150‑1 (MQ=255)
                                                                                  aagaaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGCGCGGTATTGCCTGTATTCC‑TTGCCGTGAAGAAGAGGTTTGGGGTGGGGGAGCAGCAGCGATTCggt                                                                  >  2:201947/1‑148 (MQ=255)
                                                                                   agaaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGATGGGGGGGGGGGGGGGGGGGGATTGCCTGAATTTC‑TGGGCGGGGAGGAGAGGTTTGGGGGGGGGGAGCAGCAGCGGTTCggttt                                                                >  1:1218598/1‑147 (MQ=255)
                                                                                        tAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTTGGGGGCGGGGGGGCGGGGGATTGCCTGGATTCCTTGGCGGGAAAGAAGGGTTT‑GGGGTGGCGGGGCAGCAGGGGTTCGTTA‑TTccc                                                           >  1:1290296/1‑150 (MQ=255)
                                                                                          aaaCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGGGCGGGAATTCCTGTATTTC‑TGGGCGGGAAGAAGAGGTTTGGGGTGGCGGAGGAGCAGCGACTCGTTA‑TGCCCggg                                                        >  2:1034515/1‑149 (MQ=255)
                                                                                               tgtgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGGGGGGGGGGGGGGGGGCGCGGGATTTGCCGTATTCC‑CCGCCGCGGAGAGGAGGTTTGGGGGGGGGGAGCAGCAGCGGGTCGTTG‑TTCCCGTg                                                        >  1:1051923/1‑146 (MQ=255)
                                                                                                 tgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGGTGGGGGGGGGGGGAGAGCGGGGTTGGCTGTATTTC‑TGGGCGGGGAAAAGAGGTTTGGGGGGGGGGGGCAGGAGCGGTTTGTTA‑TTCCCGTGGGCTggg                                                 >  1:1553355/1‑151 (MQ=255)
                                                                                                 tgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACAGGGGGGCGGGGGGGGGCGGGATTGCCTTTGTTTC‑TGGGCGGGGAGAAGAGGTTTGGGGTGGGGTCGCCGCAGGGGATCGTTATTGCCCggcggcggg                                                  >  1:963921/1‑143 (MQ=255)
                                                                                                  gCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGGGGCGGGGGGGCGCGGGATTGCCTGTATTTC‑CTGCCGCGGAGAAGAGGTTTGGGGTGGCGTAGCAGCAGCCGTTCGTTA‑TGTCCGCGCGCTggg                                                 >  2:838125/1‑150 (MQ=255)
                                                                                                   cTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGCGCGGGATTGCCTGTTTTTC‑TGGGCGGGGAGAAGAGGTTTGGGGTGGGGTAGCAGCAGCGATTCGTTT‑TGCCCGTGCGCTggg                                                 >  1:1261236/1‑149 (MQ=255)
                                                                                                     ggCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCC‑TGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTA‑TGCCCGTGCGCTGGGtaat                                             <  2:1013579/151‑1 (MQ=255)
                                                                                                     ggCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCC‑TGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTA‑TGCCCGTGCGCTGGGtaa                                              <  1:177257/150‑1 (MQ=255)
                                                                                                         aGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCC‑TGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTA‑TGCCCGTGCGCTGGGTAATaaaa                                         <  2:1553505/151‑1 (MQ=255)
                                                                                                          gTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACGAGGCGGGGGGGGAGCGCGGGATTGCCTGTATTCC‑CGGCCCTGAAGGAGAGGTTCGGGGAGGGGCAGCAGCACCGATTTGTTT‑TTCCCGTGCGCTGGGaaaaaaaa                                         >  2:1345586/1‑150 (MQ=255)
                                                                                                          gTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGGGCGGCGGGGGGGCGCGGTATTGCCTGTATTCC‑TGGCCGGGAAGAAGAGGTTTGGGGTGGCGGAGCAGCAGCGGTTCGTTA‑TGCCCCTGCGCTGGGGAATaaaa                                         >  2:1034727/1‑150 (MQ=255)
                                                                                                               aaaCTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCC‑TGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTA‑TGCCCGTGCGCTGGGTAATAAAAAAAtt                                    <  1:1046440/150‑1 (MQ=255)
                                                                                                               aaaCTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCC‑TGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTA‑TGCCCGTGCGCTGGGTAATAAAAAAATTa                                   <  1:619639/151‑1 (MQ=255)
                                                                                                                    cGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCC‑TGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTA‑TGCCCGTGCGCTGGGTAATAAAAAAATTAACt                                <  2:718510/149‑1 (MQ=255)
                                                                                                                     gAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGCGCGGGATTTCCCGTATTCC‑TGGCCGTGAAGAAGAGTTTTGGGGTGGGGTAGCAGCAGCGATTCGTTA‑TTCCCGGGCGCTTGGGAATAAAAAAATTaa                                  >  2:548206/1‑146 (MQ=255)
                                                                                                                      aaCATACGCTGCACAATTTATCGGCGGGGGACCGGGGGGCGGGGGGGCGCGGGAATGCCCGTATTTC‑TTGCCGTGAAGAAGAGGTTTGGGGTGGGGTAGCAGCAACGATTCGTTT‑TGCCCGGGCGCTGGGGaaaaaaaaaaat                                    >  2:929842/1‑141 (MQ=255)
                                                                                                                      aaCATACGCTGCACAATTTACCGGCGGGCGCCTGGGCGGCGGGGGGGGGGGGGATTGCCTGTTTTTC‑TCGGCCTGGAGAAGAGGTTTCGGATGGGGGGGCAGCAACGAATAGTTA‑TCCCCGGGCGCTGGGGTATAAAAAAATTaaacg                               >  2:1275652/1‑145 (MQ=255)
                                                                                                                            cGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCC‑TGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTA‑TGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCgg                      <  1:278967/151‑1 (MQ=255)
                                                                                                                            cGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCC‑TGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTA‑TGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCgg                      <  2:849474/151‑1 (MQ=255)
                                                                                                                            cGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCC‑TGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTA‑TGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCgg                      <  1:1120862/151‑1 (MQ=255)
                                                                                                                               tGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGAGCGCGGGATTGCCTGTATTCC‑TGGCCGTGAAGAAGAGGTTCGGGGTGGCGGAGCAGCAGCGATTCGTTA‑TTCCCCTGCGCTGGGGAATAAAAAAATTAAATGTCTGGTCgg                      >  2:1148123/1‑148 (MQ=255)
                                                                                                                                   caATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCC‑TGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTA‑TGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACg               <  2:5097/151‑1 (MQ=255)
                                                                                                                                   caATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCC‑TGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTA‑TGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACg               <  2:1354200/151‑1 (MQ=255)
                                                                                                                                                ggTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCC‑TGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTA‑TGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTca  <  2:618397/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
GATTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCC‑TGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTA‑TGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCA  >  NC_000913/535755‑536049

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: ‑