Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,110,560	C→T	5.4%	G317D (GGC→GAC)	PP_5361 ←	protein P47K

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,110,560	0	C	T	5.4%	227.0 / 2.4	74	G317D (GGC→GAC)	PP_5361	protein P47K
Reads supporting (aligned to +/- strand):  ref base C (33/37);  new base T (2/2);  total (35/39)
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.76e-01

ATAATTTTTCGACGTCGATATTCTGCCCAATAAACACCAACTCCTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGAT  >  NC_002947/6110414‑6110707                                                                                                                                                    |                                                                                                                                                    aTAATTTTTCGACGTCGATATTCTGCCCAATAAACACCAACTCCTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGgcca                                                                                                                                                   <  4:173306/149‑1 (MQ=255) aTAATTTTTCGACGTCGATATTCTGCCCAATAAACACCAACTCCTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGgcca                                                                                                                                                   <  4:179348/149‑1 (MQ=255) aTAATTTTTCGACGTCGATATTCTGCCCAATAAACACCAACTCCTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGgc                                                                                                                                                     >  3:24416/1‑147 (MQ=255)  tAATTTTTCGCCGTCGATCTTCTGCCCAATAAACACCAACTCCTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGgccag                                                                                                                                                  <  8:263118/149‑1 (MQ=255)  tAATTTTTCGACGTCGATATTCTGCCCAATAAACACCAACTCCTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGgccag                                                                                                                                                  >  1:40062/1‑149 (MQ=255)    aTTTTTCGACGTCGATATTCTGCCCAATAAACACCAACTCCTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGgccagc                                                                                                                                                 <  6:99992/148‑1 (MQ=255)           gACGTCGATATTCTGCCCAATAAACACCAACTCCTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGTCAGCCAGACCAt                                                                                                                                         >  7:128558/1‑149 (MQ=255)           gACGTCGATATTCTGCCCAATAAACACCAACTCCTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGTCAGCCAGACCAt                                                                                                                                         >  5:20889/1‑149 (MQ=255)              gTCGATATTCTGCCCAATAAACACCAACTCCTGTCGGCAGTGACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGgc                                                                                                                                      <  7:233589/149‑1 (MQ=255)              gTCGATATTCTGCCCAATAAACACCAACTCCTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATgg                                                                                                                                       >  1:361069/1‑148 (MQ=255)                 gATATTCTGCCCAATAAACCCCAACTCCTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CCGCGGCCAGCCAGACCATGGCGca                                                                                                                                   >  7:319923/1‑149 (MQ=255)                 gATATTCTGCCCAATAAACACCAACTCCTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGgcgc                                                                                                                                    >  2:151091/1‑148 (MQ=255)                  atatTCTGCCCAATAAACACCAACTCCTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGcat                                                                                                                                  <  7:8010/149‑1 (MQ=255)                   tatTCTGCCCAATAAACACCAACTCCTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGcat                                                                                                                                  >  5:101879/1‑148 (MQ=255)                    atTCTGCCCAATAAACACCAACTCCTGTCGGCAGTCACCGGTGTCGCCTGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGcatca                                                                                                                                <  2:71233/149‑1 (MQ=255)                     ttCTGCCCAATAAACACCAACTCCTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGcatcat                                                                                                                               <  8:245362/149‑1 (MQ=255)                       cTGCCCAATAAACACCAACTCCTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGTACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATg                                                                                                                              <  5:39850/148‑1 (MQ=255)                          cccAATAAACACCAACTCCTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATgccgc                                                                                                                          <  4:240557/149‑1 (MQ=255)                          cccAATAAACACCAACTCCTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATgccgc                                                                                                                          >  2:296796/1‑149 (MQ=255)                            cAATAAACACCAACTCCTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGccgccc                                                                                                                        >  6:30175/1‑149 (MQ=255)                             aaTAAACACCAACTCCTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGccgcccg                                                                                                                       >  3:204691/1‑149 (MQ=255)                              aTAAACACCAACTCCTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGccgcccgc                                                                                                                      <  4:204691/149‑1 (MQ=255)                                  acacCACCTCCTCTCGGCCGTCACCGGTGTCGGCGGGCCAGTGTGTCACGCTCGCAATCGTGGTCTCCACATCCTGTGGCCACTGCGCGCTTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCGGCCCGCTTGc                                                                                                                  <  8:290661/149‑1 (MQ=255)                                  acacCAACTCCTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGc                                                                                                                  >  1:215046/1‑149 (MQ=255)                                   cacCAACTCCTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCg                                                                                                                 <  7:245669/149‑1 (MQ=255)                                    acCAACTCCTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATACACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTgag                                                                                                                 >  1:138669/1‑146 (MQ=255)                                        aCTCCTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAg                                                                                                            <  7:180734/149‑1 (MQ=255)                                          tCCTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCt                                                                                                          >  6:142071/1‑149 (MQ=255)                                            cTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGTCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGc                                                                                                        <  8:128558/149‑1 (MQ=255)                                            cTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGTCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGc                                                                                                        <  6:20889/149‑1 (MQ=255)                                            cTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGc                                                                                                        >  3:364364/1‑149 (MQ=255)                                            cTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGc                                                                                                        <  6:144442/149‑1 (MQ=255)                                               tCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAg                                                                                                     >  1:357176/1‑149 (MQ=255)                                               tCGGCAGGCACCGGTGTCGGCGGTCCAGGGGTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGCAGCGCCACGCAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGc                                                                                                        <  1:296796/147‑1 (MQ=255)                                                  gCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCtt                                                                                                  <  4:364364/149‑1 (MQ=255)                                                    aGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTcc                                                                                                <  2:357176/149‑1 (MQ=255)                                                     gTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCt                                                                                               >  4:350550/1‑149 (MQ=255)                                                        aCCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGat                                                                                            <  4:24416/149‑1 (MQ=255)                                                        aCCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGa                                                                                             <  5:30175/148‑1 (MQ=255)                                                        aCCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGa                                                                                             >  5:241458/1‑148 (MQ=255)                                                                 ggcggTCCAGTGTTTCTAGATCGCATTCGTGTTATGCACATCCTGTGGCCACTGCTCTCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTg                                                                                    <  4:337369/148‑1 (MQ=255)                                                                 ggcggTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTg                                                                                    <  5:142071/148‑1 (MQ=255)                                                                  gcggTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGc                                                                                  <  6:101879/149‑1 (MQ=255)                                                                      tCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAAc                                                                              <  3:350550/149‑1 (MQ=255)                                                                       ccAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGTGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAAc                                                                              >  2:38352/1‑148 (MQ=255)                                                                       ccAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAAcc                                                                             <  6:78137/149‑1 (MQ=255)                                                                          gtgtTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAg                                                                           >  2:170979/1‑148 (MQ=255)                                                                            gtTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGaa                                                                         >  1:289566/1‑148 (MQ=255)                                                                              ttCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGaaaaa                                                                      >  6:171032/1‑149 (MQ=255)                                                                                    ttCGAATTCGTGTTATCCAGATCCTGTGGCCTCTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCGGCCCGCTTTCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCAttt                                                                <  2:372180/148‑1 (MQ=255)                                                                                    aTCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCTCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCttt                                                                <  2:361069/149‑1 (MQ=255)                                                                                    aTCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCttt                                                                <  6:241458/149‑1 (MQ=255)                                                                                       gCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTgg                                                              <  2:40062/148‑1 (MQ=255)                                                                                        cATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGgccagcca                                                                                                                                               <  5:132018/66‑1 (MQ=255)                                                                                        cATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGgccagcca                                                                                                                                               >  6:132018/1‑66 (MQ=255)                                                                                          ttCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGGGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACg                                                          >  1:273504/1‑149 (MQ=255)                                                                                            cttgtTCTCCACCTCCTGTGGACACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACg                                                          <  8:319923/145‑1 (MQ=255)                                                                                             gtgtTATGCACATCCTGTGGCCAGTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCa                                                        <  7:315508/148‑1 (MQ=255)                                                                                             gtgtTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAg                                                       >  5:170883/1‑149 (MQ=255)                                                                                              tgtTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAg                                                       >  8:28763/1‑148 (MQ=255)                                                                                                      acaTCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTcc                                              >  4:33870/1‑149 (MQ=255)                                                                                                      acaTCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTcc                                              >  8:300540/1‑149 (MQ=255)                                                                                                        aTCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGt                                            <  3:59207/149‑1 (MQ=255)                                                                                                        aTCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGt                                            >  8:262529/1‑149 (MQ=255)                                                                                                           cTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGtttt                                         >  3:124043/1‑149 (MQ=255)                                                                                                            tgtgGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGCCCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTc                                        >  5:380791/1‑149 (MQ=255)                                                                                                              tgGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTcc                                       <  6:170883/148‑1 (MQ=255)                                                                                                              tgGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTc                                        <  6:364380/147‑1 (MQ=255)                                                                                                                gCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAg                                    >  6:376441/1‑149 (MQ=255)                                                                                                                  cACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTc                                  >  8:67177/1‑149 (MQ=255)                                                                                                                   aCTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCt                                                                  <  3:318694/116‑1 (MQ=255)                                                                                                                   aCTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCt                                                                  >  4:318694/1‑116 (MQ=255)                                                                                                                       cTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTtg                              <  1:170979/148‑1 (MQ=255)                                                                                                                                 aCGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAAc                                                                              <  8:93568/90‑1 (MQ=255)                                                                                                                                 aCGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAAc                                                                              >  7:93568/1‑90 (MQ=255)                                                                                                                                   gAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGt                 <  7:28763/149‑1 (MQ=255)                                                                                                                                     aGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTaa               <  1:151091/149‑1 (MQ=255)                                                                                                                                       cgcCAC‑CATCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTaaaa             <  4:124043/149‑1 (MQ=255)                                                                                                                                        cgcgagcaGTGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTaaaa             <  6:12974/143‑1 (MQ=255)                                                                                                                                         ccac‑caGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAACgc           >  8:336617/1‑149 (MQ=255)                                                                                                                                                  ggCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGa   >  7:236646/1‑149 (MQ=255)                                                                                                                                                  ggCCAGCCAGACCATGGCGCATCACGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGa   >  2:46922/1‑149 (MQ=255)                                                                                                                                                   gccagccaGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGAt  >  4:154363/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    ATAATTTTTCGACGTCGATATTCTGCCCAATAAACACCAACTCCTGTCGGCAGTCACCGGTGTCGGCGGTCCAGTGTTTCAAGATCGCATTCGTGTTATCCACATCCTGTGGCCACTGCTCGCGTGGCACGAAGCGCCAC‑CAGCGGCCAGCCAGACCATGGCGCATCATGCCGCCCGCTTGCGACCAGCTGCCAGCTTCCTGATATTTGCTGGCCAACCAGAAAAACCCTTTGGAACGCAGTAGCCTTCCGTTTTCCCAGTCCTTGTGGATAAAGTCGTAAAAGCGCTGTGGAT  >  NC_002947/6110414‑6110707

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