Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,506,0280AC45.0% 83.2 / 18.0 61V18G (GTA→GGA) ydcOBenE family inner membrane putative transporter
Reads supporting (aligned to +/- strand):  ref base A (31/2);  new base C (0/27);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 4.61e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACGCGATACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCG‑ACAATCGCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCA  >  NC_000913/1505880‑1506172
                                                                                                                                                       |                                                                                                                                                
tACGCGATACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCG‑ACAATCGCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATT‑AATAc                                                                                                                                                 <  1:22377/150‑1 (MQ=255)
      aTACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCG‑ACAATCGCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATT‑AATACCGCCAGa                                                                                                                                          >  1:722361/1‑151 (MQ=255)
        aCCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCG‑ACAATCGCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATT‑AATACCGCCAGaaa                                                                                                                                        >  1:408398/1‑151 (MQ=255)
          cATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCG‑ACAATCGCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATT‑AATACCGCCAGAAAc                                                                                                                                       >  2:1036919/1‑150 (MQ=255)
              aTGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCG‑ACAATCGCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATT‑ACTACCGCCAGAAACCCCg                                                                                                                                   >  2:940943/1‑150 (MQ=255)
              aTGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCG‑ACAATCGCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATT‑AATACCGCCAGAAACCCCg                                                                                                                                   >  1:1153784/1‑150 (MQ=255)
                        cAGCGTACTGACGACCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCG‑ACAATCGCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATT‑AATACCGCCAGAAACCCCGCCAGCAGCgtg                                                                                                                        >  2:846890/1‑151 (MQ=255)
                          gCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCG‑ACAATCGCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCgtggg                                                                                                                      >  1:194671/1‑151 (MQ=255)
                             tACTGCCGCCAATTGCCACCCCCGGCGCCGCCACCCAGCAAGAGTTTTGTGCGGTGGTGGCCCCGGCA‑ACAATCGCGGTTTGCAAGATTATGGCCGCGGAAC‑GGGGGTGGCAAATA‑AATCCCGCCAAAACCCCCGCCAGTAGCgtgggtgg                                                                                                                   <  2:1128226/151‑1 (MQ=255)
                                 gACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCG‑ACAATCGCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTggagga                                                                                                               >  2:453806/1‑151 (MQ=255)
                                  ccccccATTCCCACCCCCAGCCCCTTCACCCAGCCAGAGTTTTGTGCGGTGGGGGCCCCGGCG‑ACAACCGCTGCTTCCCAGATTTTGGCCGCGGAAC‑GGGCGTAGCAAATT‑AATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTggagga                                                                                                               <  1:1132015/147‑1 (MQ=255)
                                       cATTGCGAGCCCACCCGCCTTCGTCGAGGCAAAGATTGGTGGAGTGGTGGCTCCGGCG‑ACACCCGCGCGTTGCAAGATTCTTGCCGGTGAAG‑GGGCGAGGCAAATA‑AACCCCGCCGGAAACACCGCCAGTAGCGTGGGTGGAGGAATAGaa                                                                                                         <  2:565411/151‑1 (MQ=255)
                                        attcccccccccGGCGCCGTCCCCCGCCCAGAGTTTTGTGCGGGGGTGGCCCGGGCG‑AAAATCGCTGCTTGCCAGATTATTCCCGCTGAAC‑TGGGGTAGCAAATA‑AATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGaaa                                                                                                        <  1:656194/143‑1 (MQ=255)
                                           gCCAGCCCCAGCGCGCTCGTCCAGCCAGAGTTTGGGGCTG‑GGTGGCCCCGGAGTAGAGTCGGCGCTGCCCAGTTTTCCGCCGCGGACG‑GGGCGGAACAAATA‑AATCCCACCAAAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAAc                                                                                                      <  2:904784/150‑1 (MQ=255)
                                           gCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCG‑ACAATCGCTGCTTGCCCGATTATTGCCGCTGAAC‑TGGCGTAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACa                                                                                                     >  2:965589/1‑151 (MQ=255)
                                               cccacggcgcCGACTCCCGGCCGGAGATTGGTGGGTGGGTGGGTCGGGCG‑ACAACCGCTGCTGGCCAGATTAGTGCCGTTGTAC‑GGGGGGAGAAAATA‑AAGCCCCCCAAAACCACCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACg                                                                                                 <  2:246345/145‑1 (MQ=255)
                                                ccccAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCG‑ACAATCGCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACg                                                                                                 >  2:947557/1‑150 (MQ=255)
                                                  ccggcgcCGCCACCCAGCCAGAGTTTTGTGCAGGGGTGGCCCGGGCG‑ACAATCGCTGCTTGCCAGATTATGGCCGCTGAAC‑GGGGGTAGCAAATT‑AATCCCCCCAAAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCAt                                                                                              <  1:56417/148‑1 (MQ=255)
                                                      ccccGTCACCCACCAAGAGTTTTGGGCAGGGGGGGCTCCGGCA‑ACAACCGCTGTTTCCAAAATTATGGCCGCTAAAC‑GGGCGTAGCAAATT‑AATCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTc                                                                                          <  1:941812/149‑1 (MQ=255)
                                                          gTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCG‑ACAATCGCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCAc                                                                                        >  1:968588/1‑149 (MQ=255)
                                                          ctCATCCGGCGAAGGATTTGGGGGGTGTTGGCCCCGGCC‑ACATTGGCGTCTCCCCAGGTTTTTCCCGCGGAAT‑TGGCG‑AACAAATTAAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGAGGCATAGTCACCt                                                                                      <  1:848992/150‑1 (MQ=255)
                                                           ccatcccgcagaagaTTGGGGCAGTGGTGGCCCGGGAA‑ACAATCCCTCTTTGCAGGATTATCGCCG‑TGAACGGGGCAAAAACAATA‑AACCCCCCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCt                                                                                      <  2:867997/138‑1 (MQ=255)
                                                                 agccagAGATTTGTGCAGTGGTGGCTCCGGCG‑ACAATCGCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGt                                                                               >  2:732642/1‑151 (MQ=255)
                                                                       agttttGGGCAGGGGGGGCTCCGGCA‑AAAATCGCTGCTTCCAAGATTATCGCCGATGAAC‑TGGCGTAGCAAATT‑AACCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAg                                                                          <  2:495163/147‑1 (MQ=255)
                                                                       agtttggtgCGGTGGGGGCCCCGGCA‑ACAATCGCTGCTTGCCAGATTTTGCCCGCTAAAG‑GGGGGTAGCAAATT‑AATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAg                                                                          <  1:567427/144‑1 (MQ=255)
                                                                        gaTTTGTGCAGTGGTGGCTCCGGCG‑ACAATCGCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCg                                                                        >  2:352678/1‑151 (MQ=255)
                                                                        gaTTTGTGCAGTGGTGGCTCCGGCG‑ACAATCGCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCg                                                                        >  1:405440/1‑151 (MQ=255)
                                                                         ttttGTGCGGTGGGGGCTCCGGCG‑ACATTCGCGGCTTGCCAGTTTATGGCCGCTGACC‑GGGTGTACCAAATA‑AATCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGt                                                                       <  2:968178/150‑1 (MQ=255)
                                                                          tttGTGCAGTGGTGGCTCCGGCG‑ACAATCGCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATT‑AATACCGCCAGCAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCCTAGTCACCTCGTGCGTTATAGCGTa                                                                      >  2:867109/1‑151 (MQ=255)
                                                                           tgttgCGGTGGTGGCCCCGGCG‑AAAACCGCTCCTTGCCAATTTTTGGCCGCTGAAC‑TGGCGTAGCAAATT‑AAACCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTAc                                                                     <  1:804452/148‑1 (MQ=255)
                                                                                cggggggggcccgggCG‑AAAACCGCTGTTTGCAAAATTTTGGCCGCTAAAC‑GGGGGTAGCAAATA‑AATCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCa                                                                <  1:690832/138‑1 (MQ=255)
                                                                                  gtggtgGCTCCGGCG‑ACAATCGCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCAc                                                               >  2:692854/1‑150 (MQ=255)
                                                                                  gtggtgGCTCCGGCG‑ACAATCGCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCAc                                                               >  2:586015/1‑150 (MQ=255)
                                                                                  gtggtgGCTCCGGCG‑ACAATCGCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCAc                                                               >  1:433120/1‑150 (MQ=255)
                                                                                  gtggtgGCTCCGGCG‑ACAATCGCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCAc                                                               >  1:394959/1‑150 (MQ=255)
                                                                                  gaggtggtGCCGGCC‑ACAGTCGCTGCTTGCCATTTTTGTGCGGGTGACG‑TGGGGGACAAAATT‑AATCCCGCCAAAAACACCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCAc                                                               <  2:436768/148‑1 (MQ=255)
                                                                                         ccccGGCG‑ACAATCGCTGCTGGCAAGATTATTGCCGCTGACC‑TGGCGTACCAAATT‑AATCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCAt                                                       <  1:526705/149‑1 (MQ=255)
                                                                                          ccgggCA‑ACACTCGTTGTTGGCAAGTTTATTGCCGCTGAAC‑GGGCGTAGCCAATT‑AATCCCGGCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCAtt                                                      <  1:979897/148‑1 (MQ=255)
                                                                                            cGGCG‑ACAACCCCTGCTTGCCAGTTTTTTCCCGCTGAAC‑TGGCGTAGCAAATT‑AATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTg                                                    <  1:732979/151‑1 (MQ=255)
                                                                                               cG‑ACAGTCGCTGCTCCCAATATTGTCGCCGGTGACG‑GGGCGTAGCAAATA‑AACCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTc                                                  <  2:595639/150‑1 (MQ=255)
                                                                                                g‑ACAATCGCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCg                                                >  1:186759/1‑151 (MQ=255)
                                                                                                     accgctctttGCCAAATTTTGGCCGCGAACC‑GGGGGTAGCCAATT‑AATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTtat                                            <  1:195091/143‑1 (MQ=255)
                                                                                                           gctTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGca                                      >  1:363888/1‑151 (MQ=255)
                                                                                                            tttGCCAGATTATGCCCGATAGAG‑TGGCGGAGCCAATA‑AACGCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGcac                                     <  2:1059730/150‑1 (MQ=255)
                                                                                                               cccAGTTTTTTGCCGCTAAAC‑GGGCGTAGCAAATT‑AATCCCCCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAg                                  <  1:918435/150‑1 (MQ=255)
                                                                                                               cccAGATTATGGCCGTGGAAC‑GGGGTTAGCAAATA‑AATCCCGCCGGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAg                                  <  2:81434/150‑1 (MQ=255)
                                                                                                               cccAGATTATGCCCGCTAAAC‑TGGCGAGGCCAATT‑AATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAg                                  <  1:501809/150‑1 (MQ=255)
                                                                                                                   gATTATTGCCGCTGAAC‑TGGCGTAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGCGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGAt                              >  2:327328/1‑151 (MQ=255)
                                                                                                                   gATTATTGCCGCTGAAC‑TGGCGTAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGAt                              >  2:77317/1‑151 (MQ=255)
                                                                                                                   gATTATTGCCGCTGAAC‑TGGCGTAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGAt                              >  1:568786/1‑151 (MQ=255)
                                                                                                                       aagagagaggAAC‑TGGCATAGCAAATA‑AATACCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAAc                          <  2:38874/142‑1 (MQ=255)
                                                                                                                              ctaaaC‑GGGGGTAGCAAATA‑AACCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGt                    <  1:1191127/147‑1 (MQ=255)
                                                                                                                              cggAAC‑GGGGGTACCAAATA‑AATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGt                    <  1:51206/148‑1 (MQ=255)
                                                                                                                              cTGACC‑TGGCGTAGCCATTT‑ATTCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGtt                   <  2:236640/151‑1 (MQ=255)
                                                                                                                                gAAC‑TGGCGTAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAAc                          >  1:38881/1‑142 (MQ=255)
                                                                                                                                gAAC‑TGGCGTAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGtttt                 >  2:7017/1‑151 (MQ=255)
                                                                                                                                      ggCGTAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCgg            >  1:1027222/1‑151 (MQ=255)
                                                                                                                                       gCGTAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCATAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCgga           >  1:737579/1‑151 (MQ=255)
                                                                                                                                            gCCAATT‑AATACTGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGt      >  2:334393/1‑151 (MQ=255)
                                                                                                                                                aTT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCa  >  1:945537/1‑151 (MQ=255)
                                                                                                                                                aTT‑AATACCGCCAGAAAAACCGCCAGCAGCGTGGGTGGAGGAATAGAAAACAGACGCATACTCCCCTCGTGCGTTACAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACTCACAAGTGATAAACTTCCGTTTTGCCGGAGGCGACGc   >  2:260306/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                
TACGCGATACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCG‑ACAATCGCTGCTTGCCAGATTATTGCCGCTGAAC‑TGGCGTAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCA  >  NC_000913/1505880‑1506172

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