Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913535,9010TG43.0% 126.0 / 21.6 86G64G (GGT→GGGhyihydroxypyruvate isomerase
Reads supporting (aligned to +/- strand):  ref base T (12/37);  new base G (37/0);  total (49/37)
Fisher's exact test for biased strand distribution p-value = 3.68e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCA  >  NC_000913/535757‑536049
                                                                                                                                                |                                                                                                                                                     
tgttCTCGCCCGTTTTGATAAAGCTGCGCAGTGTAGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCGGTCTAGTAATAAACTCGAACAGACGCTGCACAATTTACCGGCGGGTGACt                                                                                                                                                   <  2:1074844/147‑1 (MQ=255)
 tttCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTgg                                                                                                                                                 <  2:882282/150‑1 (MQ=255)
   tctcGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGggc                                                                                                                                               <  1:1131516/150‑1 (MQ=255)
    ctcGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGggcgg                                                                                                                                             >  2:1321960/1‑151 (MQ=255)
      cGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTggggg                                                                                                                                              >  2:1061449/1‑146 (MQ=255)
           gTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGa                                                                                                                                      <  2:568312/151‑1 (MQ=255)
           gTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGa                                                                                                                                      <  2:276082/151‑1 (MQ=255)
            ttttGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGAGTgggggggggggg                                                                                                                                       >  2:659346/1‑149 (MQ=255)
            ttttGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGa                                                                                                                                      <  1:552539/150‑1 (MQ=255)
             tttGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGCGGGGGAg                                                                                                                                     >  1:1359131/1‑150 (MQ=255)
             tttGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCggggggg                                                                                                                                     >  1:148190/1‑148 (MQ=255)
             tttGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCggggg                                                                                                                                       >  1:353685/1‑148 (MQ=255)
              ttGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAgc                                                                                                                                    <  2:109075/150‑1 (MQ=255)
              ttGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAgc                                                                                                                                    <  1:974414/150‑1 (MQ=255)
                gagaAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAgcgc                                                                                                                                  <  2:153013/150‑1 (MQ=255)
                 agaAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCggggggggg                                                                                                                                   >  1:1343337/1‑144 (MQ=255)
                 agaAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATTGGGGGGGGGGGGGCGCgg                                                                                                                                >  1:1041744/1‑151 (MQ=255)
                  gaAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGt                                                                                                                               <  2:387901/151‑1 (MQ=255)
                   aaaGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTTGGCGGCGGGGGAGAGCgg                                                                                                                                >  1:367225/1‑149 (MQ=255)
                           gcAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGGGGGGGGGGGCGGGGTTGCCTGt                                                                                                                     >  2:578225/1‑151 (MQ=255)
                             aGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTAt                                                                                                                   <  2:110785/151‑1 (MQ=255)
                              gtgtGGTTTTCGCGGCGTTGAATTTCTGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGTGGGATGGGCGGCGGGGGGGGGGGGTATTTTCTGTAtt                                                                                                                  >  2:190529/1‑151 (MQ=255)
                              gtgtGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGGGGGGCGGGGGGGGGCGGGGTTGCCCGTCtt                                                                                                                  >  1:274097/1‑151 (MQ=255)
                                  ggTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGGGGGGCGGGGGAGGGGGGGGATGCCTTTAtttc                                                                                                                >  2:463643/1‑147 (MQ=255)
                                   gTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTgg                                                                                                             <  1:109730/151‑1 (MQ=255)
                                    ttttCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGCGCGGGATTGCCTGTATTCCTgg                                                                                                             >  1:545715/1‑150 (MQ=255)
                                       tCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGccc                                                                                                          <  2:1191022/150‑2 (MQ=255)
                                          cggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGGCTTGGCGGCGGGGGGGCGGGGTATTTTCTGTATTTCTGGGCggg                                                                                                        >  2:1148928/1‑147 (MQ=255)
                                          cggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTTGGCGGGGGGGGGGCGCGGTATTGCCTGTATTCCt                                                                                                               >  1:969597/1‑142 (MQ=255)
                                           ggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGGGCGGCGGGGGAGGGCGGGATTGCCCGTATTTCTGGCCggg                                                                                                        >  1:711327/1‑146 (MQ=255)
                                            gcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGCGCGGGATTGGCTGTTTTTCTTGCCggggag                                                                                                     >  2:553599/1‑145 (MQ=255)
                                             cgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTTGAACATCCGCTGCACAATTTACCGGCGGGGGACTGAGCGGGGGGCGACCGCGGTAATGGCTggtatgctgggc                                                                                                            >  2:693144/1‑132 (MQ=255)
                                              gTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGaagaag                                                                                                  <  1:514963/151‑1 (MQ=255)
                                               ttGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGCGGCGGGGGGGGGGGGGATTGCCTGTATTTCTGGGCGGgaagga                                                                                                   >  1:352293/1‑147 (MQ=255)
                                               ttGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGGGGGGGGGGGCGGGATTTCCTGTCTTTCTGGGCGGGaagaag                                                                                                  >  2:1000056/1‑150 (MQ=255)
                                                tGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGGATGGGGGGCGGGGGGGCGCGGGATTGCCTGGATTTCTGGCCGTg                                                                                                        >  2:454319/1‑143 (MQ=255)
                                                tGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGGGGGGGGGCGCGGGATTTCCTGTATTTCCGGCCggg                                                                                                        >  2:923722/1‑141 (MQ=255)
                                                 gAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGGGCGGCGGGGGAGCGCGGGATTGCCTGTATTTCTTGCCGTGAAGGagag                                                                                                >  2:100862/1‑150 (MQ=255)
                                                   aTTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTTGGGGGGGGGGGGGAGCGGTATTTCGTGTCTTCCTTGCCGGGAAGAATAGgtt                                                                                             >  1:1086307/1‑151 (MQ=255)
                                                       aTGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCggg                                                                                         <  1:582428/151‑1 (MQ=255)
                                                       aTGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCggg                                                                                         <  2:327438/151‑1 (MQ=255)
                                                         gTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGa                                                                                        <  1:1283062/150‑1 (MQ=255)
                                                         gTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGATGAGGCGGCGGGGGAGGGCGCGATTGCCTGCGTTGCTGGGCGTgaggagaaggggcgggg                                                                                        >  2:779127/1‑138 (MQ=255)
                                                          tttCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGCGCGGGATTGCCTGTATTTCTGGCCGGGGAGAAGAGGTTCGGGGTg                                                                                      >  1:120891/1‑151 (MQ=255)
                                                            tCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATgg                                                                                     <  1:843691/150‑1 (MQ=255)
                                                               ttATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGGCTGGGCGGCGGGGGAGGGCGGGATTGCCTGTCTTTCTGGCCGTGGAGAAGAGATTCGGGGTGGCg                                                                                   >  1:406621/1‑149 (MQ=255)
                                                                tATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGAATGGGGGGGGGGGGGGAGCGGGATTTTCTGTATTTCTTGCCGTTAAGAAGAGGTTCtggggtgcg                                                                                   >  2:1238610/1‑142 (MQ=255)
                                                                tATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGCGCGGGTTTTCCTGTTTTTCTGGCCGTGAAGAAGAGTTTTGGGGTgggg                                                                                   >  1:1289593/1‑146 (MQ=255)
                                                                 aTGATTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACATTTTACCGGCGGGTGACTGGGCGGCGAGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTa                                                                                 <  2:25502/149‑1 (MQ=255)
                                                                  tGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGGATGGGGGGCGGGGGTGCGCGGGCTTCCCCGTTTTTCTGGCCGTGGAGTAGAGGTTTgggggggggggggg                                                                              >  1:125249/1‑140 (MQ=255)
                                                                    aCTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGGGCGGCGGGGGGGCGCGGGATTTCCTGTATTCCTGGGCGGGGAGAAGAGGTTTGGGGTgggg                                                                                   >  1:656703/1‑142 (MQ=255)
                                                                        cGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGAGCGCGGCATTTCCTGTATTACTGGGCGGGGAGAAGAGTTTCGGGAggggggtgccgc                                                                            >  1:590886/1‑138 (MQ=255)
                                                                           cATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGAtt                                                                     <  2:665088/151‑1 (MQ=255)
                                                                              tGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGGGGGGGGGCGCGGGATTGCCTGTATTTCTGGCCCTGAAGAAGAGTTTCGGGGTGGCGGAGCAGCAGCGAtt                                                                     >  1:368881/1‑148 (MQ=255)
                                                                                aagaaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGAGCGGGGGATTTTCTGTATTTCTGGGCGTGGAGGAGAGGTTCGGGGGGGGGTAGCAGCAGCggatagttt                                                                >  2:1254320/1‑143 (MQ=255)
                                                                                    aTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGcc                                                            <  2:1090743/151‑1 (MQ=255)
                                                                                      tAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGATGGGCGGCGGGGGAGCGCGGGAATGACTTTATTTCCGGCCGTGGAGAAGAGGTTCGGGGGGGCGGAGCAGCAGCGATTCGTTATGccc                                                           >  2:569794/1‑150 (MQ=255)
                                                                                       aaaaCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGGGCGGGATTTCCCGTTTTTACTGGCGTGAAGAAGAGGTTCGGGATGGGGGAACAGCACAGATTTCTTTTGcccct                                                         >  1:952227/1‑149 (MQ=255)
                                                                                         aaCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTgc                                                       <  1:867604/151‑1 (MQ=255)
                                                                                              gtgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTgg                                                  <  1:493867/151‑1 (MQ=255)
                                                                                              gtgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTg                                                   <  1:1015029/150‑1 (MQ=255)
                                                                                              gtgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGCGGGGGGTTGGCCGTATTCCTGGCCCTGAAGAAGAGTTTCGGGGTGGCGGAGCCGGAGCGGATCGTTTTTGCCGTgcgc                                                     >  1:95302/1‑148 (MQ=255)
                                                                                              gtgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGGGCGGCGGGGGAGCGCGGGATTTCCTTTATTCCTGGCCGGGAAGAAGAGGTTCGGGGTGGGGGAGGAGCAGCGATTCGtttttgcc                                                           >  2:1401376/1‑136 (MQ=255)
                                                                                                gCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGGGCGGGGGCGGGGGGCGGGATTTTCTGTATTCCTGGCCGTGGAAAAGAATTTCGGGGGGGGGTAGCCACAGGGATTCGGTATGCCCGGGCGCggg                                                  >  2:294440/1‑149 (MQ=255)
                                                                                                 cTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCAGGGCGGCGGGGGGGGGCGGGGTTGTTTTTATTTCCGGCCGCGGGGAAGAGGTTCGGGACGGCGGAGCCGCAGCGGGTCGTTATGCCCGTGCGGTggg                                                 >  2:216502/1‑149 (MQ=255)
                                                                                                   ggCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATTGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGtaat                                             <  2:53030/151‑1 (MQ=255)
                                                                                                   ggCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGta                                               <  1:1181681/149‑1 (MQ=255)
                                                                                                    gCGAGTAATAAACTCGAACATACGGTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGtaata                                            <  2:479283/151‑1 (MQ=255)
                                                                                                    gCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGtaata                                            <  2:317896/151‑1 (MQ=255)
                                                                                                     cGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGGGCGGGATTGCCCGTTTTTCTGGCCGGGGAGAAGGGGTTCGGGGTGGGGGAGCAGCAGCGGTTTGTTATGCCCGTGCGCTGGGgaata                                            >  1:1226032/1‑150 (MQ=255)
                                                                                                      gAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGGGCGGGATTTCCTGGATTTCTGGGCGGGGAAGAGAGTTTCAGGGTTGGGGAGCAGCGGGGATTCATTATGGCCgggtggcggggagtaa                                           >  1:1000366/1‑135 (MQ=255)
                                                                                                        gTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATaaaaa                                        <  2:604239/151‑1 (MQ=255)
                                                                                                        gTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATaaaaa                                        <  2:1237069/151‑1 (MQ=255)
                                                                                                             aaaCTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTa                                   <  1:749309/151‑1 (MQ=255)
                                                                                                             aaaCTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGCGCGGGATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGGTGGCGGAGCAGCAGCGATTCGTTATGCCCGGGCGCTGGGGAATAAAAAAAtt                                    >  1:74090/1‑150 (MQ=255)
                                                                                                                 tCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGGGGGGGGGGGGCGGGGTATTTCCTGTTTTCTTTGCCGTGGAGAAGAGATTCGGGGTGGGGCAGCAGGAGCGATTTGTTTTGCCCGTGGGCTGGGTTATAAAAAAAATaaat                                >  2:31602/1‑148 (MQ=255)
                                                                                                                     aCATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTg                           <  1:717391/151‑1 (MQ=255)
                                                                                                                         aCGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTc                        <  1:971102/150‑1 (MQ=255)
                                                                                                                           gCTGCACAATTTACCGGCGGGGGGCTGGGCGGCGGGGGGGCGCGGGATTTCCCGTATTCCTCGCCGCGGAGAAGAGGTTCGGGGTGGGGGAGCAGCAGCGCCTCGTTTTGCCCGTGCGCTGGGGAATAAAAAAATTAACTGTCTg                           >  1:338373/1‑145 (MQ=255)
                                                                                                                            cTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGAGCGCGGGATTGCCTGTATTCCTGGCCGGGGAGAAGAGGTTTGGGGTGGGGTAGCAGCAGCGATTCGGTATGCCCGTGCGCTGGGGAATAAAAAAAATAACTGTCTTGTCgg                      >  1:961456/1‑149 (MQ=255)
                                                                                                                                acaATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAAc                <  1:1183430/151‑1 (MQ=255)
                                                                                                                                       aCCGGCGGGTGACTGGGCGGCGGGGGGGCGCGGGCTTTCCTGTCTTTCTGGGCGGGGAGAAAAGGTTTGGGGGGGGGGAGCAGCGGCGATTTTTTATGCCCGTGGGCTGGGGGAAAAAAAAAATAAATGTCTGGTCGGGGAAAAGCCgg           >  2:450672/1‑149 (MQ=255)
                                                                                                                                        ccGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTg        <  1:1179607/151‑1 (MQ=255)
                                                                                                                                          ggcggGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGt      >  2:221650/1‑151 (MQ=255)
                                                                                                                                          ggcggGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGt      <  1:97747/151‑1 (MQ=255)
                                                                                                                                              ggTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTca  <  2:29422/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
TTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCA  >  NC_000913/535757‑536049

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