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

TTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACG  >  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)
      aaGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTAt                                                                                                                                                <  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)
          tGCGCAGTGTGGTTTCCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGCGCGGGATTGcc                                                                                                                                           >  1:1263650/1‑151 (MQ=255)
           gcgcAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAAGTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCt                                                                                                                                          <  1:661661/151‑1 (MQ=255)
            cgcAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTg                                                                                                                                         <  1:751375/151‑1 (MQ=255)
             gcAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGt                                                                                                                                        <  1:378895/151‑1 (MQ=255)
                  gtgGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGAATGGGGGGGGGGGGGGGGCGGTGttttccttattt                                                                                                                                     >  2:1438786/1‑140 (MQ=255)
                    ggTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTg                                                                                                                                 <  2:730448/151‑1 (MQ=255)
                    ggTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGGGGGGGGGCGCGGGATTTCCTTTAtttc                                                                                                                                   >  1:594699/1‑147 (MQ=255)
                          cgcgGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGt                                                                                                                            <  2:1240927/150‑1 (MQ=255)
                            cggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTgaa                                                                                                                         <  2:892854/151‑1 (MQ=255)
                                  tGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGGCTGGGGGGGGGGGGAGCGCGGGATTGCCTGTATTTTTGGGCGGggagaa                                                                                                                      >  1:1295479/1‑148 (MQ=255)
                                     aTTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGGGGGGGGGGGCGGGATTTCCTGTATTTCCGGGCGGGGAGAAGAGgtt                                                                                                                >  2:92379/1‑151 (MQ=255)
                                      tttATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGGGCGGGGGGGGGGCGGGGGAGTGTCCGTATTTCTTGCCGTggaggaaggatt                                                                                                                >  2:262379/1‑142 (MQ=255)
                                       ttATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTc                                                                                                               <  2:1319119/150‑1 (MQ=255)
                                          tGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTAGGGGGGGGGGGAGGGCGGGATTTTCTGTATTTCTTGCCGTGGAGAAAAGttat                                                                                                               >  2:402027/1‑145 (MQ=255)
                                            tttCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACCTGGGGGCGGGGGGGCGCCGGGTTACCCGTGCTACTCGCCGGGGAGAAGAGgtt                                                                                                                >  1:749826/1‑144 (MQ=255)
                                            tttCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGAGCGCGGGATTGCCTGTATTTCTGGCCGGGAAGAAGGGGTTCggggt                                                                                                          >  1:321808/1‑148 (MQ=255)
                                              tCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATgg                                                                                                        <  2:1178867/150‑1 (MQ=255)
                                                cTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGGCGGGGGGGGCGGGATTTTCCTGATTTCTGTCCGGGGAGAAGGGGTTCTGGGTgggg                                                                                                      >  2:890101/1‑148 (MQ=255)
                                                  tATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTa                                                                                                    <  2:1019318/150‑1 (MQ=255)
                                                  tATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAGCGGGGAGCGGGGGGGGGCGGGTTTACCTCTATTTCTGGGCGGGAAAAAGAGTTTTGGGGTGGCGt                                                                                                     >  2:1118750/1‑149 (MQ=255)
                                                   aTGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTTGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTag                                                                                                   <  2:1323090/150‑1 (MQ=255)
                                                   aTGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTagc                                                                                                  <  1:242439/151‑1 (MQ=255)
                                                     gACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTATGGGGTGGGGGGGCGCGGGATTGTCTTGATTCCTTGGCGTTGAGAAAGGGTTCGGGGTGGGGTag                                                                                                   >  1:1296322/1‑148 (MQ=255)
                                                       cTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGGGGGGGGGCGCGGTATTGCCTGTATTTCTGGCCGTGGAGAAGGGGTTTGGGGTGGCGTagca                                                                                                 >  1:1342729/1‑148 (MQ=255)
                                                       cTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTTGGCGGCGGGGGGGCGCGGTATTTCCTGGTTTCCTGGCCCGGGAGAAGAGGTTCGGGATGGc                                                                                                       >  1:632706/1‑142 (MQ=255)
                                                         aCGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTagcagca                                                                                              <  2:642691/149‑1 (MQ=255)
                                                            aCATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGGGGGGGGGGGGGGCGGGATTGCCTGTATTTCTGGCCGTGAAGAAGAGTTTTGGGGTGGCGTAGCAGCAGCGa                                                                                         >  1:1185151/1‑150 (MQ=255)
                                                            aCATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGAt                                                                                        <  1:628581/151‑1 (MQ=255)
                                                              aTTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGGGGGGGGGCGCGGGATTGCCTGTATTTCTTGCCGGGGAGAAGAGGTTTGGGGGGGGGGAGCAGCAGCGAtt                                                                                       >  2:972249/1‑150 (MQ=255)
                                                                tGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGCGCGGGATTGCCTGTATTTCCGGCCGTGGAGAAGAGGTTTGGGGTGGGGGAGGAGCAGCGAtt                                                                                       >  1:1034680/1‑148 (MQ=255)
                                                                 gaagaaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGt                                                                                    <  1:361512/150‑1 (MQ=255)
                                                                      aTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGGTGGGGGGGGGGGGGGGGCGGGAATGGCTGTATTTCTGGCCGGGGAGAAGAGGTTTGGGGGGGGGTAGGAGCAGCGAGTCGtt                                                                                   >  2:393044/1‑146 (MQ=255)
                                                                      aTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGCGGGGGGGCGCGGGGTTGCCCTGTTTTCCTGGCGTGGAGAAGTGGTTCGGGGTGGCGGTGCAGCAGCggattgtta                                                                                  >  2:501266/1‑139 (MQ=255)
                                                                       ttAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACCGGGCGGCGGGGGGGCGCGGGATTGCCTGTATTTCTTGCCGTGGAGAAGAGTTTTGGGATGGGGGAGCAGCAGCGAATCGTTAtttcc                                                                             >  1:280097/1‑147 (MQ=255)
                                                                       ttAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGCGCCGTATTGCCTGTATTTCTGGCCGTGAAGAAGAGGTTTGGGGTGGCGGAGCAGCCAGGATTCGTTATGccc                                                                             >  1:733727/1‑151 (MQ=255)
                                                                         aaaaCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCg                                                                            <  2:143403/150‑1 (MQ=255)
                                                                         aaaaCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGt                                                                           <  1:1327278/151‑1 (MQ=255)
                                                                         aaaaCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGt                                                                           <  1:612675/151‑1 (MQ=255)
                                                                           aaCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGGGGGGGGGGGGGGCGGGATTTGCTGTATTCCTGGGCGGGGAGAAGAGGTTCGGGGTGGGGGAGCAGCAACGATTCGTTATGCCCGTg                                                                          >  1:1237604/1‑150 (MQ=255)
                                                                           aaCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTgc                                                                         <  2:1155158/151‑1 (MQ=255)
                                                                            aCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGCGCGGGATTTCCTGTATTTCTGGGCGTGAAGAAGAGGTTTGGGATGGGGGAGGAGGAGCGATTCGTTATTCCCGTgcg                                                                        >  1:1166643/1‑151 (MQ=255)
                                                                            aCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGGGCGGGATTGCCTCTATTTCTGGCCGGGGAGAAGAGGTTCGGGGGGGGGGGGCAGGAGCGGTTCGGTATGCCCCTg                                                                          >  2:1442286/1‑149 (MQ=255)
                                                                              aTGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCt                                                                      <  1:232078/151‑1 (MQ=255)
                                                                               tgtgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGCGCGGGATTGCCTGTATTTCTGGCCGTGGAGAAGAGTTTTGGGGTGGGGGAGCAGCAGCGATTCGTTATGCCCGTgcg                                                                        >  1:882105/1‑148 (MQ=255)
                                                                                gtgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTgg                                                                    <  2:421255/151‑1 (MQ=255)
                                                                                gtgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTgg                                                                    <  1:567145/151‑1 (MQ=255)
                                                                                     ggCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGGCTGGGCGGCGGGGGGGCGCGGGATTGTCTGTTTTTCTGGCCGGGGAGAAGGGGTTTGGGGGGGGGGAGCAGGAACGGTTTGTTAATCCCCTGCGCTgggga                                                                 >  1:237733/1‑147 (MQ=255)
                                                                                       cGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACCGGGGGGCGGGGGGGCGCGGGATTGCCTGTATTTCCGGCCGGGAAGAAGAGTTTCGGGAGGGCGTAGAAGGAGCGATTTGTTATGCCCCGGCGCTggggcata                                                              >  1:950746/1‑145 (MQ=255)
                                                                                           taataaACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGGGCGGCGGGGGGGGGCGGGATTGCCTGTATTCCTTGGCGGGGAGAAAAGGTTTGGGGGGGGGTCGCAGCCGCCAATTGTTTTTCCCCTGCGCTGGGGAATaaaaa                                                          >  2:1063355/1‑150 (MQ=255)
                                                                                            aataaACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGGGGGGCGGGGGGGGGCGGGATTGCCTGTATTCCTGGCCGGGAAGAAGAGGTTCGGGATGGGGGAGCAGCAGCGATTCGTTATGGCCGGGCGCTGGGGAATaaaaaa                                                         >  1:1066410/1‑150 (MQ=255)
                                                                                               aaaCTCGAACATACGCTGCACAATTTACCGGCGGGGGACTCGGCGGCGGGGGGGGGCGGGATTTGCCGTATTTCTTGCCGTGGAGAAAAATTTTGGGGTGGCGGAGGAGCCGCGGATTGTTTTTCCCGTGCGCCGGGTTATAAAAAAAtt                                                      >  1:752426/1‑150 (MQ=255)
                                                                                                aaCTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGCGCGGTATTGCCTGTATTTCTCGCCCGGAAGAAGAGGTTTGGGGTGGGGGAGGAGCAGCGATTCGTTATGCCCGGGCGCTGGGGGATaaaaaaa                                                        >  2:256952/1‑147 (MQ=255)
                                                                                                     gAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTc                                               <  1:920391/151‑1 (MQ=255)
                                                                                                     gAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTc                                               <  1:619104/151‑1 (MQ=255)
                                                                                                        cATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTgg                                            <  2:999912/151‑1 (MQ=255)
                                                                                                             gCTGCACAATTTACCGGCGGGTGAGTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGTGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGt                                       <  2:1310260/151‑1 (MQ=255)
                                                                                                              cTGCACAATTTACCGGCGGGTGAATGGGGGGGGGGGGGGTGCGGTGTTGCCTGCCTTTCCGGGCGCGGAGGAAAGGTTTGGGGTGGGGGGGGAGGAGCCGATTGGTTTGTTCACGCGGGGGGGAATAAAAAAATTAAATGTCTGTTc                                          >  2:1187336/1‑147 (MQ=255)
                                                                                                                gCACAATTTACCGGCGGGTGACTGGACGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTa                                      <  2:1299535/149‑1 (MQ=255)
                                                                                                                  acaATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAAc                                  <  1:723133/151‑1 (MQ=255)
                                                                                                                       ttACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCg                              <  1:1264979/150‑1 (MQ=255)
                                                                                                                            ggcggGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGt                        >  1:469417/1‑151 (MQ=255)
                                                                                                                                  tGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTcagc                  >  2:14083/1‑151 (MQ=255)
                                                                                                                                    aCTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTcagc                  >  1:138093/1‑149 (MQ=255)
                                                                                                                                      tGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTg              >  2:950149/1‑151 (MQ=255)
                                                                                                                                         gcggcgGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAAc           <  1:14816/151‑1 (MQ=255)
                                                                                                                                             cggGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGAt       <  2:1266004/151‑1 (MQ=255)
                                                                                                                                               ggggAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGAtt      >  1:1261460/1‑150 (MQ=255)
                                                                                                                                               ggggAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTc     <  1:1213192/151‑1 (MQ=255)
                                                                                                                                                  gAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACg  >  1:1084670/1‑151 (MQ=255)
                                                                                                                                                  gAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACg  <  2:986261/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                        
TTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACG  >  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: