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→GGGhyihydroxypyruvate 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

GATTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCA  >  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)
                                    ggTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCt                                                                                                                <  2:28476/150‑1 (MQ=255)
                                    ggTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTg                                                                                                               <  2:590025/151‑1 (MQ=255)
                                     gTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGCGCGGGATTGCCTGTATTTCttg                                                                                                              >  1:1045149/1‑149 (MQ=255)
                                         tCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCg                                                                                                           <  2:472558/150‑1 (MQ=255)
                                             ggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGGGGCGGGGGAGCGCGGGGTTGCCTGTATTTCTGGCCGt                                                                                                          >  2:1503382/1‑147 (MQ=255)
                                             ggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTgaag                                                                                                      <  1:655081/151‑1 (MQ=255)
                                               cgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGTGGGGGGGGGGGGGATTTCCTGTATTTCTTGGCGTgaag                                                                                                      >  1:1288419/1‑149 (MQ=255)
                                                gTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGGGCGGGATTGCCTGTATTCCTGGCCg                                                                                                           >  2:702759/1‑143 (MQ=255)
                                                  tGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGaagaag                                                                                                   <  1:46381/149‑1 (MQ=255)
                                                   gAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGt                                                                                                <  1:1489295/151‑1 (MQ=255)
                                                     aTTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGGATGGGGGGGGGGGGGGGGCGGGATTGCCTGTTTTTCTTGCCGGGGAGAAAAGgtt                                                                                              >  1:618319/1‑151 (MQ=255)
                                                            tttCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGCGGGGGGGCGCGGGATTGCCTGTTTTTCTGGCCGGGGAGAAGAGGGTCggg                                                                                          >  1:898295/1‑148 (MQ=255)
                                                            tttCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGGCGGGGAGCGCGGGTTTGGCTGTATTTCCGGCCGCGGAGAAAAGtttttg                                                                                           >  2:149346/1‑144 (MQ=255)
                                                               ccTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGGGCGGGATTGGCTGTATTTCTGGGCGGGGAGGAGAGGTTTGGGGTgggg                                                                                    >  1:874276/1‑149 (MQ=255)
                                                                 ttATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGTATGGGGGGGCGGGGAAGGCCGTATTTTCTGTATTCCTTGACGGGGTGAAAAGTTtggggttggggt                                                                                   >  1:306084/1‑150 (MQ=255)
                                                                   aTGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTag                                                                                 <  2:38331/150‑1 (MQ=255)
                                                                   aTGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTag                                                                                 <  1:204520/150‑1 (MQ=255)
                                                                     gACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTagcag                                                                              <  1:665100/151‑1 (MQ=255)
                                                                     gACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTagca                                                                               <  1:190372/150‑1 (MQ=255)
                                                                         aCGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGCGCGGGATTGCCTGTATTCCTGGCCGTGGAGAAGAGTTTTGGGGTGGGGTagcagc                                                                             >  2:1253568/1‑148 (MQ=255)
                                                                         aCGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTagcagcag                                                                           <  2:1493955/150‑1 (MQ=255)
                                                                         aCGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTagcagca                                                                            <  2:601960/149‑1 (MQ=255)
                                                                            aCATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGAt                                                                       <  2:94710/151‑1 (MQ=255)
                                                                            aCATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGCGCGGGGGGGCGCGGGATTGGCCGTTTTCCTGGGCGTGTAGAAGAGTTTTGTGGGGGCGTaggagcag                                                                           >  2:1416093/1‑147 (MQ=255)
                                                                              aTTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGGGGCGGGGGAGCGCGGGATTGCCTGTATTCCCTGCCGTGGAGAAGAGTTTTGGGATGGCGTAGCAGCAGCGCtt                                                                      >  2:1165800/1‑150 (MQ=255)
                                                                                 gaagaaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGt                                                                   <  2:504634/150‑1 (MQ=255)
                                                                                  aagaaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGCGCGGTATTGCCTGTATTCCTTGCCGTGAAGAAGAGGTTTGGGGTGGGGGAGCAGCAGCGATTCggt                                                                  >  2:201947/1‑148 (MQ=255)
                                                                                   agaaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGATGGGGGGGGGGGGGGGGGGGGATTGCCTGAATTTCTGGGCGGGGAGGAGAGGTTTGGGGGGGGGGAGCAGCAGCGGTTCggttt                                                                >  1:1218598/1‑147 (MQ=255)
                                                                                        tAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTTGGGGGCGGGGGGGCGGGGGATTGCCTGGATTCCTTGGCGGGAAAGAAGGGTTTGGGGTGGCGGGGCAGCAGGGGTTCGTTATTccc                                                           >  1:1290296/1‑150 (MQ=255)
                                                                                          aaaCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGGGCGGGAATTCCTGTATTTCTGGGCGGGAAGAAGAGGTTTGGGGTGGCGGAGGAGCAGCGACTCGTTATGCCCggg                                                        >  2:1034515/1‑149 (MQ=255)
                                                                                               tgtgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGGGGGGGGGGGGGGGGGCGCGGGATTTGCCGTATTCCCCGCCGCGGAGAGGAGGTTTGGGGGGGGGGAGCAGCAGCGGGTCGTTGTTCCCGTg                                                        >  1:1051923/1‑146 (MQ=255)
                                                                                                 tgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGGTGGGGGGGGGGGGAGAGCGGGGTTGGCTGTATTTCTGGGCGGGGAAAAGAGGTTTGGGGGGGGGGGGCAGGAGCGGTTTGTTATTCCCGTGGGCTggg                                                 >  1:1553355/1‑151 (MQ=255)
                                                                                                 tgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACAGGGGGGCGGGGGGGGGCGGGATTGCCTTTGTTTCTGGGCGGGGAGAAGAGGTTTGGGGTGGGGTCGCCGCAGGGGATCGTTATTGCCCggcggcggg                                                  >  1:963921/1‑143 (MQ=255)
                                                                                                  gCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGGGGCGGGGGGGCGCGGGATTGCCTGTATTTCCTGCCGCGGAGAAGAGGTTTGGGGTGGCGTAGCAGCAGCCGTTCGTTATGTCCGCGCGCTggg                                                 >  2:838125/1‑150 (MQ=255)
                                                                                                   cTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGCGCGGGATTGCCTGTTTTTCTGGGCGGGGAGAAGAGGTTTGGGGTGGGGTAGCAGCAGCGATTCGTTTTGCCCGTGCGCTggg                                                 >  1:1261236/1‑149 (MQ=255)
                                                                                                     ggCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGtaat                                             <  2:1013579/151‑1 (MQ=255)
                                                                                                     ggCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGtaa                                              <  1:177257/150‑1 (MQ=255)
                                                                                                         aGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATaaaa                                         <  2:1553505/151‑1 (MQ=255)
                                                                                                          gTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACGAGGCGGGGGGGGAGCGCGGGATTGCCTGTATTCCCGGCCCTGAAGGAGAGGTTCGGGGAGGGGCAGCAGCACCGATTTGTTTTTCCCGTGCGCTGGGaaaaaaaa                                         >  2:1345586/1‑150 (MQ=255)
                                                                                                          gTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGGGCGGCGGGGGGGCGCGGTATTGCCTGTATTCCTGGCCGGGAAGAAGAGGTTTGGGGTGGCGGAGCAGCAGCGGTTCGTTATGCCCCTGCGCTGGGGAATaaaa                                         >  2:1034727/1‑150 (MQ=255)
                                                                                                               aaaCTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAAtt                                    <  1:1046440/150‑1 (MQ=255)
                                                                                                               aaaCTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTa                                   <  1:619639/151‑1 (MQ=255)
                                                                                                                    cGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACt                                <  2:718510/149‑1 (MQ=255)
                                                                                                                     gAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGCGCGGGATTTCCCGTATTCCTGGCCGTGAAGAAGAGTTTTGGGGTGGGGTAGCAGCAGCGATTCGTTATTCCCGGGCGCTTGGGAATAAAAAAATTaa                                  >  2:548206/1‑146 (MQ=255)
                                                                                                                      aaCATACGCTGCACAATTTATCGGCGGGGGACCGGGGGGCGGGGGGGCGCGGGAATGCCCGTATTTCTTGCCGTGAAGAAGAGGTTTGGGGTGGGGTAGCAGCAACGATTCGTTTTGCCCGGGCGCTGGGGaaaaaaaaaaat                                    >  2:929842/1‑141 (MQ=255)
                                                                                                                      aaCATACGCTGCACAATTTACCGGCGGGCGCCTGGGCGGCGGGGGGGGGGGGGATTGCCTGTTTTTCTCGGCCTGGAGAAGAGGTTTCGGATGGGGGGGCAGCAACGAATAGTTATCCCCGGGCGCTGGGGTATAAAAAAATTaaacg                               >  2:1275652/1‑145 (MQ=255)
                                                                                                                            cGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCgg                      <  1:278967/151‑1 (MQ=255)
                                                                                                                            cGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCgg                      <  2:849474/151‑1 (MQ=255)
                                                                                                                            cGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCgg                      <  1:1120862/151‑1 (MQ=255)
                                                                                                                               tGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGAGCGCGGGATTGCCTGTATTCCTGGCCGTGAAGAAGAGGTTCGGGGTGGCGGAGCAGCAGCGATTCGTTATTCCCCTGCGCTGGGGAATAAAAAAATTAAATGTCTGGTCgg                      >  2:1148123/1‑148 (MQ=255)
                                                                                                                                   caATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACg               <  2:5097/151‑1 (MQ=255)
                                                                                                                                   caATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACg               <  2:1354200/151‑1 (MQ=255)
                                                                                                                                                ggTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTca  <  2:618397/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
GATTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCA  >  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: