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

AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATG  >  NC_000913/993770‑994064
                                                                                                                                                             |                                                                                                                                                   
aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCccacca                                                                                                                                                  <  2:931716/151‑1 (MQ=255)
   aGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCg                                                                                                                                               >  2:521674/1‑151 (MQ=255)
       gcgcATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCcac                                                                                                                                           >  1:345549/1‑151 (MQ=255)
                gAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGAcc                                                                                                                                  >  2:602068/1‑151 (MQ=255)
                 aaaCATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCt                                                                                                                                 >  1:1009163/1‑151 (MQ=255)
                  accatcatTCGCGTTTCTCCCTAAATTTCAGCCGGGGGTGGGGCCCCCGCTAACACCCCGCGTGGGGTTGGCTTTCCGAGCCCTGCCAGCGGGCGCAGCGGGGTCTTTTTGCCCCCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTg                                                                                                                                <  2:599729/149‑1 (MQ=255)
                     atcatTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCg                                                                                                                             >  2:860251/1‑151 (MQ=255)
                         tccgggtgtCTCCTGGATTTCCGCCCGGGGGTGGGGCCCCCGCAACCACCTCGCTGGGGTTGGGCTTTCCCGGCCCGGCCAGCGGGCGCGGCAGGTTTTTTTCGCCACCCCCGGTGGGCCCCCCCACCGCCACAAACTGACCTGc                                                                                                                               <  2:1233351/140‑1 (MQ=255)
                         tccgcgTTTCTTCCTAATTTCCAGCCGGGGGTGGGGCCCCCGCAAACCCCCCGCCGGGGGTGGGCGTTCCCAGCCCTGCCAGCGGGCGCAGCGGGTTCTTTTCCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCg                                                                                                                             <  1:848368/145‑1 (MQ=255)
                             cgTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGt                                                                                                                     >  1:339459/1‑151 (MQ=255)
                               tttCTTCCTGATTTCCAGCCAGGGGTGTGGTCCCCGCTAAAACGCCGCCTGGGGTGGGCTTTTCCAGCCCTGCCGGCAGGCGCAGCAGGGTTTTTTGCCCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAata                                                                                                                        <  1:483015/144‑1 (MQ=255)
                                gtCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTaa                                                                                                                   >  2:789890/1‑150 (MQ=255)
                                 tcttctTGATTTCCAGCGGGGGGTGGGGTCCCCCCAAACACCCCGCCGGGGTTGGGCTTTCCCACCCCGCCCAGCGGGCGCAGCGGGTTCTTTTGCCCCCCCCCGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGt                                                                                                                     <  2:1075870/147‑1 (MQ=255)
                                   ttcgtaattttCGGCCAGGGGTTGGGTGCCCTTAACCACTCCGCCGGCGTTGGCCTTTCCCGGCCCGGCCAGAAGGCGAAGAAGGGTTCTTTTCCCCCCCCCGTGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTa                                                                                                                    <  1:499284/138‑1 (MQ=255)
                                   ttCTTGAATTTCACCCGGGGGTGGGGTCCCCGAAAACCCACCGCCTGCGGTGGCCTTTTCCACCCCTCCCGGCGGGCGCAGCAGGTTATTTTCCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGaa                                                                                                                          <  1:1258189/140‑1 (MQ=255)
                                     ctgaaaTTCCACCCGGGGGGGTGGCCCCCGCAAACCCACCGCCGGGGTTGGGCTTTCCCAGCCCTGCCAGCAGGCGCGGCGGGGTATTTTTCCCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAatat                                                                                                                       <  1:769547/137‑1 (MQ=255)
                                     cctaattctCAGCCGGGGGTGTGGCCCCCCCAAACCCGCCGCCTGGGGTGGCCTTTCCCACCCCGCCAAGCAGGCGCAGCAGGGTCTTTTTCCCCCCCCCGATGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCa                                                                                                              <  2:852230/142‑1 (MQ=255)
                                     ccTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                               >  1:966108/1‑149 (MQ=255)
                                      cTGAATTTCAGCCAGCGTTTTGTGCCCGGCTACCCCAACCCCGCCGGTTGCGGTCACGAGGCCAACCAGCCGGCCAACAGGGGTCTTTTTGCCCCCCCGGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTaaa                                                                                                                  <  1:256008/145‑1 (MQ=255)
                                          aTTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGt                                                                                                        >  2:188410/1‑151 (MQ=255)
                                           tttCACCCAGGGGTGTGGTGCCCGCTAACACATCGCCTGGGGTTGGCTTTCCCAGCCCTGCCAGCGGGCGCAGCAGGGTTTTTTTGCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGt                                                                                                        <  1:1086737/150‑1 (MQ=255)
                                            ttCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                      >  1:621418/1‑151 (MQ=255)
                                              cagccggcGGTGGGGGGCCCGTAACCACCCCCCCGGGGTTGGGGTTTTCCGGCCCGCCCGGCAGGCGCAGCGGGTTTTTTTGCCCCCCGCCGCGGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAg                                                                                                    <  1:865070/151‑1 (MQ=255)
                                                   cagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACTGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                                <  1:94501/151‑1 (MQ=255)
                                                     gggggggggcgCCCGCAAACACCCCCCCTGGGGTGGGCTTTCCCGGCCCGCCCAGCGGGCGCGGAAGTTTTTTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGat                                                                                              <  2:1114200/141‑1 (MQ=255)
                                                       ggTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                            >  2:300428/1‑151 (MQ=255)
                                                       ggTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                            >  1:606542/1‑151 (MQ=255)
                                                       ggTGGGGTCCCCGTAAACCCCCCGCCTGGGTTTGGTTTTCCCAGCCCTGCCAGGGGGCGCAGCAGGTTTTTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                                <  1:1178042/147‑1 (MQ=255)
                                                         tgtgGTCCCCGCTACCACACCGCCTGGGTTTGGTTTTCCCAGCCCGCCCAGCAGGCGCAGCGGGTTTTTTTTCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtt                                                                                           <  1:1174164/150‑1 (MQ=255)
                                                            ggcccccGATACCACACCGCCGGCGGTGGGCGTTCCCACCCCCCCCAGCGGGCGCAGCAGTGTTCTTTTCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTcc                                                                                       <  2:1123415/147‑1 (MQ=255)
                                                               gCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgcg                                                                                    >  1:1129388/1‑151 (MQ=255)
                                                                 ccGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGt                                                                                   >  1:627894/1‑150 (MQ=255)
                                                                     cTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATg                                                                               >  2:786649/1‑151 (MQ=255)
                                                                      tAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGt                                                                              >  1:1126745/1‑151 (MQ=255)
                                                                       aaccacgccgcgtgggTTGGCCTTTCCCAACCCCCCCACCAGGCGCAGCGGGTTCTTTTCCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTaa                                                                                 <  2:26508/134‑1 (MQ=255)
                                                                        accgcgccgcTGGGGGTGGCCTTTCCAAGCCTGCCCAGGGGGCGCAGCGGGTTACTTTGCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGt                                                                              <  2:136513/141‑1 (MQ=255)
                                                                          acaTCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGc                                                                          >  1:219852/1‑151 (MQ=255)
                                                                            acccccggggtttGGCTTTCCCGGCCCGGCCGGCGGGCGCAGGAGGGTATTTTTCCCCCCCCCGGTGGGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATg                                                                               <  1:371128/133‑1 (MQ=255)
                                                                              cGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTAc                                                                      >  1:671072/1‑151 (MQ=255)
                                                                              cGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTAc                                                                      >  2:668587/1‑151 (MQ=255)
                                                                              cGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTAc                                                                      >  1:875137/1‑151 (MQ=255)
                                                                               gCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACt                                                                     >  2:120911/1‑151 (MQ=255)
                                                                                 cTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc                                                                   >  1:926518/1‑151 (MQ=255)
                                                                                   ggggTGGGCTTTCCCACCCCGGCCGGCGGGCGCAGCGGGTTTTTTTCCCCCCCCCCGTTGCGCCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAt                                                                 <  2:663639/149‑1 (MQ=255)
                                                                                     ggTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                               >  1:262694/1‑151 (MQ=255)
                                                                                          ggCTTTCCCAGCCCTGCCAGCAGGCGCAGCAGGTTTCTTTCCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGc                                                            <  1:938761/150‑1 (MQ=255)
                                                                                           gCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGca                                                           >  2:658567/1‑150 (MQ=255)
                                                                                             gTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaac                                                         >  1:299262/1‑150 (MQ=255)
                                                                                               tttCCAGCCCCCCCAGCGGGCGCAGCGGGGTTCTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaa                                                          <  2:854437/148‑1 (MQ=255)
                                                                                                ttCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt                                                      >  2:651694/1‑151 (MQ=255)
                                                                                                   cccaccccgcccAGCGGGCGCAGGGGGTTTTTTTTCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATg                                                     <  2:929515/141‑1 (MQ=255)
                                                                                                    cAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGc                                                   >  2:768445/1‑151 (MQ=255)
                                                                                                     ccaccccgcccggcgggCGCAGCAGGGTTTTTTTGCCCCCCCCGCGGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATgg                                                    <  2:780292/136‑1 (MQ=255)
                                                                                                             ccAGCAGGCGCAGCGGGTTTCTTTCCCCCCCCCCGCGGCCCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCt                                            <  2:560489/150‑1 (MQ=255)
                                                                                                             ccAGCAGGCGCAGCAGGTTTCTTTCCCCCCCCCCGCGGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGccc                                             <  2:1126731/149‑1 (MQ=255)
                                                                                                                agcagGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGt                                         >  1:649962/1‑151 (MQ=255)
                                                                                                                  cagGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                       >  1:1267407/1‑151 (MQ=255)
                                                                                                                   agGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                       >  2:976294/1‑150 (MQ=255)
                                                                                                                    ggCGCAGCGGGTTTTTTTTCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                       <  1:1264362/149‑1 (MQ=255)
                                                                                                                         agcagGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                                  >  2:862594/1‑149 (MQ=255)
                                                                                                                                  aCTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                        >  2:170884/1‑151 (MQ=255)
                                                                                                                                   cTTTCGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGCTCAGACGAGCAGTATTc                        >  2:216185/1‑150 (MQ=255)
                                                                                                                                        tccccccccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATa                     <  1:352460/141‑1 (MQ=255)
                                                                                                                                          ccccccccccccgcgcCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCg                  <  2:621416/139‑1 (MQ=255)
                                                                                                                                           ccccccccGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGt                 <  1:305224/145‑1 (MQ=255)
                                                                                                                                           ccccccccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGt                 <  1:1118786/145‑1 (MQ=255)
                                                                                                                                           ccccccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                    <  1:105955/142‑1 (MQ=255)
                                                                                                                                              cccccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTgg               <  1:521676/147‑1 (MQ=255)
                                                                                                                                               cgcggttGCGCCCCCCCTCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCAt                      <  1:1265311/136‑1 (MQ=255)
                                                                                                                                                     tGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTc       >  1:1244670/1‑151 (MQ=255)
                                                                                                                                                          ccccccccCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCaaa    <  2:219851/149‑1 (MQ=255)
                                                                                                                                                          ccccccGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCaaa    <  1:941582/149‑1 (MQ=255)
                                                                                                                                                          ccccccACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTc       <  1:528612/146‑1 (MQ=255)
                                                                                                                                                          ccccccACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCaaa    <  1:897370/149‑1 (MQ=255)
                                                                                                                                                          cccACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATg  >  2:88377/1‑151 (MQ=255)
                                                                                                                                                             |                                                                                                                                                   
AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATG  >  NC_000913/993770‑994064

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: