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

AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCC‑GCTAACACATCGCCTGCGGTT‑GGCG‑TTT‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATG  >  NC_000913/993770‑994064
                                                                                                                                                             |                                                                                                                                                   
aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCC‑GCTAACACATCGCCTGCGGTT‑GGCG‑TTT‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCccacca                                                                                                                                                  <  2:931716/151‑1 (MQ=255)
   aGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCC‑GCTAACACATCGCCTGCGGTT‑GGCG‑TTT‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCg                                                                                                                                               >  2:521674/1‑151 (MQ=255)
       gcgcATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCC‑GCTAACACATCGCCTGCGGTT‑GGCG‑TTT‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCcac                                                                                                                                           >  1:345549/1‑151 (MQ=255)
                gAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCC‑GCTAACACATCGCCTGCGGTT‑GGCG‑TTT‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGAcc                                                                                                                                  >  2:602068/1‑151 (MQ=255)
                 aaaCATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCC‑GCTAACACATCGCCTGCGGTT‑GGCG‑TTT‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCt                                                                                                                                 >  1:1009163/1‑151 (MQ=255)
                  accatcatTCGCGTTTCTCCCTAAATTTCAG‑CCGGGGGTGGGGCCCCC‑GCTAACACCCCGCGTGGGGTT‑GGCT‑TTC‑CGAGCCCTG‑CC‑AGCGGGCGCAGCG‑GGGTCTTTTT‑G‑CCCCCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTg                                                                                                                                <  2:599729/149‑1 (MQ=255)
                     atcatTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCC‑GCTAACACATCGCCTGCGGTT‑GGCG‑TTT‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCg                                                                                                                             >  2:860251/1‑151 (MQ=255)
                         tccgggtgtC‑TCCTGGATTTCCGCCCGGGGGTGGGGCCCCC‑GCAACCACCTCG‑CTGGGGTTGGGCT‑TTC‑CCGGCCCGG‑CC‑AGCGGGCGCGGCA‑GGTTTTTTTC‑G‑CCACCCCCGGTGGGCCCCCCCACCGCCACAAACTGACCTGc                                                                                                                               <  2:1233351/140‑1 (MQ=255)
                         tccgcgTTTCTTCCTAATTTCCAG‑CCGGGGGTGGGGCCCCC‑GCAAACCCCCCGCCGGGGGTG‑GGCG‑TTC‑CCAGCCCTG‑CC‑AGCGGGCGCAGC‑‑GGGTTCTTTTCC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCg                                                                                                                             <  1:848368/145‑1 (MQ=255)
                             cgTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCC‑GCTAACACATCGCCTGCGGTT‑GGCG‑TTT‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGt                                                                                                                     >  1:339459/1‑151 (MQ=255)
                               tttCTTCCTGATTTCCAG‑CCAGGGGTGTGGTCCCC‑GCTAAAACGCCGCCTGGGGTG‑GGCT‑TTT‑CCAGCCCTG‑CC‑GGCAGGCGCAGCA‑GGGTTTTTTG‑C‑CCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAata                                                                                                                        <  1:483015/144‑1 (MQ=255)
                                gtCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCC‑GCTAACACATCGCCTGCGGTT‑GGCG‑TTT‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTaa                                                                                                                   >  2:789890/1‑150 (MQ=255)
                                 tcttctTGATTTCCAG‑CGGGGGGTGGGGTCCCC‑CCAAACACCCCGCCGGGGTTG‑GGCT‑TTC‑CCACCCCGC‑CC‑AGCGGGCGCAGC‑‑GGGTTCTTTT‑GCCCCCCCCCGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGt                                                                                                                     <  2:1075870/147‑1 (MQ=255)
                                   ttcgtaattttCGG‑CCAGGGGTTGGGTGCCC‑TTAACCACTCCGCCGGC‑GTT‑GGCC‑TTTCCCGGCCCGG‑CC‑AGAAGGCGAAGAA‑GGGTTCTTTT‑C‑CCCCCCCCGTGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTa                                                                                                                    <  1:499284/138‑1 (MQ=255)
                                   ttCTTGAATTTCAC‑CCGGGGGTGGGGTCCCC‑GAAAACCCACCGCCTGCGGTG‑GCCT‑TTT‑CCACCCCTC‑CC‑GGCGGGCGCAGCA‑GGTTATTTTC‑C‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGaa                                                                                                                          <  1:1258189/140‑1 (MQ=255)
                                     ctgaaaTTCCAC‑CCGGGGGGGTGGCCCCC‑GCAAACCCACCGCCGGGGTTG‑GGCT‑TTC‑CCAGCCCTG‑CC‑AGCAGGCGCGGCG‑GGGTATTTTT‑C‑CCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAatat                                                                                                                       <  1:769547/137‑1 (MQ=255)
                                     cctaattctCAG‑CCGGGGGTGTGGCCCCC‑CCAAACCCGCCGCCTG‑GGGT‑GGCC‑TTT‑CCCACCCCG‑CCAAGCAGGCGCAGCA‑GGGTCTTTTT‑C‑CCCCCCCCGATGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCa                                                                                                              <  2:852230/142‑1 (MQ=255)
                                     ccTGAATTTCAG‑CCAGCGGTGTGGTGCCC‑GCTAACACATCGCCTGCGGTT‑GGCG‑TTT‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                               >  1:966108/1‑149 (MQ=255)
                                      cTGAATTTCAG‑CCAGC‑GTTTTGTGCCCGGCTACCCCAACCCCGCCGGTT‑GCGG‑TCA‑CGAGGCCAA‑CC‑AGCCGGC‑CAACAGGGGTCTTTTT‑G‑CCCCCCCGGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTaaa                                                                                                                  <  1:256008/145‑1 (MQ=255)
                                          aTTTCAG‑CCAGCGGTGTGGTGCCC‑GCTAACACATCGCCTGCGGTT‑GGCG‑TTT‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGt                                                                                                        >  2:188410/1‑151 (MQ=255)
                                           tttCAC‑CCAGGGGTGTGGTGCCC‑GCTAACACATCGCCTGGGGTT‑GGCT‑TTC‑CCAGCCCTG‑CC‑AGCGGGCGCAGCA‑GGGTTTTTTT‑G‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGt                                                                                                        <  1:1086737/150‑1 (MQ=255)
                                            ttCAG‑CCAGCGGTGTGGTGCCC‑GCTAACACATCGCCTGCGGTT‑GGCG‑TTT‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                      >  1:621418/1‑151 (MQ=255)
                                              cag‑ccggcGGTGGGGGGCCC‑G‑TAACCACCCCCCCGGGGTT‑GGGGTTTT‑CCGGCCCGC‑CC‑GGCAGGCGCAGC‑‑GGGTTTTTTT‑GCCCCCCGCCGCGGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAg                                                                                                    <  1:865070/151‑1 (MQ=255)
                                                   cagcGGTGTGGTGCCC‑GCTAACACATCGCCTGCGGTT‑GGCG‑TTT‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACTGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                                <  1:94501/151‑1 (MQ=255)
                                                     gggggggggcgCCC‑GCAAACACCCCCCCTGGGGTG‑GGCT‑TTC‑CCGGCCCGC‑CC‑AGCGGGCGCGGAA‑GTTTTTTTTT‑C‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGat                                                                                              <  2:1114200/141‑1 (MQ=255)
                                                       ggTGTGGTGCCC‑GCTAACACATCGCCTGCGGTT‑GGCG‑TTT‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                            >  2:300428/1‑151 (MQ=255)
                                                       ggTGTGGTGCCC‑GCTAACACATCGCCTGCGGTT‑GGCG‑TTT‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                            >  1:606542/1‑151 (MQ=255)
                                                       ggTGGGGTCCCC‑GTAAACCCCCCGCCTGGGTTT‑GGTT‑TTC‑CCAGCCCTG‑CC‑AGGGGGCGCAGCA‑GGTTTTTTTT‑C‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                                <  1:1178042/147‑1 (MQ=255)
                                                         tgtgGTCCCC‑GCTACCACACCGCCTGGGTTT‑GGTT‑TTC‑CCAGCCCGC‑CC‑AGCAGGCGCAGCG‑GGTTTTTTTT‑C‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtt                                                                                           <  1:1174164/150‑1 (MQ=255)
                                                            ggccccc‑GATACCACACCGCCGGCGGTG‑GGCG‑TTC‑CCACCCCCC‑CC‑AGCGGGCGCAGCA‑GTGTTCTTTT‑C‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTcc                                                                                       <  2:1123415/147‑1 (MQ=255)
                                                               gCCC‑GCTAACACATCGCCTGCGGTT‑GGCG‑TTT‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgcg                                                                                    >  1:1129388/1‑151 (MQ=255)
                                                                 cc‑GCTAACACATCGCCTGCGGTT‑GGCG‑TTT‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGt                                                                                   >  1:627894/1‑150 (MQ=255)
                                                                     cTAACACATCGCCTGCGGTT‑GGCG‑TTT‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATg                                                                               >  2:786649/1‑151 (MQ=255)
                                                                      tAACACATCGCCTGCGGTT‑GGCG‑TTT‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGt                                                                              >  1:1126745/1‑151 (MQ=255)
                                                                       aaccacgccgcgtgggTT‑GGCC‑TTT‑CCCAACCCCCCC‑ACCAGGCGCAGC‑‑GGGTTCTTTTCC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTaa                                                                                 <  2:26508/134‑1 (MQ=255)
                                                                        accgcgccgcTGGGGGT‑GGCC‑TTT‑CCAAGCCTGCCC‑AGGGGGCGCAGCG‑GGTTACTTTG‑C‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGt                                                                              <  2:136513/141‑1 (MQ=255)
                                                                          acaTCGCCTGCGGTT‑GGCG‑TTT‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGc                                                                          >  1:219852/1‑151 (MQ=255)
                                                                            acccccggggttt‑GGCT‑TTC‑CCGGCCCGG‑CC‑GGCGGGCGCAGGA‑GGGTATTTTT‑C‑CCCCCCCCGGTGGGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATg                                                                               <  1:371128/133‑1 (MQ=255)
                                                                              cGCCTGCGGTT‑GGCG‑TTT‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTAc                                                                      >  1:671072/1‑151 (MQ=255)
                                                                              cGCCTGCGGTT‑GGCG‑TTT‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTAc                                                                      >  2:668587/1‑151 (MQ=255)
                                                                              cGCCTGCGGTT‑GGCG‑TTT‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTAc                                                                      >  1:875137/1‑151 (MQ=255)
                                                                               gCCTGCGGTT‑GGCG‑TTT‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACt                                                                     >  2:120911/1‑151 (MQ=255)
                                                                                 cTGCGGTT‑GGCG‑TTT‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc                                                                   >  1:926518/1‑151 (MQ=255)
                                                                                   ggggTG‑GGCT‑TTC‑CCACCCCGG‑CC‑GGCGGGCGCAGCG‑GGTTTTTTTC‑C‑CCCCCCCCGTTGCGCCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAt                                                                 <  2:663639/149‑1 (MQ=255)
                                                                                     ggTT‑GGCG‑TTT‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                               >  1:262694/1‑151 (MQ=255)
                                                                                          ggCT‑TTC‑CCAGCCCTG‑CC‑AGCAGGCGCAGCA‑GGTTTCTTTC‑C‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGc                                                            <  1:938761/150‑1 (MQ=255)
                                                                                           gCG‑TTT‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGca                                                           >  2:658567/1‑150 (MQ=255)
                                                                                             g‑TTT‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaac                                                         >  1:299262/1‑150 (MQ=255)
                                                                                               ttt‑CCAGCCCCC‑CC‑AGCGGGCGCAGCG‑GGGTTCTTTT‑C‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaa                                                          <  2:854437/148‑1 (MQ=255)
                                                                                                tt‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt                                                      >  2:651694/1‑151 (MQ=255)
                                                                                                   cccaccccgccc‑AGCGGGCGCAGGG‑GGTTTTTTTT‑C‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATg                                                     <  2:929515/141‑1 (MQ=255)
                                                                                                    cAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGc                                                   >  2:768445/1‑151 (MQ=255)
                                                                                                     ccaccccgcccggcgggCGCAGCA‑GGGTTTTTTT‑G‑CCCCCCCCGCGGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATgg                                                    <  2:780292/136‑1 (MQ=255)
                                                                                                             cc‑AGCAGGCGCAGCG‑GGTTTCTTTC‑C‑CCCCCCCCGCGGCCCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCt                                            <  2:560489/150‑1 (MQ=255)
                                                                                                             cc‑AGCAGGCGCAGCA‑GGTTTCTTTC‑C‑CCCCCCCCGCGGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGccc                                             <  2:1126731/149‑1 (MQ=255)
                                                                                                                agcagGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGt                                         >  1:649962/1‑151 (MQ=255)
                                                                                                                  cagGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                       >  1:1267407/1‑151 (MQ=255)
                                                                                                                   agGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                       >  2:976294/1‑150 (MQ=255)
                                                                                                                    ggCGCAGCG‑GGTTTTTTTT‑C‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                       <  1:1264362/149‑1 (MQ=255)
                                                                                                                         agca‑gGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                                  >  2:862594/1‑149 (MQ=255)
                                                                                                                                  aCTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                        >  2:170884/1‑151 (MQ=255)
                                                                                                                                   cTTTC‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGCTCAGACGAGCAGTATTc                        >  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)
                                                                                                                                                             |                                                                                                                                                   
AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCC‑GCTAACACATCGCCTGCGGTT‑GGCG‑TTT‑CCAGACCTG‑CC‑AGCAGGCGCAGCA‑GGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATG  >  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: ‑