Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	807,900	A→T	6.5%	L45Q (CTG→CAG)	PP_0694 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	807,900	0	A	T	6.5%	187.6 / 2.2	61	L45Q (CTG→CAG)	PP_0694	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (32/25);  new base T (2/2);  total (34/27)
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.80e-01

TGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCATGGGCTGGGCTCGCGT  >  NC_002947/807756‑808048                                                                                                                                                 |                                                                                                                                                      tGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑Gt                                                                                                                                                  >  1:228482/1‑149 (MQ=255) tGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑Gt                                                                                                                                                  >  2:324463/1‑149 (MQ=255)      aGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCgg                                                                                                                                              <  1:173696/148‑1 (MQ=255)      aGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCgg                                                                                                                                              <  2:254562/148‑1 (MQ=255)        cAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTAc                                                                                                                                           >  3:65299/1‑149 (MQ=255)         aaGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACt                                                                                                                                          <  6:264093/149‑1 (MQ=255)          aGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTg                                                                                                                                         <  2:228482/149‑1 (MQ=255)            gTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGt                                                                                                                                       >  8:207/1‑149 (MQ=255)            gTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGt                                                                                                                                       >  1:40205/1‑149 (MQ=255)               cTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCt                                                                                                                                    >  4:245153/1‑149 (MQ=255)                tCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTg                                                                                                                                   >  3:324278/1‑149 (MQ=255)                   aCCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGaa                                                                                                                                <  1:188009/149‑1 (MQ=255)                      aGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGtt                                                                                                                             <  4:324278/149‑1 (MQ=255)                            aaGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAg                                                                                                                       >  6:242516/1‑149 (MQ=255)                                tGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTgg                                                                                                                   <  3:245153/149‑1 (MQ=255)                                 gAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGa                                                                                                                                 >  8:168479/1‑134 (MQ=255)                                 gAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGa                                                                                                                                 <  7:168479/134‑1 (MQ=255)                                 gAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTggg                                                                                                                  >  4:173303/1‑149 (MQ=255)                                             gcgcTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAg                                                                                                      >  5:123857/1‑149 (MQ=255)                                              cgcTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAgg                                                                                                     <  2:40205/149‑1 (MQ=255)                                                        tCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCgg                                                                                           <  2:234612/149‑1 (MQ=255)                                                          ttgttgACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCt                                                                                         >  6:181838/1‑149 (MQ=255)                                                                aCCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGt                                                                                    <  4:147269/148‑1 (MQ=255)                                                                aCCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGt                                                                                    <  1:242729/148‑1 (MQ=255)                                                                  cTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTgg                                                                                  >  4:125892/1‑148 (MQ=255)                                                                  cTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTgg                                                                                  >  6:228628/1‑148 (MQ=255)                                                                   tGCGCGGCGGTCATGTCCGGCGCGATGTAGTAGCAACGACGGGTATCGCCAACCTGTTCGGTGGGCAGGTGCGTGATGAGG‑GTGGGGTACTGGTCCTGGGTGATAGACAGCTGGGGGTGGAGGTGGGGGTCGGGGCGGCTGACGTgg                                                                                  >  2:223027/1‑147 (MQ=255)                                                                        ggcggTCATGTCCGGCACGATGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCtt                                                                           >  2:122645/1‑149 (MQ=255)                                                                              cATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGaa                                                                       >  7:8241/1‑147 (MQ=255)                                                                                tGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGc                                                                   <  1:324463/149‑1 (MQ=255)                                                                                  tCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTc                                                                 >  5:92366/1‑149 (MQ=255)                                                                                            agtagtAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCgg                                                             >  4:61533/1‑143 (MQ=255)                                                                                               agtagAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGgc                                                    >  3:217843/1‑149 (MQ=255)                                                                                                gtagAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAAc                                                                      <  3:127179/130‑1 (MQ=255)                                                                                                gtagAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAAc                                                                      >  4:127179/1‑130 (MQ=255)                                                                                                gtagAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGgcg                                                   <  6:92366/149‑1 (MQ=255)                                                                                                  agAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGgcgc                                                  >  5:255597/1‑148 (MQ=255)                                                                                                   gAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGgcgc                                                  >  3:145406/1‑147 (MQ=255)                                                                                                      aCCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACTGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCt                                             <  4:65299/149‑1 (MQ=255)                                                                                                      aCCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCt                                             >  5:23810/1‑149 (MQ=255)                                                                                                        cACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTgg                                           <  3:173303/149‑1 (MQ=255)                                                                                                        cACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTACTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAATATCGGCGCGAGCTgg                                           >  6:261871/1‑149 (MQ=255)                                                                                                           gggTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGa                                         >  7:224855/1‑148 (MQ=255)                                                                                                           gggTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGAc                                        <  6:74687/149‑1 (MQ=255)                                                                                                           cagTACCGCAC‑CCTGTTCAAGGCCCTGGGCTTCGAACTGGTGTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGcgcg                                                 >  8:10509/3‑140 (MQ=25)                                                                                                           cagTACCGCAC‑CCTGTTCAAGGCCCTGGGCTTCGAACTGGTGTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGcgcg                                                 <  7:10509/138‑1 (MQ=25)                                                                                                           cagTACCGCAC‑CCTGTTCAAGGCCCTGGGCTTCGAACTGGTGTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGa                                         >  3:140471/3‑148 (MQ=25)                                                                                                                  gCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCgg                                 <  4:217843/149‑1 (MQ=255)                                                                                                                   cacaCCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCgg                                 >  2:110660/1‑148 (MQ=255)                                                                                                                   cacaCCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGa                                <  6:23810/149‑1 (MQ=255)                                                                                                                    acacCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGaa                               <  5:181838/149‑1 (MQ=255)                                                                                                                    acacCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCCAGCTGGGACACTACGGa                                >  6:236010/1‑148 (MQ=255)                                                                                                                     cacCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGaaa                              <  8:8241/149‑1 (MQ=255)                                                                                                                       ccAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTgg                                           <  7:25003/134‑1 (MQ=255)                                                                                                                       ccAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTgg                                           >  8:25003/1‑134 (MQ=255)                                                                                                                              gAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGAt                     >  7:173536/1‑149 (MQ=255)                                                                                                                               aggCCCTGGGCTTCGAACTGGTGTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGAt                     <  4:140471/147‑1 (MQ=255)                                                                                                                                  cccAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGTCAGCTGTGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCTAGCTGGGACAGTTCGGAAAGATCAAGATGCAt                 <  7:207/149‑1 (MQ=255)                                                                                                                                  cccAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCAt                 <  3:61533/149‑1 (MQ=255)                                                                                                                                    cAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCAtgg               >  4:14780/1‑149 (MQ=255)                                                                                                                                              aaCAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCATGGGCTGGGCTcg     >  6:290824/1‑149 (MQ=255)                                                                                                                                                 aGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCATGGGCTGGGCTCGCGt  >  6:92378/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                      TGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCATGGGCTGGGCTCGCGT  >  NC_002947/807756‑808048

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