Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,508,6570AG46.2% 77.7 / 22.2 65F386S (TTC→TCC) yhfWphosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base A (8/27);  new base G (30/0);  total (38/27)
Fisher's exact test for biased strand distribution p-value = 1.69e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCC  >  NC_000913/3508508‑3508806
                                                                                                                                                       |                                                                                                                                                     
cAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGaa                                                                                                                                                      <  1:366674/151‑1 (MQ=255)
   tCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGgaa                                                                                                                                                     >  1:820506/1‑149 (MQ=255)
    ccACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGaaaaa                                                                                                                                                   >  2:947770/1‑150 (MQ=255)
        gttgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGCGCGTGGgaaaaa                                                                                                                                                  >  2:197749/1‑146 (MQ=255)
          tgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGGAAAAAttaaa                                                                                                                                              >  1:1132320/1‑147 (MQ=255)
          tgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATcaaa                                                                                                                                             >  2:716774/1‑148 (MQ=255)
            tttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCAACCATTTTGTGGCGGTGGCGCGCGGGAAAAATcacaca                                                                                                                                           >  2:852480/1‑150 (MQ=255)
                tAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACgg                                                                                                                                        <  1:258349/149‑1 (MQ=255)
                tAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGt                                                                                                                                       <  1:719665/150‑1 (MQ=255)
                 aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTa                                                                                                                                      <  2:609867/150‑1 (MQ=255)
                 aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAg                                                                                                                                     <  2:386976/151‑1 (MQ=255)
                  aaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACAAACGGt                                                                                                                                       >  1:729800/1‑148 (MQ=255)
                       ggAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGGAGccccc                                                                                                                                >  2:952005/1‑150 (MQ=255)
                        gagagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCcac                                                                                                                              <  2:208146/151‑1 (MQ=255)
                            gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCa                                                                                                                          <  1:1197623/151‑1 (MQ=255)
                             aTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAACACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGCGGGCAAAATTACCCACAGTCGGCCCCCCAACa                                                                                                                          >  2:703164/1‑150 (MQ=255)
                             aTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACAAACGGGAGccccccc                                                                                                                              >  2:573231/1‑144 (MQ=255)
                                aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAATTTACACACGGTAGCCCCCCCATCaaaa                                                                                                                       >  2:478469/1‑150 (MQ=255)
                                 atatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACAAACGGGAGCCCCCCCCTCAGcaa                                                                                                                      >  1:485691/1‑150 (MQ=255)
                                   atGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGc                                                                                                                    <  1:57481/150‑1 (MQ=255)
                                   atGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGc                                                                                                                    <  1:101970/150‑1 (MQ=255)
                                           gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTgcg                                                                                                            <  1:286206/150‑1 (MQ=255)
                                             tCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTTACACACGGGAGCCCCCCCATCAGAAAACgtgggggg                                                                                                            >  2:1161258/1‑144 (MQ=255)
                                              cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACGGGAGCCCCCCCATCAGAAAGAGGGGGGGGcac                                                                                                         >  2:808937/1‑150 (MQ=255)
                                               aTGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGGAAAAATTACACACGGTAGCCCCCCCATCAGAAAACGCGGGgcgcgc                                                                                                         >  1:982292/1‑147 (MQ=255)
                                                tGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCAAACACGGGAGCCCCCACATCAGAAAACGGGGGGGGGACAcc                                                                                                      >  2:96865/1‑151 (MQ=255)
                                                    cGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGc                                                                                                  <  1:154327/151‑1 (MQ=255)
                                                    cGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGc                                                                                                  <  2:429779/151‑1 (MQ=255)
                                                        ggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCACATCAGAAAGCGGGGGGCGCACACCGAGCTg                                                                                                >  2:113218/1‑149 (MQ=255)
                                                        ggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACCGTAGCCCCCACATCAGAAAGCGTGGGGCGCACACCCAGCTGc                                                                                               >  2:1178490/1‑150 (MQ=255)
                                                               cccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACAAAACGTAGCCCCCACCACAGAAAACGTGGTGCGCACCCCGGGCTg                                                                                                >  2:721341/1‑142 (MQ=255)
                                                                  gTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGCGGGAAAAATTAAACACAGTAGCCCCCACACCAGAAAACGTGGGGCGCCCACCCAgcggcggggcg                                                                                         >  1:618310/1‑141 (MQ=255)
                                                                   tCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGCGCGCGGGAAAAATTAAACACGGGAGCCCCCCCACCAAAAAACGGGGGGGGCCCACCCAGGCGCGGAGCGatc                                                                                      >  1:1196059/1‑148 (MQ=255)
                                                                       aTCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGGAGCCCCCCCATCAGAAAACGGGGGGGGCCCACCCAGCCGCGGAGCGGTCATcccct                                                                               >  1:393395/1‑149 (MQ=255)
                                                                       aTCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCtt                                                                               <  1:524566/151‑1 (MQ=255)
                                                                       aTCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGCGCGGGAAAAATTAAACACGGGAGCCCCCCCCACAGAAAACGGGGTGGGCACACCCAGCCGCGTAGCGATCAAccc                                                                                 >  1:1035118/1‑149 (MQ=255)
                                                                         cAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGc                                                                             <  2:677719/151‑1 (MQ=255)
                                                                          aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACACGGGAGCCCCCCCATCAGAAAGCGTGGTGCGCACACCCAGCTGCGTAGCGATCATCCCTTTCt                                                                            >  2:760983/1‑151 (MQ=255)
                                                                              cTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACAACGGAGCCCCCCCATCAGAAAACGTGGGGGGCACACCCAGCCGCGGAGCGATCATCCCTTTCTGAt                                                                         >  1:1043559/1‑150 (MQ=255)
                                                                                 aTAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGGGGGAAAAATTAAACACGGGAGCCCCCCCAACAGAAAACGGGGGGGGGGCAACCAGCTGCGGAGCGAGCATCCCCTTCTGATaaa                                                                      >  1:38048/1‑150 (MQ=255)
                                                                                   aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCCCCTCAAAAAACGTGGTGGGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATaaaac                                                                    >  1:485623/1‑148 (MQ=255)
                                                                                   aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCa                                                                   <  1:737835/151‑1 (MQ=255)
                                                                                       tgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGGAAAAATTACACAACGGAGCCCCCCCACCAGAAAAGGGGGGGCGCACACCCAGGTGCGGAGCGATCAACCCTTGCTGATaaa                                                                      >  1:1171905/1‑144 (MQ=255)
                                                                                         tctcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAATTTAAACACGGGAGCCCCCCAATAAGAAAACGTGGTGGGCCCACCGAGCTGCGTAGCGATCACCCCTTGCTGTTAAAACa                                                                   >  2:108362/1‑145 (MQ=255)
                                                                                           tcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACACACGGGAGCCCCCCCATCAGAAAACGGGGGGCGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCACTgg                                                            >  2:1203341/1‑150 (MQ=255)
                                                                                            ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCCCATCAGAAAGCGGGGGGCGCACACCGAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCACCGGc                                                           >  2:340319/1‑150 (MQ=255)
                                                                                            ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTgg                                                            <  1:33938/149‑1 (MQ=255)
                                                                                            ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTg                                                             <  2:233266/148‑1 (MQ=255)
                                                                                                       gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCgg                                               <  2:309081/151‑1 (MQ=255)
                                                                                                       gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCgg                                               <  1:197749/151‑1 (MQ=255)
                                                                                                        agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACGGGAGCCCCCCCATCAGAAAACGTGGGGGGCACACCCAGCTGCGTAGCGATCATCCCTTTCTGATAAACCAACACTGGCACCACTTCGCggg                                              >  1:1112258/1‑151 (MQ=255)
                                                                                                          ggaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCCGGAAAAATTACACACCGGAGCCCCCCCATACGAAAGAGTGGGGGGCACACCCCGCCGCGGAGCGATCATCCCCTTCTGGTAAAACAACACTGGgacccccttccgggg                                             >  1:134578/1‑139 (MQ=255)
                                                                                                            aggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACAAACCGGAGCCCCCCCATCAGAAAACGGGGGGGGCACACCCGACTGCGTGGCGGTCATCCCCTGCTGATAAAACAGCACCGGCACCACCTCGCggg                                              >  1:592205/1‑147 (MQ=255)
                                                                                                              gaaagaaaAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtggt                                        <  1:1090975/151‑1 (MQ=255)
                                                                                                               aaagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATAAAACACGGTAGGCCCCCCATTAGAAAACATGGGGCGCACACCAAGATGTGGAGCGAATATGTCTTGCTGGTTAACGACAACCGGCACCCCCTCGCGGGTAtgggg                                       >  2:624908/1‑149 (MQ=255)
                                                                                                                aagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGGAACCCCCCCATCAGAAAACGGGGTGCGCACACCGAGCTGCGTAGCGATCATCCCCTGCTGATAAACCAACACTGGCACCCCCTCGCGGGGAtggggg                                      >  2:609851/1‑151 (MQ=255)
                                                                                                                 agaaaagaGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACGGTAGCCCCCACATCAGAAAACGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCACTGGCACCACTTCGCGGGGAtggggg                                      >  2:39926/1‑150 (MQ=255)
                                                                                                                             ccATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATgg                         <  2:276374/151‑1 (MQ=255)
                                                                                                                                 tttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCg                      <  2:1115736/150‑1 (MQ=255)
                                                                                                                                  ttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCgg                     <  2:894108/150‑1 (MQ=255)
                                                                                                                                   tgtgGCGGTGGGGCGCGGAAAAAATCACACACCGGAGCCCCCCCATCAGAAAACGGGGTGGGCACAACCAGCTGCGGAGCGATCATCCCTTGCTGATAAAACAGCACTGGCACCACCTCGCGGGGATGGTGGGTGTGACCAATGGTCgg                     >  1:663218/1‑149 (MQ=255)
                                                                                                                                    gtggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGTTGTGACCAATGGTCGGATc                  <  2:986414/151‑1 (MQ=255)
                                                                                                                                        cggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTg              <  1:156049/151‑1 (MQ=255)
                                                                                                                                            ggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAATCCAGCACTGGCAACACTTCGCGGGTATGGTGGCTGTGACCAATGGGCGGATCGTTGCCGt          <  1:777475/151‑1 (MQ=255)
                                                                                                                                            ggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCg           <  1:654602/150‑1 (MQ=255)
                                                                                                                                              cgcgcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTg         >  2:530265/1‑150 (MQ=255)
                                                                                                                                                   cgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg    >  2:1168696/1‑151 (MQ=255)
                                                                                                                                                     ggAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGc   <  2:29960/151‑1 (MQ=255)
                                                                                                                                                      gAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGcc  >  2:171749/1‑151 (MQ=255)
                                                                                                                                                       aaaaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGcc  >  2:555071/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
CAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCC  >  NC_000913/3508508‑3508806

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: