Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913802,1180AC47.6% 61.1 / 17.9 63T78P (ACC→CCC) ybhIputative DASS family tricarboxylate or dicarboxylate transporter
Reads supporting (aligned to +/- strand):  ref base A (22/11);  new base C (0/30);  total (22/41)
Fisher's exact test for biased strand distribution p-value = 3.78e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGCTGGGTCTGG  >  NC_000913/801977‑802268
                                                                                                                                              |                                                                                                                                                       
cTGGCGTGGGGGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCACGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCAccgcc                                                                                                                                                >  2:758266/1‑150 (MQ=255)
   gCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTAt                                                                                                                                            >  1:135790/1‑151 (MQ=255)
      tGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTa                                                                                                                                          >  1:881928/1‑150 (MQ=255)
          tGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCgg                                                                                                                                     >  2:307805/1‑151 (MQ=255)
             tttttGTTATTTCCCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTCCCCGGACCCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAAGGGGGGTGGTCGGTAATTTACCCGACGGTGCGTTAAAACCCACCGCCGTATTAAGCGGtt                                                                                                                                   <  2:54790/150‑1 (MQ=255)
              tttgGGTATTTCCCTGGCGGCCTTTGTGGGGCGGGTTATCAACCTTTCCCCGAACCCGGTCTTATTGTTAATGGCCTTTGCTGCCTCAATGGTGGGGGTCGGTAATTTACCCGCCGGTGCTTTTAAACCCACCGCCGTATTAAGCGGtt                                                                                                                                   <  1:234232/149‑1 (MQ=255)
                 tggtttttccgggggggcAATTGTGGGCGTGTTAACCAAGCTTTTCCCGGAACTTTTCGTCTTTTTAATTCCCGTTGCTCCCTCAATGGTGTTGGTGGTTAACTTAGCCGACGGGGCGTTTAAACCCACCGCCGTATTAAGCGGTTActc                                                                                                                               <  2:917679/135‑1 (MQ=255)
                 ggttttttCCCGGGGGGCCTTTGGGGGGCTTGTTATAAGGCCTTCCCCGGACCTTTTCGTACTTTTAATTGCCGTTCCTCCCTAATGGGTGGTGGTGGGTACTTTTCCCGACGGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTActc                                                                                                                               <  2:1171811/145‑1 (MQ=255)
                    tATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTc                                                                                                                            >  1:115864/1‑150 (MQ=255)
                       ttcccgggcggcCATTGTGGGGCTGTTTAACAAGCTTTCCCCGGAACCTGTCGTACTGTTAATGGCCGTTGCTGCCTAAATGGGGGGGGTGGGTAACTTACCCGACGGTGCGTTAAAAACCACCGCCGTATTAAGCGGTTACTCTTCAgg                                                                                                                         <  1:155165/144‑1 (MQ=255)
                       ttACCTGGCGGCGTTTGGGGGGTTGTTTATCAAGCTTTTCCCGGACCCTTTGGTATTGTAAATTGCTGTTGCTGCCTAAAGGGTGGTGGTCGTTACCTTACCGGACGTTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAgg                                                                                                                         <  1:582837/150‑1 (MQ=255)
                          ccgggcggcCTTTGTGGGGCTGGTTATCAAGCTTTCCCCGGACCCTGTCGTACTGTTAATTGCCGTTGCTCCCTCAATGGTGGTGGTCGGTAATTTACCCGACGGGGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTa                                                                                                                       <  1:1116233/146‑1 (MQ=255)
                          ccTGGCGCCCTTTTGGGGGCTGTTTAACAAGCTTTTCCGGAACCCTGTCGTCTTGTTAATTGCGTTTGCTGCCCCAATGGTGGTGTTCGGTAATTTACCCGCCGGTGTTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGt                                                                                                                        <  1:574692/148‑1 (MQ=255)
                            tGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTaccac                                                                                                                   >  2:882617/1‑151 (MQ=255)
                             ggcggcCTTTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTTTCGTACTTTTAATTGCCGTTGCTGCCTAAAGGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTaccac                                                                                                                   <  1:1048852/150‑1 (MQ=255)
                                       gtgGGGCTGGTTATCAAGCTTTCCCCGGACCCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAATTTACCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGgt                                                                                                          <  1:256230/149‑1 (MQ=255)
                                       gtgGGGCTGGTTATAAACCTTTTCCCGGACCCTGTCTTATTTTTATTTGCCTTTGCTGCCTCAATGGTGGTGGTCGGTACCTTACCCGACGGTGCTTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGgtgt                                                                                                        <  1:996349/151‑1 (MQ=255)
                                        tgGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTCCCACCTGGCTGGTGtt                                                                                                       >  2:242876/1‑151 (MQ=255)
                                         ggggTTGGTTATCAAGCTTTCCCCGGAACTTGTCTTATTGTTAATTGCCGTTGCTGCCTAAATGGTGGTGGTCGGTAATTTATCCGACGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGtt                                                                                                       <  1:135459/150‑1 (MQ=255)
                                              tGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGACGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTtctc                                                                                                    >  2:795645/1‑148 (MQ=255)
                                                   aTCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGttt                                                                                             >  2:399578/1‑150 (MQ=255)
                                                      atgcttttCCGAGATCCTGTCGTATTGTTAATTACCTTCGCTGCCACAATGGTGGTGGTGGAAATCTTAGCCGTGGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTAcc                                                                                          <  2:268800/145‑1 (MQ=255)
                                                      aatccTTTCCCGAAACCGGCCGTTTTTTAAAGCGCCGTTGCGGCCTCATTGGGGTGGGTCGGTAACTTATCCGACGGTGCTTTAAAACCCTCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGttt                                                                                             <  1:117902/144‑1 (MQ=255)
                                                         ccTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTa                                                                                       >  2:133098/1‑150 (MQ=255)
                                                           ttccccGGAACCTTTCGTTCTTTAAATTGCCTTTGCTCCCTCAATGGTGGTGGTGGGAAACTTACCCGCCGGTGCTTTAAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCtt                                                                                        <  1:576432/144‑1 (MQ=255)
                                                             tCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgcgc                                                                                  >  1:385241/1‑151 (MQ=255)
                                                              cccGGACCTTGTCGTACTTTTAATTCCCTTTGCTCCCTAAAGGGTGGTGGTCGGTAACTTATCCGACGGGGTTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCtt                                                                                        <  2:219045/144‑1 (MQ=255)
                                                              cccGGACCCTGTGGTACTTTTATTTGCCTTTGCCGCCTCAAGGGTGGTGGTGGGTATTTTATCGGCCGGTGTTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgc                                                                                    <  2:1058433/148‑1 (MQ=255)
                                                               ccGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGcgccg                                                                                >  2:584112/1‑151 (MQ=255)
                                                                    acccTTTCGTACTGTTATTGGCCTTTGCTGCCACAATGGGGGTGGTCGGTATCTTATCCGAGGGTCTTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCa                                                                              <  1:539995/147‑1 (MQ=255)
                                                                       cTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTtgt                                                                         >  1:849749/1‑151 (MQ=255)
                                                                            gTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGacc                                                                     >  1:501278/1‑150 (MQ=255)
                                                                             tACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGaccac                                                                   >  1:726645/1‑151 (MQ=255)
                                                                              actttaatttGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAATTAACCCGCCGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGaccac                                                                   <  1:509861/142‑1 (MQ=255)
                                                                                tGTAAATGGCCGTTGCGCCCTAAAGGGGGGGGTGCGGTAACTTTCCCCACGTTGCTTTAAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGacca                                                                    <  2:236949/147‑1 (MQ=255)
                                                                                         ccGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTa                                                        >  2:968003/1‑150 (MQ=255)
                                                                                         ccGTTGCTGCCTCAAGGGTGGGGGTCGTTAACTTACCCGACGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGt                                                         <  1:290163/149‑1 (MQ=255)
                                                                                            ttGCTCCCTCAAGGGGGTTGGTCGGTAATTTATCCGGGGGTGCTTTAAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAAc                                                     <  2:583930/150‑1 (MQ=255)
                                                                                             tgctgcCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGt                                                   >  1:241880/1‑151 (MQ=255)
                                                                                                tgcCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAtt                                                >  1:474396/1‑151 (MQ=255)
                                                                                                 gcCTAAAGGGGGGTGTTGGTTAACTTATCCGACGGTGCGTTAAAAACCTCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGTCCACCGGTTTAGGTAAAc                                                     <  1:627812/145‑1 (MQ=255)
                                                                                                  ccTAAATGGTGGTGGTGGGTACTTTCTCCGAGGGGGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAtt                                                <  2:387393/149‑1 (MQ=255)
                                                                                                   ctaattggtgggggTGGTTAACTTACCCGAGGGTGCGTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGt                                                   <  1:897521/140‑1 (MQ=255)
                                                                                                      aaTGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTAt                                          >  2:627420/1‑151 (MQ=255)
                                                                                                      aaTGGTGGTGGTCGGTAACTTACCCGACGGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTa                                           <  1:201588/150‑1 (MQ=255)
                                                                                                         ggtggtggtCGGTATCTTACCCGACGGTGCGTTAAAACCCGCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATct                                        <  1:882677/150‑1 (MQ=255)
                                                                                                         gggggtggtCGTTATCTTACCCGACGGTGGTTTAAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctg                                       <  2:1154830/148‑1 (MQ=255)
                                                                                                               ggttggtttcTTACCCGAGGGTGGTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctgc                                      <  2:1105178/137‑1 (MQ=255)
                                                                                                                gtgggtattttACCCGGCGGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTg                                <  1:584140/142‑1 (MQ=255)
                                                                                                                  cGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGt                              >  2:862028/1‑151 (MQ=255)
                                                                                                                      aaCTTTCCCGTCGGTGCTTAAAAACCCCCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGAt                                  <  2:1177818/143‑1 (MQ=255)
                                                                                                                         ttACCCGACGGTGTTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGt                              <  1:339756/144‑1 (MQ=255)
                                                                                                                          ttacccgcgggtgttttAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATc                        <  2:209215/136‑1 (MQ=255)
                                                                                                                           acccGACGGTGCGTTAAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCgg                      <  1:806405/148‑1 (MQ=255)
                                                                                                                           aTCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCgg                      >  1:114505/1‑150 (MQ=255)
                                                                                                                             ccgacgtgtcGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTa                    <  2:137733/141‑1 (MQ=255)
                                                                                                                             ccGAGGGTGTTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTaa                   <  2:1169137/151‑1 (MQ=255)
                                                                                                                             ccGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGt                     <  2:455227/149‑1 (MQ=255)
                                                                                                                             ccGACGGTGCGTTTAAAACCGCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTa                    <  2:778510/150‑1 (MQ=255)
                                                                                                                                aCGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAacac                >  2:1139694/1‑151 (MQ=255)
                                                                                                                                  ggTGCGTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAAcacca              <  2:447237/151‑1 (MQ=255)
                                                                                                                                     gctttgaaacCCCCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTaa                   <  1:793139/137‑1 (MQ=255)
                                                                                                                                      gctttAAAACCCTCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACg            <  2:1168071/148‑1 (MQ=255)
                                                                                                                                        tttAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACg            <  2:201313/147‑1 (MQ=255)
                                                                                                                                          aaaacCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggt      <  2:683213/150‑1 (MQ=255)
                                                                                                                                           aaaaCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggtc     >  1:394129/1‑151 (MQ=255)
                                                                                                                                              accaccGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggtctgg  >  1:1055368/1‑151 (MQ=255)
                                                                                                                                              |                                                                                                                                                       
CTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGCTGGGTCTGG  >  NC_000913/801977‑802268

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