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

AAGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGCTG  >  NC_000913/801970‑802261
                                                                                                                                                     |                                                                                                                                               
aaGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAac                                                                                                                                                <  2:1256918/150‑1 (MQ=255)
 aGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAacca                                                                                                                                              <  2:787066/151‑1 (MQ=255)
    gAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCAccg                                                                                                                                           >  2:593522/1‑151 (MQ=255)
      aCTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGATCCTGTCGTCCTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCAccgcc                                                                                                                                         <  2:213420/151‑1 (MQ=255)
              gggTGTTTTTTGGTATTTACCGGGCGGCTATTGTGGGGCTGTTTATCAAGCTTTTCCCGGAACCTGTCGTACTTTTAATTGCCGTTGCTGCCTCAATGGGGGTGGTGGGTAACTTACCCGACGGTGCTTTAAAACCCACCGCCGTATTaa                                                                                                                                  <  1:913768/150‑1 (MQ=255)
                           tttttACCGGGCGGCCTTTGTGGGGCGGTTTACCAAGCCTTCCCGGGACCCGCTCGTACTGTTAATTGCCGTTGCTCCCTAAATGGTGGTGGTCGGTAATTTATCCGCGGGTGCTTTAAAAACCACCGCCGTATTAAGCGGTTACTCtt                                                                                                                      <  1:615223/147‑1 (MQ=255)
                             tttACCGGGCGGCCTTTGTGGGGTTGGTTATAAACCTTTCCCCGGACCCTGTCGTACTGTTATTTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGCCGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAg                                                                                                                   <  2:1113327/150‑1 (MQ=255)
                                aCCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTac                                                                                                               >  1:176310/1‑151 (MQ=255)
                                 ccTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCTTATTGTTAATTGCCGTTGCTGCCTCAATGGGGGTGGTCGGTAACTTATCCGACGGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTa                                                                                                                <  2:732943/149‑1 (MQ=255)
                                 ccTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTacc                                                                                                              >  2:846935/1‑151 (MQ=255)
                                 ccTGCCGGCCATTGGGGGGCTGGTTATCAAGCCTTCCCCGGACCCTGTCGTCTTTTTAATTGCCGTTGCTGCCCCAAGGGTGGTGGTCGGTAACTTATCCGAGGGTGCGTTTAAACCAACCGCCGTATTAAGCGGTTACTCTTCAGGt                                                                                                                 <  2:134489/148‑1 (MQ=255)
                                  cTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTacc                                                                                                              >  1:1196132/1‑150 (MQ=255)
                                    ggcggcCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGACCCTGTCGTACTTTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTACCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTacca                                                                                                             <  2:134683/149‑1 (MQ=255)
                                        gcCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctggc                                                                                                       >  1:683011/1‑151 (MQ=255)
                                          cATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctggctg                                                                                                     >  1:980070/1‑151 (MQ=255)
                                            ttGTGGGTCTGTTTATCTGCCCTTCCCTGGAACCTGTCGTTTTGTAACTTGCCTGTGCTCCCTCAGTGTGGTGGGGGGTTATCTTCCCCGCCGGTCCGTTAAAGCCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctg                                                                                                         <  1:709491/145‑1 (MQ=255)
                                            ttGGGGGGCTGTTTATCAAGCCTTCCCGGGACCCTTTCTTACTTTTAATTGCCGTTGCTCCCTCAATGGGGGGGGTCGTTAATTTTTCCGACGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctggctg                                                                                                     <  1:704242/149‑1 (MQ=255)
                                              gtgGGGGGGTTTATCAACCTTTCCCCGGACCCTTTGTTACTTTTAATCGGCGTGCCTGCCTCAGGGTGGGGGGTCGTTACTTTATCCGTCGGTGCGTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctggct                                                                                                      <  1:624908/146‑1 (MQ=255)
                                                  ggCTGTTTATCAAGCCTTCCCGGGACCCTGTCTTACTTTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAATTTTCCCGACGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTtctc                                                                                             <  2:1150336/151‑1 (MQ=255)
                                                  ggCTGGTTATCAACCTTTTCCCGGACCCTGTCGTACTGTTAATTGCCGTTGCTGCCTAAATGGGGGTGGTCGGTAATTTATCCGACGGTGCGTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTtctc                                                                                             <  2:378303/151‑1 (MQ=255)
                                                     tGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGc                                                                                          >  1:855264/1‑151 (MQ=255)
                                                             acgcttttCCCGGACCCTGTCGTACTTTTAATTGCCGTTGCTGCCTAAATGGGGGTGGTCGTTAACTTTCCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTAcc                                                                                   <  1:119388/145‑1 (MQ=255)
                                                             aatccTTTCCCGGACCTTGTCGTCTTGTAAATTCCCGTTGCTGCCTAAAGGGTGGTGGTCGTTACGTGACCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTa                                                                                     <  2:1193019/145‑1 (MQ=255)
                                                             aacccTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTa                                                                                     <  2:1068441/145‑1 (MQ=255)
                                                             aaGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCt                                                                                  >  1:573010/1‑151 (MQ=255)
                                                             aaGCCTTCCCCGTACCCTGTCGTATTGTTTGTCCCGGTGGCTCCCAAAATGGGGGGGTTCGATACCTCAACCGAGGGTGTTATTAAACCCTCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGttt                                                                                      <  2:532587/147‑1 (MQ=255)
                                                                ccTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTa                                                                                >  2:271806/1‑150 (MQ=255)
                                                                  ttccccGGACCCTGTCGTACTTTTAATTGCCGTTGCTGCCTAAATGGTGGTGGTCGTTATCTTACCCGCGGGTGTTTTAAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTaa                                                                               <  2:473200/146‑1 (MQ=255)
                                                                       cgcaccgtgccGTACTCTAAATTGCGGTGGCTGCCTCAAGGGGGGTGGTGGTTAATTTGTCCGGCGGGGCTTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGcgccg                                                                         <  1:54983/140‑1 (MQ=255)
                                                                        ggAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGcgccgc                                                                        >  2:1285850/1‑150 (MQ=255)
                                                                          aaCCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTa                                                                                                                                   >  2:143908/1‑89 (MQ=255)
                                                                            ccTGCCGTACTTTTAATTGCCGTTGCTCCCTCAATGGTGGTGGTCGGTAACTTACCCGACGGTGCTTTAAAACCCACCGCCGTATTa                                                                                                                                   <  1:143906/87‑1 (MQ=255)
                                                                             cTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTtgt                                                                  >  2:23874/1‑151 (MQ=255)
                                                                                 cGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGacc                                                              >  2:280117/1‑151 (MQ=255)
                                                                                  gtaattgaaaTTGCCGTTCATGGCAAAAGGGTGTTGGTCGATATCTTCCCCTACGGTTTGTATAAACCCCCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGa                                                                <  1:1226472/140‑1 (MQ=255)
                                                                                      ttttAATGGCTGTTGCTGCCTAAATGGTGGTCGCCGTAAAATTCACGGCCGGGCCGTTTAAAACCTCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCg                                                          <  2:475117/148‑1 (MQ=255)
                                                                                         tAAGTGCCGTTGCTGCCTAAAGGGTGGTGCTGGGTATCTTACCCGACGGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGtt                                                       <  2:527093/150‑1 (MQ=255)
                                                                                          atttGCCGTTGCTGCCTCAAGGGTGGTGGTCGGTATTTTATCCGAGGGTGCGTTAAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGttt                                                      <  2:826921/148‑1 (MQ=255)
                                                                                          atttGCCGTTGCGCCCTAAATGGTGGTGGTCGGTAACTTATCCGACGGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGttt                                                      <  2:910700/148‑1 (MQ=255)
                                                                                          aaTTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGttt                                                      >  1:738548/1‑150 (MQ=255)
                                                                                          aaTTGCCGTTGCTGCCTCAATGGTGGAGGTCGGTAATTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGttt                                                      >  2:982667/1‑150 (MQ=255)
                                                                                          aaTTGCCGTTGCTCCCTAAATGGGGGTGGTCGGTAACTTACCCGACGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGttt                                                      <  2:789586/150‑1 (MQ=255)
                                                                                           aTTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGttt                                                      >  1:165809/1‑149 (MQ=255)
                                                                                             tGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGt                                                  >  1:662093/1‑151 (MQ=255)
                                                                                               cctttGCTGCCTCAATGGTGGTGGTCGGTAACTTACCCGACGGTGTTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAgg                                                   <  2:974489/145‑1 (MQ=255)
                                                                                               ccGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTa                                                 >  1:596217/1‑150 (MQ=255)
                                                                                                  ttGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACg                                             >  1:1086962/1‑151 (MQ=255)
                                                                                                  ttGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACg                                             >  2:918975/1‑151 (MQ=255)
                                                                                                   tgctgcCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACg                                             >  1:258788/1‑150 (MQ=255)
                                                                                                    gctccCTAAATGGTGGTGGTCGGTAATTTACCCGCGGGTGCTTTAAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACg                                             <  2:954348/149‑1 (MQ=255)
                                                                                                        ccgcaggggtggtggtGGGTAACTTACCCGGCGGTGTTTTAAAACCCGCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACg                                             <  1:1095562/138‑1 (MQ=255)
                                                                                                        ccTCAATGGTGGTGGCCGGTACCTTATCCGAGGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAtt                                         <  2:549386/149‑1 (MQ=255)
                                                                                                             aTGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATc                                  >  1:335732/1‑151 (MQ=255)
                                                                                                               ggtggtgttCTGTACTTTCTCGGCGGGTGTGTTCAAAACCTCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGTCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAt                                          <  2:919633/141‑1 (MQ=255)
                                                                                                                 tggtggtCGTAAATTTTCCCGTCGGGGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctg                                <  1:1298351/149‑1 (MQ=255)
                                                                                                                  ggtggtCGTTATCTTATCCGACGGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctgct                              <  2:68044/150‑1 (MQ=255)
                                                                                                                  ggtggtCGGTAACTTACCCGCGGGTGCGTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctgct                              <  2:125874/150‑1 (MQ=255)
                                                                                                                      gtgggTAACTTATCCGCCGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGAtt                          <  2:768753/147‑1 (MQ=255)
                                                                                                                         ggTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGt                       >  2:704452/1‑150 (MQ=255)
                                                                                                                         cgttactttACCCGAGGGTGGTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTg                         <  1:3876/142‑1 (MQ=255)
                                                                                                                          gTAACTAACCCGCCGGTGCTTTAAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGt                       <  1:1276798/149‑1 (MQ=255)
                                                                                                                             aaCTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTaaa                    >  2:845417/1‑150 (MQ=255)
                                                                                                                             aaCTTACCGGACGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTaa                     <  1:530223/149‑1 (MQ=255)
                                                                                                                                ttATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCg                >  2:380048/1‑151 (MQ=255)
                                                                                                                                 tATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCgg               >  1:922981/1‑151 (MQ=255)
                                                                                                                                  aTCCGACGGTGCGTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATc                 <  2:341989/148‑1 (MQ=255)
                                                                                                                                    ccGAGGGTGTGTTAAACCCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAAt                  <  1:695165/145‑1 (MQ=255)
                                                                                                                                     cGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAac           >  1:259664/1‑151 (MQ=255)
                                                                                                                                       aCGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAacac         >  1:535744/1‑151 (MQ=255)
                                                                                                                                        cGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAAcacc        >  2:299781/1‑151 (MQ=255)
                                                                                                                                        cGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAAcacc        >  1:424664/1‑151 (MQ=255)
                                                                                                                                         ggTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAaca          <  1:956775/148‑1 (MQ=255)
                                                                                                                                           tGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACg     >  2:885246/1‑151 (MQ=255)
                                                                                                                                            gCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGc    >  1:713567/1‑151 (MQ=255)
                                                                                                                                              gTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctg  >  2:214843/1‑151 (MQ=255)
                                                                                                                                              gTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctg  >  2:1277866/1‑151 (MQ=255)
                                                                                                                                               ttaaaacCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAAcaccac      <  1:577147/143‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
AAGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGCTG  >  NC_000913/801970‑802261

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