Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,337,461	C→G	6.1%	A34G (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,461	0	C	G	6.1%	206.9 / 3.3	65	A34G (GCC→GGC)	PP_3812	lipase
Reads supporting (aligned to +/- strand):  ref base C (29/32);  new base G (2/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

CAGGCAAGACGCCCACTACCAACACTTGCACAGAAAAGGATTCGCCCATGAACATTGTCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAG  >  NC_002947/4337314‑4337603                                                                                                                                                      |                                                                                                                                                   cAGGCATGACGCCCACTACCAACACTTGCACAGAAAAGGATTCGCCCATGAACATTGTCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGcc                                                                                                                                                   >  3:73937/1‑149 (MQ=255)     cAAGACGCCCACTACCAACACTTGCACAGAAAAGGATTCGCCCATGAACATTGTCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCTGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTcc                                                                                                                                               <  1:284260/149‑1 (MQ=255)      aaGACGCCCACTACCAACACTTGCACAGAAAAGGATTCGCCCATGAACATTGTCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCt                                                                                                                                              <  7:92564/149‑1 (MQ=255)       aGACGCCCACTACCAACAATTGCACAGACACGGATTCGCCCATGAACATTGTCCACAAAGCC‑‑CTCACCAACTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCATCCCTTTTGTCGAAAAACACACCCAGTCCTTCCTc                                                                                                                                             <  1:233080/149‑1 (MQ=255)          cGCCCACTACCAACACTTGCACAGAAAAGGATTCGCCCATGAACATTGTCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCTGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGaa                                                                                                                                          <  1:89618/149‑1 (MQ=255)          cGCCCACTACCAACACTTGCACAGAAAAGGATTCGCCCATGAACATTGTCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGaa                                                                                                                                          <  4:73937/149‑1 (MQ=255)             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                                                                                                                                <  3:327556/121‑2 (MQ=255)                                              ccATGAACATTGTCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGGCTTCGAGGAAGGCCTg‑‑gg                                                                                                                                >  4:327556/1‑120 (MQ=255)                                               cATGAACATTGTCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGc                                                                                                   >  4:323125/1‑149 (MQ=255)                                                   aaCATTGTCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTca                                                                                                <  3:290451/148‑1 (MQ=255)                                                     cATTGTCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCacac                                                                                             <  7:221962/149‑1 (MQ=255)                                                     cATTGTCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCacac                                                                                             <  2:87753/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)                                                                           ccaccTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTc                                                                                                                                                <  5:342195/80‑1 (MQ=255)                                                                           ccaccTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTc                                                                                                                                                >  6:342195/1‑80 (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                                      <  3:117681/149‑1 (MQ=255)                                                                                                            ccTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGTCGAc                                      <  5:128018/149‑1 (MQ=255)                                                                                                              ttCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTGGAGAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGTCGAc                                      >  1:108822/1‑149 (MQ=255)                                                                                                                  ccgcAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCtt                                                    >  3:71777/1‑129 (MQ=255)                                                                                                                  ccgcAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCtt                                                    <  4:71777/129‑1 (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             <  1:371234/144‑1 (MQ=255)                                                                                                                                          aacaCACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAAc             >  2:371234/1‑144 (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)                                                                                                                                                      |                                                                                                                                                   CAGGCAAGACGCCCACTACCAACACTTGCACAGAAAAGGATTCGCCCATGAACATTGTCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAG  >  NC_002947/4337314‑4337603

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: ‑