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

ACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGCTGGGTC  >  NC_000913/801976‑802265
                                                                                                                                                 |                                                                                                                                                   
aCTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCAccgcc                                                                                                                                             <  2:1023157/151‑1 (MQ=255)
  tGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGt                                                                                                                                           >  2:348783/1‑151 (MQ=255)
     cGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTAtt                                                                                                                                        >  1:498658/1‑151 (MQ=255)
      gTGGGTGTTTTTTGTTTTTTCCCTGGGGGCCATTGTGGGGCTGGTTATCAAGCTTTTCCCGGACCCTGTCTTACTTTTAATTGCCGTTGCTGCCTAAATGGTGGTGGTCGGTAACTTACCCGACGGTGCTTTTAAACCCACCGCCGTAtt                                                                                                                                        <  1:1005678/150‑1 (MQ=255)
       tGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTa                                                                                                                                       >  1:111629/1‑150 (MQ=255)
        gggTGTTTTTTGGTATTTCCCTGGCGGCCTTTGTGGGGCTGTTTATCAACCTTTTCCCGGACCCTTTCGTACTGTTAATTGCCGTTGCTGCCTCAAGGGTGGTGGTCGGTAATTTACCCGACGGTGCTTTTAAACCCACCGCCGTATTaa                                                                                                                                      <  1:634898/150‑1 (MQ=255)
              tttttGGTATTTCCCTGGGGGCCTTTGGGGGGCTGGTTACCAAGCCTTTCCCGGACCCTGTCGTATTGTTAATTGCCGTTGCTCCCTCAAGGGTGGTGGTCGGTACTTTACCCGACGGTGCTTTAAAAACCACCGCCGTATTAAGCgg                                                                                                                                  <  1:712758/148‑1 (MQ=255)
                tttGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTAc                                                                                                                              >  1:437734/1‑150 (MQ=255)
                   ggTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCtt                                                                                                                          >  1:700543/1‑151 (MQ=255)
                     tATTTCCCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTCCCCGGACCCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTACCCGAGGGTGCTTTTAAAACCACCGCCGTATTAAGCGGTTActct                                                                                                                           <  2:347718/148‑1 (MQ=255)
                       tttCCCGGGCGCCCATTGTGGGGCTGGTTATCAACCTTTCCCGGGACCCTGTCGTCCTGTTAATTGCCGTTGCTGCCTCAAGGGTGGTGGTCGGTAACTTACCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCtt                                                                                                                          <  2:135290/147‑1 (MQ=255)
                       tttACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAgg                                                                                                                      >  2:20026/1‑151 (MQ=255)
                       tttACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAgg                                                                                                                      >  2:550163/1‑151 (MQ=255)
                              ggcggcCATTGGGGGGCTGGTTATCAAGCCTTCCCCGGACCCTGTCGTTCTTTTAATTCCCGTTGCTGCCTCAAGGGTGGTGGTCGGTATCTTACCCGACGGTGCGTTAAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTaccac                                                                                                                <  2:237463/150‑1 (MQ=255)
                              ggcgccCATTGTGGGGCTGTTTACAAAGCCTTCCCCGGACCCTGTCTTACTTTTAATTGCCGTTGCTCCCTCAATGGGGGTGGCGGGTAATTTACCCGGCGGTGCGTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTacca                                                                                                                 <  1:535165/149‑1 (MQ=255)
                               gcggcattgggggggCGGTTTATCAAGCTTTCCCCGGACCCTGTGGTCCTTTTAATTGCCTTTGCTGCCTAAAGGGGGGTGGTGGGTACCTTACCCGACGTTGCGTTTAAACCCCCCGCCGTATTAAGCGGTTACTCTTCAGGTacca                                                                                                                 <  2:1011723/137‑1 (MQ=255)
                                  gcCATTGGGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTACCCGACGGTGCGTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctgg                                                                                                            <  1:542777/150‑1 (MQ=255)
                                                ggTTATCAAGCCTTCCCCGGAACCTGTCGTACTTTTAATTGCGGTTGCTGCCTCAATGGGGGTGGTCGGTATTTTATCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGc                                                                                              <  1:867202/150‑1 (MQ=255)
                                                 gTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGt                                                                                            >  2:430889/1‑151 (MQ=255)
                                                  ttATCAACCTTTTCCCGGACCCTGTCTTACTTTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCg                                                                                             <  1:940158/149‑1 (MQ=255)
                                                     tCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTAc                                                                                        >  1:759536/1‑151 (MQ=255)
                                                       aacccTTTCCGGGACCCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAAGGGTGGTGGTCGTTAACTTATCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTAc                                                                                        <  2:606060/146‑1 (MQ=255)
                                                        aaccTTTCCCCGGAACCTGTTTTACTTTTATTTGCCTTTGCTGCCTCAATGGTGGTGGTCGGTAATTTATCCGACGGTGCGTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCt                                                                                      <  2:440385/149‑1 (MQ=255)
                                                         gCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCtt                                                                                     >  1:512539/1‑150 (MQ=255)
                                                           cTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTa                                                                                    >  2:708732/1‑149 (MQ=255)
                                                            tttCCGGGACCCTGTCTTCTTGTTATTTGCCGTTGCTGCCTCAAGGGTGGTGGTCGTTAACTTATCCGACGTTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTaa                                                                                   <  2:1161791/149‑1 (MQ=255)
                                                             ttCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgcg                                                                                >  2:185845/1‑151 (MQ=255)
                                                             ttCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgcg                                                                                >  1:759937/1‑151 (MQ=255)
                                                             ttCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgcg                                                                                >  2:943575/1‑151 (MQ=255)
                                                                cccGGCTCCTGTCGTTTTGTTAATGCCGGTTGCTTCCTCAATGGGGGGGATAGTAACTTTACCCGGTTGGGCGTTTAAACCCGCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgct                                                                                <  1:1056129/149‑2 (MQ=255)
                                                                cccGGACCCTGTCGTCCTTTTAATTGCCGTTGCTCCCTCAATGGTGGTGGTCGGTAATTTACCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgcgc                                                                               <  1:558967/150‑1 (MQ=255)
                                                                cccGGACCCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTATTTTACCCGACGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgc                                                                                 <  1:1082680/148‑1 (MQ=255)
                                                                cccGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAAGGGGGGTGGTCGGTAATTTACCCGAGGTTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgcgc                                                                               <  2:579547/150‑1 (MQ=255)
                                                                      acccTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAATTTATCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCAt                                                                          <  1:335666/148‑1 (MQ=255)
                                                                      aaCCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCAtt                                                                         >  1:246054/1‑151 (MQ=255)
                                                                      aaCCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCAtt                                                                         >  2:738897/1‑151 (MQ=255)
                                                                           gTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTACCCGACGGGGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTtgt                                                                      <  2:1199890/149‑1 (MQ=255)
                                                                             cGTTCTTTTAATGGCCGTTGCTGCCTCACTGGGGGTGTTCGTTATTTTATCGGCCGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTtgt                                                                      <  1:703172/147‑1 (MQ=255)
                                                                                    ttATTTCCCTTTGCTCCCCAAATGGTGTTGGTCGTTATCTTACCCGACGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGaccac                                                                <  1:358781/146‑1 (MQ=255)
                                                                                    ttAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGacca                                                                 <  1:965518/145‑1 (MQ=255)
                                                                                    ttAATTCCCGTGGCTGCCTCAAGGGTGGTGGTCGGTAACTTTTCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCgg                                                             <  2:543468/149‑1 (MQ=255)
                                                                                    ttAATGGCCGTTGCTGCCTAAAGGGTGGTGGTCGGTAACTTACCCGACGGTGCTTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCgg                                                             <  2:104842/149‑1 (MQ=255)
                                                                                     tAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGttt                                                          >  2:36085/1‑151 (MQ=255)
                                                                                              tttccgcCTAAAGGGTGGTGGTCGTAAACTTATCCGACGTTGTTTTAAACACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTaaa                                                   <  1:995596/144‑1 (MQ=255)
                                                                                              ttgcgcccTAAAGGGTGGTGGTCGGTAACTTACCCGAGGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAg                                                        <  1:1182832/138‑1 (MQ=255)
                                                                                                gctgcCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGt                                                <  2:1164327/150‑1 (MQ=255)
                                                                                                   gcCTCAATGGTGGGGGTGGGTAACTTATCCGGCGGTGCTTTTAAACCCCCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAt                                              <  1:916011/149‑1 (MQ=255)
                                                                                                     cTCATTGGGGGTGGTCGTTACTTTACCCGAGGGTGCTTTAAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAtt                                             <  1:260905/148‑1 (MQ=255)
                                                                                                        agtggtggtggTTGGAAATTAATCCGTCGGGGCTTTAAAAACCTCCGCCGTATTAAGCGGTTACGCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTAt                                       <  1:730225/149‑1 (MQ=255)
                                                                                                        aaTGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTAt                                       >  2:785828/1‑151 (MQ=255)
                                                                                                        aaTGGGGGTGGTGGGTCATTTACCCGGCGGTGTTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTg                                            <  1:833494/146‑1 (MQ=255)
                                                                                                        aaTGGGGGTGGTCGGTAACTTACCCGACGGTGCGTTAAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTAt                                       <  1:550137/151‑1 (MQ=255)
                                                                                                           ggtggtggtTGGTACCTTATCCGGCGGTGCGTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTa                                        <  1:711773/147‑1 (MQ=255)
                                                                                                           ggtggtggtGGGTAACTTACCCGACGGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCt                                         <  2:613376/146‑1 (MQ=255)
                                                                                                           ggggggggtCGTTACCTTATCCGACGGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATct                                     <  2:1150321/144‑1 (MQ=255)
                                                                                                            gtggtggtCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctgc                                   >  1:1947/1‑151 (MQ=255)
                                                                                                              ggtggtGGGTACCTTATCCGACGGTCCGTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctgc                                   <  2:842658/149‑1 (MQ=255)
                                                                                                              ggtggtCGGTAATTTACCCGACGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctgc                                   <  1:904237/149‑1 (MQ=255)
                                                                                                              ggtggtCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctgctg                                 >  1:748229/1‑151 (MQ=255)
                                                                                                                  ggcGGTAACTTATCCGACGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGAt                               <  1:148952/147‑1 (MQ=255)
                                                                                                                   tCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTgg                            >  2:194728/1‑151 (MQ=255)
                                                                                                                    cGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGt                           >  1:387733/1‑151 (MQ=255)
                                                                                                                     ggTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGt                           >  2:127116/1‑150 (MQ=255)
                                                                                                                       tAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTaa                         >  2:163227/1‑150 (MQ=255)
                                                                                                                         acctttcccGACGGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTaaa                        <  1:147218/143‑1 (MQ=255)
                                                                                                                         aCTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAAt                      >  1:24915/1‑151 (MQ=255)
                                                                                                                              acccGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGt                  <  1:908525/149‑1 (MQ=255)
                                                                                                                               tCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGt                  >  2:546166/1‑150 (MQ=255)
                                                                                                                                ccgccgGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGt                  <  2:608744/149‑1 (MQ=255)
                                                                                                                                ccGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTa                 <  1:1060150/150‑1 (MQ=255)
                                                                                                                                ccGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTACGCGCCGCATTTGTGACCACCGGTTTAGGTAAAGGTATTGCCTATCTGCTGATTGGGAAAATCGGTa                 >  2:629790/1‑150 (MQ=255)
                                                                                                                                     ggTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAacac             <  2:111846/149‑1 (MQ=255)
                                                                                                                                     ggTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAAcacc            <  2:191937/150‑1 (MQ=255)
                                                                                                                                       tGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACg         >  1:981098/1‑151 (MQ=255)
                                                                                                                                        gCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAAcaccac          <  1:696834/149‑1 (MQ=255)
                                                                                                                                           tttAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACg         <  2:512563/147‑1 (MQ=255)
                                                                                                                                           tttAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgg     >  2:549470/1‑151 (MQ=255)
                                                                                                                                           ttaaaaaCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGct       <  2:673995/146‑1 (MQ=255)
                                                                                                                                              aaaaCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggtc  >  1:703821/1‑151 (MQ=255)
                                                                                                                                              aaaaCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggtc  >  2:243826/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGCTGGGTC  >  NC_000913/801976‑802265

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: