Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,337,473	C→G	6.4%	A38G (GCC→GGC)	PP_3812 →	lipase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,337,473	0	C	G	6.4%	189.6 / 3.7	62	A38G (GCC→GGC)	PP_3812	lipase
Reads supporting (aligned to +/- strand):  ref base C (29/29);  new base G (2/2);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

CCACTACCAACACTTGCACAGAAAAGGATTCGCCCATGAACATTGTCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGA  >  NC_002947/4337326‑4337619                                                                                                                                                      |                                                                                                                                                       ccACTACCAACACTTGCACAGAAAAGGATTCGCCCATGAACATTGTCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGcc                                                                                                                                                       >  8:327793/1‑149 (MQ=255)       ccAACACTTGCACAGAAAAGGATTCGCCCATGAACATTGTCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑ga                                                                                                                                                <  7:324300/148‑1 (MQ=255)        cAACACTTGCACAGAAAAGGATTCGCCCATTAACATTGTCCACAAAGCCTACTTACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑ga                                                                                                                                                <  1:122815/149‑1 (MQ=255)                 cacaGAAAAGGATTCGCCCATGAACATTGTCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGt                                                                                                                                      <  8:133119/148‑1 (MQ=255)                 cacaGAAAAGGATTCGCCCATGAACAATGTCCACAACGCC‑‑CTCACCACCTCCCTCTTGGCCATGACCGTTTCCAGTGCCTTCGCCGCAGCCAGCCCTGGTGTCGAACAACACACCCACGCCTTACTCGAAGCCCTT‑‑GAGAAAGGGCGTg                                                                                                                                     >  6:76865/1‑149 (MQ=255)                        aaGGATTCGCCCATGAACATTGTCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGcc                                                                                                                              >  4:290451/1‑149 (MQ=255)                         aGGATTCGCCCATGAACATTGTCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCCAGCCg                                                                                                                             >  8:151049/1‑149 (MQ=255)                                gCCCATGAACATTGTCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAAc                                                                                                                      <  3:230823/149‑1 (MQ=255)                                  ccATGAACATTGTCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGGCTTCGAGGAAGGCCTg‑‑gg                                                                                                                                                >  4:327556/1‑120 (MQ=255)                                  ccATGAACATTGTCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGGCTTCGAGGAAGGCCTg‑‑gg                                                                                                                                                <  3:327556/121‑2 (MQ=255)                                   cATGAACATTGTCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGc                                                                                                                   >  4:323125/1‑149 (MQ=255)                                       aaCATTGTCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTca                                                                                                                <  3:290451/148‑1 (MQ=255)                                         cATTGTCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCacac                                                                                                             <  2:87753/149‑1 (MQ=255)                                         cATTGTCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCacac                                                                                                             <  7:221962/149‑1 (MQ=255)                                          aTTGTCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCacac                                                                                                             >  5:26757/1‑148 (MQ=255)                                              tCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGGCTTCGAGGAAGGCCTg‑‑ggtgtg                                                                                                                                            <  2:314893/113‑6 (MQ=38)                                              tCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGGCTTCGAGGAAGGCCTg‑‑ggtgtg                                                                                                                                            >  1:314893/1‑108 (MQ=38)                                               ccACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAac                                                                                                        >  6:128018/1‑148 (MQ=255)                                                cacaAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAacgt                                                                                                      >  2:211184/1‑148 (MQ=255)                                                cacaAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAacga                                                                                                      >  8:171611/1‑149 (MQ=255)                                                 acaAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAacgac                                                                                                     >  2:51232/1‑149 (MQ=255)                                                 acaAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAacgac                                                                                                     <  1:36278/149‑1 (MQ=255)                                                   aaaGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGc                                                                                                   <  7:327793/149‑1 (MQ=255)                                                   aaaGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGc                                                                                                   >  4:117681/1‑149 (MQ=255)                                                        c‑‑cTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTg                                                                                                                          >  8:61034/1‑121 (MQ=255)                                                        c‑‑cTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTg                                                                                                                          <  7:61034/121‑1 (MQ=255)                                                              accaccTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCgg                                                                                           >  3:286037/1‑148 (MQ=255)                                                                              cccTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGTCGCCCGCGCGGTACTGA‑CCG‑GTGCCCAg                                                                        <  1:211184/149‑1 (MQ=255)                                                                                cTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGc                                                                      >  6:136591/1‑149 (MQ=255)                                                                                     cGTTTCCAGTGCCTTCGCCTCATGCAGCCCTGGTGTCGTACAACACACCCAGGCCTTCGTCGAATCCCTG‑‑GAGAAATGGGTTGGCAAGCCGCTGTAACAGCTCACACCAAACTACGCCCGCGCGGTACTGA‑CCT‑GTGCCCAGGCTTCg                                                                  <  2:62061/148‑1 (MQ=255)                                                                                     cGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCg                                                                  <  7:151049/148‑1 (MQ=255)                                                                                      gTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGt                                                                <  8:108345/149‑1 (MQ=255)                                                                                       tttCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGt                                                                >  6:2605/1‑148 (MQ=255)                                                                                          ccAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAg                                                            >  4:344392/1‑149 (MQ=255)                                                                                           cAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCACGCCGCGGGAACAGCTCCCACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAgg                                                           >  7:159182/1‑149 (MQ=255)                                                                                            aGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGt                                                          >  4:320319/1‑149 (MQ=255)                                                                                                ccTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGTCGAc                                                      <  5:128018/149‑1 (MQ=255)                                                                                                ccTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGTCGAc                                                      <  3:117681/149‑1 (MQ=255)                                                                                                  ttCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTGGAGAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGTCGAc                                                      >  1:108822/1‑149 (MQ=255)                                                                                                      ccgcAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCtt                                                                    <  4:71777/129‑1 (MQ=255)                                                                                                      ccgcAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCtt                                                                    >  3:71777/1‑129 (MQ=255)                                                                                                         caggcagCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGTCGACCTATCgggg                                             <  7:171611/149‑1 (MQ=255)                                                                                                               gCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAg                                       <  1:51232/149‑1 (MQ=255)                                                                                                                        tCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAgg                                                           <  2:80683/120‑1 (MQ=255)                                                                                                                        tCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAgg                                                           >  1:80683/1‑120 (MQ=255)                                                                                                                          gAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCGTGTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAAc                             <  1:94576/149‑1 (MQ=255)                                                                                                                          gAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACg                            >  4:290944/1‑149 (MQ=255)                                                                                                                           aacaacACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACg                            >  1:369159/1‑148 (MQ=255)                                                                                                                              aacaCACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGACCCG‑GTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAgg                                      <  1:263617/136‑1 (MQ=255)                                                                                                                              aacaCACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAgg                                      >  2:263617/1‑135 (MQ=255)                                                                                                                              aacaCACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAAc                             >  2:371234/1‑144 (MQ=255)                                                                                                                              aacaCACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAAc                             <  1:371234/144‑1 (MQ=255)                                                                                                                                 acacCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGAt                      >  6:93322/1‑148 (MQ=255)                                                                                                                                   acCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGAtt                     >  1:143823/1‑147 (MQ=255)                                                                                                                                   acCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCa                   <  5:76865/149‑1 (MQ=255)                                                                                                                                    cccAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGGCCGCGCGGTACTGA‑CCG‑GTGCCGAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAg                  <  3:144349/149‑1 (MQ=255)                                                                                                                                      cAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGaa                                  >  7:296477/1‑131 (MQ=255)                                                                                                                                      cAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGaa                                  <  8:296477/131‑1 (MQ=255)                                                                                                                                           cTTCCTCGAAGCCCTT‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCACGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATg           <  8:167497/149‑1 (MQ=255)                                                                                                                                               cTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGc         >  5:328328/1‑147 (MQ=255)                                                                                                                                               cTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCa       <  6:26757/149‑1 (MQ=255)                                                                                                                                                tCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAg      >  7:250195/1‑149 (MQ=255)                                                                                                                                                    aGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGa  <  6:328328/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                       CCACTACCAACACTTGCACAGAAAAGGATTCGCCCATGAACATTGTCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAGTCGA  >  NC_002947/4337326‑4337619

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