Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913802,1180AC46.0% 112.1 / 20.7 74T78P (ACC→CCC) ybhIputative DASS family tricarboxylate or dicarboxylate transporter
Reads supporting (aligned to +/- strand):  ref base A (31/9);  new base C (0/34);  total (31/43)
Fisher's exact test for biased strand distribution p-value = 7.39e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGCTGGGTCTGG  >  NC_000913/801970‑802268
                                                                                                                                                      |                                                                                                                                                       
aaGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAacc                                                                                                                                                      <  1:1240964/151‑1 (MQ=255)
   cGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAAACTGTCGTACTGTTAATTGCCGTTGCTGCCTCAAGGGTGGTGGTCGGTAACTTATCCGCCGGTGCGTTTAAAtccacc                                                                                                                                                   <  2:1265980/151‑1 (MQ=255)
     aaCTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCAccgc                                                                                                                                                 <  2:985393/151‑1 (MQ=255)
          gCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTAt                                                                                                                                            <  2:984004/151‑1 (MQ=255)
           cGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTAtt                                                                                                                                           >  1:1241601/1‑151 (MQ=255)
           cGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTAtt                                                                                                                                           >  1:1429644/1‑151 (MQ=255)
            gTGGGTGTTTTTTGGTATTTCCCTGGCGGCCATGGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTCCCTCAAGGGTGGTGGTCGGTAACTTACCCGACGGTGCGTTAAAAACCACCGCCGTAtt                                                                                                                                           <  2:1293515/150‑1 (MQ=255)
              gggTGTTTTTTGGTATTTACCTGGCGGCCTTTGTGGGGCTGGTTATCAAGCCTTCCCCGGAACCTGTCTTACTGTTAATTGCCGTTGCTCCCTCAATGGTGGTGGTCGGTAACTTACCCGCCGGTGCGTTAAAACCCACCGCCGTATTa                                                                                                                                          <  2:832444/149‑1 (MQ=255)
              gggTGTTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTCCCCGGACCCTGTCGTACTGTTAATTGCCGTGGCTGCCTCAATGGTGGTGGTCGGTAACTTACCCGACGGTGCGTTTAAACCCACCGCCGTAtt                                                                                                                                           <  1:1320135/148‑1 (MQ=255)
                  gCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGt                                                                                                                                    >  1:563129/1‑151 (MQ=255)
                      tttGTTTTTTCCCTGGCGGCCTTTGTGGGGCTGGTTATCAAGCTTTCCCGGGACCCTGTCGTACTGTTAATTCCCGTTGCTCCCACAATGGTGGTGTTGGGTAATTTATCCGAGGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTAc                                                                                                                                 <  2:129249/150‑1 (MQ=255)
                        tgttttttCCCTGGGGGCTTTTGTGGGGCTGGTTATCAAGCTTTCCCCGAACCTTTTCGTACTTTTATTTGCTTTTGCTCCCTCAATGGTGGTGGTCGGTAATTTATCCGACGGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTActc                                                                                                                               <  2:1317118/145‑1 (MQ=255)
                        tGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTActct                                                                                                                              >  1:1416271/1‑151 (MQ=255)
                        tGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTActct                                                                                                                              >  1:327782/1‑151 (MQ=255)
                        tGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTActct                                                                                                                              >  1:228012/1‑151 (MQ=255)
                          gTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTc                                                                                                                            >  1:774379/1‑151 (MQ=255)
                             tttCCCGGGCGCCCTTTGTGGGGTTGTTTACCAAGCTTTCCCCGGACCCGGTCTTACTGTAATTTGCCGTTGCTGCCCCAATGGGGGTGGTCGGTAATTTACCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCa                                                                                                                           <  2:813497/149‑1 (MQ=255)
                              ttACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGt                                                                                                                        >  2:778106/1‑151 (MQ=255)
                                  cTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTacc                                                                                                                     >  2:291440/1‑150 (MQ=255)
                                           aTTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctggctgg                                                                                                           >  2:6908/1‑151 (MQ=255)
                                             tgtgGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGgtg                                                                                                         >  2:510141/1‑151 (MQ=255)
                                                ggggTTGGTTATCAAGCCTTTCCCGGACCCTGTCGTATTGTTAATTGCCTTTGCTGCCTCAATGGTGGTGGTCGGTAACTTACCCGACGGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGtt                                                                                                       <  1:97203/150‑1 (MQ=255)
                                                  ggCTGGTTATAAAGCCTTCCCCGGACCCTGTCGTATTGTTAATTGCGGTTGCTGCCTCAATGGGGGTGGTGGGTAACTTACCCGGGGGGGTTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTtctc                                                                                                    <  2:526971/151‑1 (MQ=255)
                                                    cTGGTTATCAGGCCTTTCCCGGACCCTGTCGTACTGTTAATTGCCGTTGCTCCCTAAATGGGGGTGGTCGGTACCTTACCCGACGGTGCTTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCg                                                                                                   <  1:342589/150‑1 (MQ=255)
                                                      ggTTATCAAGGCTTTCCCGTAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGCTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCATCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCg                                                                                                >  2:607565/1‑151 (MQ=255)
                                                        ttATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGtt                                                                                              >  2:628947/1‑151 (MQ=255)
                                                             aaGCCTTCCCCGGACCCTGTCGTACTGTTATTGGCCGTTGCTCCCTCAATGGTGGTGGTCGGTAACTTACCCGACGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCt                                                                                         <  1:496187/151‑1 (MQ=255)
                                                                 cTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTa                                                                                       >  1:332493/1‑149 (MQ=255)
                                                                  tttCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAg                                                                                     >  2:913607/1‑150 (MQ=255)
                                                                   ttCCCGGACCCTGTCTTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgc                                                                                    <  1:145492/150‑1 (MQ=255)
                                                                   ttCCCGGAACCTGTCGTATTGTTAATGGCCGTTGCTGCCTAAAGGGGGGTGGTCGTTAATTTATCCGACGGTGCGTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgc                                                                                    <  2:1376120/150‑1 (MQ=255)
                                                                   tcccgggACCCTGTCTTACTGTTAATTGCCGTTGGTGCCTCAAGGGTGGGGGTGGGTAACTTACCCGAGGGTGCGTTAAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgc                                                                                    <  1:936046/145‑1 (MQ=255)
                                                                         gacccTGTCGTACTGTTAATGGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAACCCGCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCAt                                                                             <  1:6908/148‑1 (MQ=255)
                                                                          acccTTTCGTCCTTTAAATGGCCGTTGTTGCCTCAATGGTGGTGGTCGTAAATTTACCCGCCGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCAtt                                                                            <  1:420387/149‑1 (MQ=255)
                                                                          aaCCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCAtt                                                                            >  1:814956/1‑151 (MQ=255)
                                                                                gTCGTATTGTTATTTGCCGTTGCTGCCTAAAGGGGGGTGGTCGTTAACTTATCCGACGGTGCGTTAAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTtg                                                                          <  2:1110359/148‑1 (MQ=255)
                                                                                    tACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGaccac                                                                   >  1:1277914/1‑151 (MQ=255)
                                                                                     aCTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGaccacc                                                                  >  2:1146517/1‑151 (MQ=255)
                                                                                      cTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCg                                                                 >  1:1073185/1‑151 (MQ=255)
                                                                                         ttAATTGCCGTTGCTCCCTAAAGGGTGGTGGTCGGTAACTTACCCGACGGGGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGacca                                                                    <  2:670471/145‑1 (MQ=255)
                                                                                          tAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAATTTATCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGtt                                                              <  1:1079961/150‑1 (MQ=255)
                                                                                           atttGCCGTGGCTGCCTCAATGGTGGTGGTGGGTAACTTATCCGACGGTGCGTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTa                                                            <  1:731068/149‑1 (MQ=255)
                                                                                              ggcggtggctgcCCCTAGGGTGGTGGCGGTAAACTGCCCCGACGGTGCTTTAAACCCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAggc                                                         <  2:608940/144‑2 (MQ=255)
                                                                                                ccGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGCCGGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGt                                                         <  1:1315686/149‑1 (MQ=255)
                                                                                                ccGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTa                                                        >  1:269046/1‑150 (MQ=255)
                                                                                                    tgctgcCTCAGGGGGGGTGTTCGGTAATTTACCCGCCGGTGCTTTAAAACCCGCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTaaa                                                      <  2:904262/148‑1 (MQ=255)
                                                                                                    tgctgcCTCAATGGTGGTGGTCGGTAATTTACCCGGCGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAgg                                                          <  2:332488/144‑1 (MQ=255)
                                                                                                     gctgcCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGt                                                   >  2:1165987/1‑150 (MQ=255)
                                                                                                      ctgcCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACg                                                    <  1:778121/148‑1 (MQ=255)
                                                                                                      ctgcCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAt                                                 >  1:514492/1‑151 (MQ=255)
                                                                                                        gcCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTg                                               >  1:266436/1‑151 (MQ=255)
                                                                                                         ccTAAATGGGGGTGGTGGGTAACTTACCCGGCGGTGCTTTTAACCCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAtt                                                <  2:893895/149‑1 (MQ=255)
                                                                                                           tCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCt                                            >  1:534888/1‑151 (MQ=255)
                                                                                                             aaggggggtggtCGGTAACTTACCCGGCGGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAt                                                 <  1:1384230/138‑1 (MQ=255)
                                                                                                             aaTGGTGGTGGTGGTTAATTTATCCGCCGGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAtt                                                <  2:904962/145‑1 (MQ=255)
                                                                                                             aaTGGTGGTGGTCGTTATCTTACCCGACGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTg                                               <  2:1133747/146‑1 (MQ=255)
                                                                                                             aaTGGTGGTGGTCGGTAACTTATCCGAGGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAt                                                 <  2:1077662/144‑1 (MQ=255)
                                                                                                                ggtggtggtCGGTATTTTACCCGACGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATct                                        <  2:781380/150‑1 (MQ=255)
                                                                                                                ggtggtggtCGGTAACTTATCCGACGGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATct                                        <  2:1242754/150‑1 (MQ=255)
                                                                                                                ggtgggggtCGGTAACTTACCCGAGGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCt                                            <  2:1371710/146‑1 (MQ=255)
                                                                                                                      gtccGGTACTTTATCCGCCGGTCTTTTAAAACCCTCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATct                                        <  1:628955/141‑1 (MQ=255)
                                                                                                                         cGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGt                              >  1:1277862/1‑151 (MQ=255)
                                                                                                                         cGGTAACTTACCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTgg                               <  2:123517/150‑1 (MQ=255)
                                                                                                                                 ttaccggACGGTGCGTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATc                        <  2:1085651/144‑1 (MQ=255)
                                                                                                                                 ttATCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTaaaa                          <  2:1165082/148‑1 (MQ=255)
                                                                                                                                  tATCCGACGGGCCGTTAAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATc                        <  1:417950/149‑1 (MQ=255)
                                                                                                                                   aTCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGt                     >  1:1403167/1‑151 (MQ=255)
                                                                                                                                    tCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGt                     >  2:1045031/1‑150 (MQ=255)
                                                                                                                                    cccgagggcggttttAACCCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTa                    <  2:1145847/139‑1 (MQ=255)
                                                                                                                                     ccaagggTGCGTTAAAACCCTCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTa                    <  2:213558/145‑1 (MQ=255)
                                                                                                                                          ggTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAAcacc               <  1:204467/150‑1 (MQ=255)
                                                                                                                                          ggTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAAcacc               >  1:442154/1‑150 (MQ=255)
                                                                                                                                            tGCGTTTAAAACCACCGCTGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACg            >  1:1017530/1‑151 (MQ=255)
                                                                                                                                              cGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGct          >  1:377175/1‑151 (MQ=255)
                                                                                                                                              cGTTTAAAACCACCGCCATATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGct          >  1:84148/1‑151 (MQ=255)
                                                                                                                                               ttttAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACg            <  1:946516/147‑1 (MQ=255)
                                                                                                                                                   aaacCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggtc     <  2:220093/151‑1 (MQ=255)
                                                                                                                                                   aaaaCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggtc     >  2:1377332/1‑151 (MQ=255)
                                                                                                                                                      accaccGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggtctgg  >  2:1315503/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                       
AAGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGCTGGGTCTGG  >  NC_000913/801970‑802268

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