Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,337,467	T→A	6.2%	L36H (CTC→CAC) ‡	PP_3812 →	lipase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,337,467	0	T	A	6.2%	193.2 / 2.6	64	L36H (CTC→CAC) ‡	PP_3812	lipase
Reads supporting (aligned to +/- strand):  ref base T (29/31);  new base A (2/2);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01

AAGACGCCCACTACCAACACTTGCACAGAAAAGGATTCGCCCATGAACATTGTCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAG  >  NC_002947/4337319‑4337615                                                                                                                                                       |                                                                                                                                                         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                                                                                                                                            >  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                         <  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)                                                                                                                                                  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)                                                                                                                                                       |                                                                                                                                                         AAGACGCCCACTACCAACACTTGCACAGAAAAGGATTCGCCCATGAACATTGTCCACAAAGCC‑‑CTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTG‑‑GAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGA‑CCG‑GTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGGGGTCGAGGTGAAGGAACGCACGATTCAGGCCAATGGCCAG  >  NC_002947/4337319‑4337615

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