Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913993,9170AC46.2% 92.7 / 18.2 65V43G (GTG→GGG) ssuBaliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (31/4);  new base C (0/30);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 2.07e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTT  >  NC_000913/993768‑994058
                                                                                                                                                           |                                                                                                                                             
gTAGCAGACGCGCATCCTGGAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCccac                                                                                                                                              >  1:789752/1‑151 (MQ=255)
 tAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCccacc                                                                                                                                             >  2:1072015/1‑151 (MQ=255)
  aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCccacca                                                                                                                                            <  1:162730/151‑1 (MQ=255)
           gcATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCcaca                                                                                                                                    >  2:770971/1‑150 (MQ=255)
            cATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACaa                                                                                                                                   >  2:139623/1‑150 (MQ=255)
              tCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACt                                                                                                                                >  1:237866/1‑151 (MQ=255)
                      catcttCCGCGTTTTTTCTTGATTTCCACCGGGCGGTGTGGTCCCCGCTAACCCCTCCCCGGGGTTGGGCTTTCCCACACCTGCCAGCGGGCGCAGCAGGTTCTTTTCGCCCCCCCCGCTGGGCCCCACCACCGCCACAAACTGACCTg                                                                                                                          <  2:26134/149‑1 (MQ=255)
                                 tttCTTCCTAATTTTCAGCCAGGGGTGTGGTCCCCGTAAACTCCGCCCCTGGTTTGGGCTTTCCCAGCCCCGCCAGCAGGCGCAGCAGGTTTTTTTGCCCCCCCCCCCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATg                                                                                                                <  2:435109/146‑1 (MQ=255)
                                 tgtCTTCCGAATTTCCGCCCAGGGGTGTGGCCCCCCATACCCCATCCCCGGGGGTTGGCTTTCCCGACCCGCCAAGCGGGCGCAGCAGGGTCTTTTTACCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGt                                                                                                               <  2:468891/149‑1 (MQ=255)
                                      tcttaattttCAGCCAGGGGTGTGGTGCCCGCTAACACATCGCCTGGGGTTGGCGTTCCCAGCCCTGCCAGCAGGCGCAGCAGGGTACTTTTCCCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAAt                                                                                                           <  2:82900/141‑1 (MQ=255)
                                        cTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGt                                                                                                      >  2:106517/1‑151 (MQ=255)
                                         tGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGtt                                                                                                     >  2:484641/1‑151 (MQ=255)
                                          gAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg                                                                                                    >  2:649572/1‑151 (MQ=255)
                                           gaatcccccccagcGGTGGGGCCCCCCCAAACACCGCCCCTGGGGTGGGCTTTCCCAGCCCCGCCAGCGGCCGCAGAGGGTTATTTTCCCCACCCCCGCTGCCCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGt                                                                                                      <  1:384533/140‑1 (MQ=255)
                                           attttCAGCCAGGGGTGTGGTCCCCGCAAACACACCGCCTGGGTTTGGCTTTCCCAGCCCGCCCAGCGGGCGCAGCGGGGTATTTTTCCCCCCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg                                                                                                    <  2:1027950/148‑1 (MQ=255)
                                           aaTTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTgg                                                                                                   >  1:137410/1‑151 (MQ=255)
                                               tCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                >  1:639552/1‑150 (MQ=255)
                                                    ccgggggTGGGGTCCCCCCAACCACCCCGCCGGCGGTGGCCTTTCCCAGCCCTGCAAGCGGGCGCGGAGGGGTCTTTTGGCCCCCCCCGCGCCGCCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGa                                                                                            <  1:370412/145‑1 (MQ=255)
                                                     cgggggTGTGGTCCCCACTACCCCACCGCCGGCGTTGGGCTTTTCCAGACCTCCCGGCAGGCGCGGCGGGGTATTTTTCCCCCCCCCGCCGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg                                                                                                    <  2:401629/137‑1 (MQ=255)
                                                     cagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                          >  2:29282/1‑151 (MQ=255)
                                                         ggggtggggCCCCCCCACACACCGCCCTGGGGGTGGCTTTTCCCACCCCCGCCAGCAGGCCCACAGGTGTACTTTCCCCGCCCCCGCCCCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAgg                                                                                             <  1:198435/143‑1 (MQ=255)
                                                         ggTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                      >  2:752746/1‑151 (MQ=255)
                                                         ggTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                      >  1:154613/1‑151 (MQ=255)
                                                            gtggtgCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtttt                                                                                   >  2:268535/1‑151 (MQ=255)
                                                              ggtgCCCGCTAACCCCCCGCCTGGGTTTGGCTTTCCCAGCCCCCCCAGCGGCCGCAGCGGGTTACTTTCCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAttt                                                                                    <  2:359472/148‑1 (MQ=255)
                                                                  cccGCTAACACATCGCGTGGGTTGGGCTTTCCCAGCCCTGCCAGCAGGCGCAGCAGGGTTCTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgcg                                                                              <  2:249986/150‑1 (MQ=255)
                                                                  cccGATAACCCACCGCCGGGGGTGGGCTTTCCCAGCCCTGCCAGCGGGCGCAGCGGGTTTTTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgc                                                                               <  2:740326/149‑1 (MQ=255)
                                                                     gCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAAt                                                                          >  2:550134/1‑151 (MQ=255)
                                                                        aaCACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGtt                                                                       >  1:24072/1‑151 (MQ=255)
                                                                        aaCACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGtt                                                                       >  1:892392/1‑151 (MQ=255)
                                                                         accccccccccggCGTTTGGCTTTCCCAGCCCTCCCAGCGGGCGCAGCGGGTTATTTTCCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGt                                                                        <  1:47481/138‑1 (MQ=255)
                                                                         acacATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGttt                                                                      >  1:394783/1‑151 (MQ=255)
                                                                             aTCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                  >  1:1056043/1‑151 (MQ=255)
                                                                                   cTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc                                                             >  2:448926/1‑151 (MQ=255)
                                                                                     gCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAt                                                           >  2:671167/1‑151 (MQ=255)
                                                                                       gttgggCTTTCCCAGCCCGGCCAGCAGGCGCAGCAGGGTACTTTTCCCACCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                  <  2:623119/138‑1 (MQ=255)
                                                                                         tgggCTTTCCCAGCCCTGCCAGCGGGCGCGGAGGGGTTTTTTTCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAtt                                                          <  2:404587/146‑1 (MQ=255)
                                                                                           ggCTTTCCCAGCCCTGCCAGCGGGCGCAGCAGGGTTCTTTTCCCCCCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGca                                                     <  2:389433/151‑1 (MQ=255)
                                                                                                  cccggccccgCCAGGGGGCGCAGCGGGTTTCTTTCCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt                                                <  1:576271/141‑1 (MQ=255)
                                                                                                   ccagcccccccAGCGGGCGCAGAAGGGTTTTTTTCCCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGCTTGCTTACTGCATTGAGCAACAAt                                                <  1:619433/140‑1 (MQ=255)
                                                                                                   ccAGCCCGGCCAGCGGGGGCAGCGGGTTATTTGCGCCCCCCCGGCGGCCCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGAAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt                                                <  1:632896/149‑1 (MQ=255)
                                                                                                     aGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg                                            >  1:34891/1‑151 (MQ=255)
                                                                                                        ccgcccAGCGGGCGCGGCAGGTTCCTTTGCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                          <  2:140681/146‑1 (MQ=255)
                                                                                                         ccgccccgcggccgcAGCGGGGTCCTTTTCCCCCCCCCCCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGc                                         <  2:175596/139‑1 (MQ=255)
                                                                                                            ccAGCGGGCGAAGCAGGTTCCTTTCCCCCCCCCCGCGGCCCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCt                                      <  1:612804/150‑1 (MQ=255)
                                                                                                             ccggcgggCGCAGCGGGGTTTTTTTCCCCCCCCCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCt                                      <  1:253042/144‑1 (MQ=255)
                                                                                                               agcagGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGt                                   >  1:804279/1‑151 (MQ=255)
                                                                                                                gcagGCGCAGCGGGGTTTTTTTCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGt                                   <  2:249009/150‑1 (MQ=255)
                                                                                                                   ggcgcggaggggTATTTTTCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGt                                   <  1:845333/138‑1 (MQ=255)
                                                                                                                   ggCGCAGCAGGTTTCTTTTCCCCCCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCa                                <  2:1057044/150‑1 (MQ=255)
                                                                                                                       cagcagGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                            >  2:852072/1‑150 (MQ=255)
                                                                                                                       cagcagGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                            >  1:14652/1‑150 (MQ=255)
                                                                                                                        agcagGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                            >  2:753749/1‑149 (MQ=255)
                                                                                                                               gttattttGCCCACCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTa                     <  1:1091049/147‑1 (MQ=255)
                                                                                                                                     ttccccccccccgccgCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCAt                <  1:293568/140‑1 (MQ=255)
                                                                                                                                     ttcccccccaccGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAtt                   <  1:1083723/140‑1 (MQ=255)
                                                                                                                                      tgcccacccCCGCTCCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                  <  1:810143/144‑1 (MQ=255)
                                                                                                                                        cccccccccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCg            <  2:645248/144‑1 (MQ=255)
                                                                                                                                         ccccccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAcc             <  2:1046397/143‑1 (MQ=255)
                                                                                                                                           accaccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTgg         >  1:73289/1‑151 (MQ=255)
                                                                                                                                           accaccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTgg         >  2:507115/1‑151 (MQ=255)
                                                                                                                                             caccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcc       >  1:985356/1‑151 (MQ=255)
                                                                                                                                               ccgccgCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcct      <  1:211029/150‑1 (MQ=255)
                                                                                                                                               ccGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctc     <  1:385340/151‑1 (MQ=255)
                                                                                                                                                cGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcc    >  1:97424/1‑151 (MQ=255)
                                                                                                                                                 gCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcct   >  2:1048284/1‑151 (MQ=255)
                                                                                                                                                  cTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCtt  >  1:468583/1‑151 (MQ=255)
                                                                                                                                                  cTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCtt  >  1:762031/1‑151 (MQ=255)
                                                                                                                                                           |                                                                                                                                             
GTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTT  >  NC_000913/993768‑994058

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