Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913535,9180AG42.5% 112.3 / 28.6 80E70G (GAG→GGG) hyihydroxypyruvate isomerase
Reads supporting (aligned to +/- strand):  ref base A (9/37);  new base G (34/0);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 2.18e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAG‑CGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACG  >  NC_000913/535771‑536067
                                                                                                                                                   |                                                                                                                                                        
ttGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAgc                                                                                                                                                        <  1:1165990/150‑1 (MQ=255)
 tGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAgcg                                                                                                                                                       <  2:1096564/150‑1 (MQ=255)
  gagaAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAgcgc                                                                                                                                                      <  2:580357/150‑1 (MQ=255)
   agaAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCgg                                                                                                                                                    <  2:1473627/151‑1 (MQ=255)
    gaAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTACACAATTTACCGGCGGGTGGATGGGGGGGGGGGGGGCgggg                                                                                                                                                    >  1:1426927/1‑150 (MQ=255)
      aaGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTA‑t                                                                                                                                                <  2:1092623/151‑1 (MQ=255)
       aGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGGGGGGGGGGGCggg                                                                                                                                                   >  1:572006/1‑147 (MQ=255)
          tGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGCGGGGGGGGGCggggat                                                                                                                                                >  2:619090/1‑144 (MQ=255)
          tGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGGGGGGGGGGGCgg                                                                                                                                                    >  2:847648/1‑144 (MQ=255)
          tGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGGGCGGGGGGGGGGCGCgggagt                                                                                                                                                >  2:146524/1‑144 (MQ=255)
          tGCGCAGTGTGGTTTCCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGCGCGGGA‑TTG‑cc                                                                                                                                           >  1:1263650/1‑151 (MQ=255)
           gcgcAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAAGTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCt                                                                                                                                          <  1:661661/151‑1 (MQ=255)
            cgcAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTg                                                                                                                                         <  1:751375/151‑1 (MQ=255)
             gcAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGt                                                                                                                                        <  1:378895/151‑1 (MQ=255)
                  gtgGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGAATGGGGGGGGGGGGGGGGCGGTG‑ttttccttattt                                                                                                                                     >  2:1438786/1‑140 (MQ=255)
                    ggTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTg                                                                                                                                 <  2:730448/151‑1 (MQ=255)
                    ggTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGGGGGGGGGCGCGGGA‑TTT‑CCTTTAtttc                                                                                                                                   >  1:594699/1‑147 (MQ=255)
                          cgcgGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGt                                                                                                                            <  2:1240927/150‑1 (MQ=255)
                            cggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTgaa                                                                                                                         <  2:892854/151‑1 (MQ=255)
                                  tGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGGCTGGGGGGGGGGGGAGCGCGGGA‑TTG‑CCTGTATTTTTGGGCGGggagaa                                                                                                                      >  1:1295479/1‑148 (MQ=255)
                                     aTTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGGGGGGGGGGGCGGGA‑TTT‑CCTGTATTTCCGGGCGGGGAGAAGAGgtt                                                                                                                >  2:92379/1‑151 (MQ=255)
                                      tttATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGGGCGGGGGGGGGGCGGGGGA‑GTG‑TCCGTATTTCTTGCCGTggaggaaggatt                                                                                                                >  2:262379/1‑142 (MQ=255)
                                       ttATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTc                                                                                                               <  2:1319119/150‑1 (MQ=255)
                                          tGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTAGGGGGGGGGGGAGGGCGGGA‑TTT‑TCTGTATTTCTTGCCGTGGAGAAAAGttat                                                                                                               >  2:402027/1‑145 (MQ=255)
                                            tttCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACCTGGGGGCGGGGGGGCGCCGGG‑TTA‑CCCGTGCTACTCGCCGGGGAGAAGAGgtt                                                                                                                >  1:749826/1‑144 (MQ=255)
                                            tttCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGAGCGCGGGA‑TTG‑CCTGTATTTCTGGCCGGGAAGAAGGGGTTCggggt                                                                                                          >  1:321808/1‑148 (MQ=255)
                                              tCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATgg                                                                                                        <  2:1178867/150‑1 (MQ=255)
                                                cTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGGCGGGGGGGGCGGGATTTT‑CCTG‑ATTTCTGTCCGGGGAGAAGGGGTTCTGGGTgggg                                                                                                      >  2:890101/1‑148 (MQ=255)
                                                  tATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTa                                                                                                    <  2:1019318/150‑1 (MQ=255)
                                                  tATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAGCGGGGAGCGGGGGGGGGCGGGT‑TTA‑CCTCTATTTCTGGGCGGGAAAAAGAGTTTTGGGGTGGCGt                                                                                                     >  2:1118750/1‑149 (MQ=255)
                                                   aTGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTTGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTag                                                                                                   <  2:1323090/150‑1 (MQ=255)
                                                   aTGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTagc                                                                                                  <  1:242439/151‑1 (MQ=255)
                                                     gACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTATGGGGTGGGGGGGCGCGGGA‑TTG‑TCTTGATTCCTTGGCGTTGAGAAAGGGTTCGGGGTGGGGTag                                                                                                   >  1:1296322/1‑148 (MQ=255)
                                                       cTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGGGGGGGGGCGCGGTA‑TTG‑CCTGTATTTCTGGCCGTGGAGAAGGGGTTTGGGGTGGCGTagca                                                                                                 >  1:1342729/1‑148 (MQ=255)
                                                       cTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTTGGCGGCGGGGGGGCGCGGTA‑TTT‑CCTGGTTTCCTGGCCCGGGAGAAGAGGTTCGGGATGGc                                                                                                       >  1:632706/1‑142 (MQ=255)
                                                         aCGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTagcagca                                                                                              <  2:642691/149‑1 (MQ=255)
                                                            aCATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGGGGGGGGGGGGGGCGGGA‑TTG‑CCTGTATTTCTGGCCGTGAAGAAGAGTTTTGGGGTGGCGTAGCAGCAG‑CGa                                                                                         >  1:1185151/1‑150 (MQ=255)
                                                            aCATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAG‑CGAt                                                                                        <  1:628581/151‑1 (MQ=255)
                                                              aTTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGGGGGGGGGCGCGGGA‑TTG‑CCTGTATTTCTTGCCGGGGAGAAGAGGTTTGGGGGGGGGGAGCAGCAG‑CGAtt                                                                                       >  2:972249/1‑150 (MQ=255)
                                                                tGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGCGCGGGA‑TTG‑CCTGTATTTCCGGCCGTGGAGAAGAGGTTTGGGGTGGGGGAGGAGCAG‑CGAtt                                                                                       >  1:1034680/1‑148 (MQ=255)
                                                                 gaagaaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAG‑CGATTCGt                                                                                    <  1:361512/150‑1 (MQ=255)
                                                                      aTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGGTGGGGGGGGGGGGGGGGCGGGA‑ATG‑GCTGTATTTCTGGCCGGGGAGAAGAGGTTTGGGGGGGGGTAGGAGCAG‑CGAGTCGtt                                                                                   >  2:393044/1‑146 (MQ=255)
                                                                      aTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGCGGGGGGGCGCGGGG‑TTGCCCTGTTTTCCTGG‑CGTGGAGAAGTGGTTCGGGGTGGCGGTGCAGCAG‑Cggattgtta                                                                                  >  2:501266/1‑139 (MQ=255)
                                                                       ttAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACCGGGCGGCGGGGGGGCGCGGGA‑TTG‑CCTGTATTTCTTGCCGTGGAGAAGAGTTTTGGGATGGGGGAGCAGCAG‑CGAATCGTTAtttcc                                                                             >  1:280097/1‑147 (MQ=255)
                                                                       ttAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGCGCCGTA‑TTG‑CCTGTATTTCTGGCCGTGAAGAAGAGGTTTGGGGTGGCGGAGCAGCCA‑GGATTCGTTATGccc                                                                             >  1:733727/1‑151 (MQ=255)
                                                                         aaaaCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAG‑CGATTCGTTATGCCCg                                                                            <  2:143403/150‑1 (MQ=255)
                                                                         aaaaCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAG‑CGATTCGTTATGCCCGt                                                                           <  1:1327278/151‑1 (MQ=255)
                                                                         aaaaCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAG‑CGATTCGTTATGCCCGt                                                                           <  1:612675/151‑1 (MQ=255)
                                                                           aaCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGGGGGGGGGGGGGGCGGGA‑TTT‑GCTGTATTCCTGGGCGGGGAGAAGAGGTTCGGGGTGGGGGAGCAGCAA‑CGATTCGTTATGCCCGTg                                                                          >  1:1237604/1‑150 (MQ=255)
                                                                           aaCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAG‑CGATTCGTTATGCCCGTgc                                                                         <  2:1155158/151‑1 (MQ=255)
                                                                            aCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGCGCGGGA‑TTT‑CCTGTATTTCTGGGCGTGAAGAAGAGGTTTGGGATGGGGGAGGAGGAG‑CGATTCGTTATTCCCGTgcg                                                                        >  1:1166643/1‑151 (MQ=255)
                                                                            aCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGGGCGGGA‑TTG‑CCTCTATTTCTGGCCGGGGAGAAGAGGTTCGGGGGGGGGGGGCAGGAG‑CGGTTCGGTATGCCCCTg                                                                          >  2:1442286/1‑149 (MQ=255)
                                                                              aTGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAG‑CGATTCGTTATGCCCGTGCGCt                                                                      <  1:232078/151‑1 (MQ=255)
                                                                               tgtgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGCGCGGGA‑TTG‑CCTGTATTTCTGGCCGTGGAGAAGAGTTTTGGGGTGGGGGAGCAGCAG‑CGATTCGTTATGCCCGTgcg                                                                        >  1:882105/1‑148 (MQ=255)
                                                                                gtgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAG‑CGATTCGTTATGCCCGTGCGCTgg                                                                    <  2:421255/151‑1 (MQ=255)
                                                                                gtgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAG‑CGATTCGTTATGCCCGTGCGCTgg                                                                    <  1:567145/151‑1 (MQ=255)
                                                                                     ggCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGGCTGGGCGGCGGGGGGGCGCGGGA‑TTG‑TCTGTTTTTCTGGCCGGGGAGAAGGGGTTTGGGGGGGGGGAGCAGGAA‑CGGTTTGTTAATCCCCTGCGCTgggga                                                                 >  1:237733/1‑147 (MQ=255)
                                                                                       cGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACCGGGGGGCGGGGGGGCGCGGGA‑TTG‑CCTGTATTTCCGGCCGGGAAGAAGAGTTTCGGGAGGGCGTAGAAGGAG‑CGATTTGTTATGCCCCGGCGCTggggcata                                                              >  1:950746/1‑145 (MQ=255)
                                                                                           taataaACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGGGCGGCGGGGGGGGGCGGGA‑TTG‑CCTGTATTCCTTGGCGGGGAGAAAAGGTTTGGGGGGGGGTCGCAGCCG‑CCAATTGTTTTTCCCCTGCGCTGGGGAATaaaaa                                                          >  2:1063355/1‑150 (MQ=255)
                                                                                            aataaACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGGGGGGCGGGGGGGGGCGGGA‑TTG‑CCTGTATTCCTGGCCGGGAAGAAGAGGTTCGGGATGGGGGAGCAGCAG‑CGATTCGTTATGGCCGGGCGCTGGGGAATaaaaaa                                                         >  1:1066410/1‑150 (MQ=255)
                                                                                               aaaCTCGAACATACGCTGCACAATTTACCGGCGGGGGACTCGGCGGCGGGGGGGGGCGGGA‑TTT‑GCCGTATTTCTTGCCGTGGAGAAAAATTTTGGGGTGGCGGAGGAGCCG‑CGGATTGTTTTTCCCGTGCGCCGGGTTATAAAAAAAtt                                                      >  1:752426/1‑150 (MQ=255)
                                                                                                aaCTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGCGCGGTA‑TTG‑CCTGTATTTCTCGCCCGGAAGAAGAGGTTTGGGGTGGGGGAGGAGCAG‑CGATTCGTTATGCCCGGGCGCTGGGGGATaaaaaaa                                                        >  2:256952/1‑147 (MQ=255)
                                                                                                     gAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAG‑CGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTc                                               <  1:920391/151‑1 (MQ=255)
                                                                                                     gAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAG‑CGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTc                                               <  1:619104/151‑1 (MQ=255)
                                                                                                        cATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAG‑CGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTgg                                            <  2:999912/151‑1 (MQ=255)
                                                                                                             gCTGCACAATTTACCGGCGGGTGAGTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAG‑TGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGt                                       <  2:1310260/151‑1 (MQ=255)
                                                                                                              cTGCACAATTTACCGGCGGGTGAATGGGGGGGGGGGGGGTGCGGTG‑TTG‑CCTGCCTTTCCGGGCGCGGAGGAAAGGTTTGGGGTGGGGGGGGAGGAGCCGATTGGTT‑TGTTCACGCGGGGGGGAATAAAAAAATTAAATGTCTGTTc                                          >  2:1187336/1‑147 (MQ=255)
                                                                                                                gCACAATTTACCGGCGGGTGACTGGACGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAG‑CGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTa                                      <  2:1299535/149‑1 (MQ=255)
                                                                                                                  acaATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAG‑CGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAAc                                  <  1:723133/151‑1 (MQ=255)
                                                                                                                       ttACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAG‑CGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCg                              <  1:1264979/150‑1 (MQ=255)
                                                                                                                            ggcggGTGACTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAG‑CGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGt                        >  1:469417/1‑151 (MQ=255)
                                                                                                                                  tGACTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAG‑CGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTcagc                  >  2:14083/1‑151 (MQ=255)
                                                                                                                                    aCTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAG‑CGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTcagc                  >  1:138093/1‑149 (MQ=255)
                                                                                                                                      tGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAG‑CGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTg              >  2:950149/1‑151 (MQ=255)
                                                                                                                                         gcggcgGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAG‑CGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAAc           <  1:14816/151‑1 (MQ=255)
                                                                                                                                             cggGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAG‑CGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGAt       <  2:1266004/151‑1 (MQ=255)
                                                                                                                                               ggggAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAG‑CGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGAtt      >  1:1261460/1‑150 (MQ=255)
                                                                                                                                               ggggAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAG‑CGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTc     <  1:1213192/151‑1 (MQ=255)
                                                                                                                                                  gAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAG‑CGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACg  >  1:1084670/1‑151 (MQ=255)
                                                                                                                                                  gAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAG‑CGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACg  <  2:986261/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                        
TTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTA‑TTG‑CCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAG‑CGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACG  >  NC_000913/535771‑536067

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