Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,691,381	A→C	7.1%	E110A (GAG→GCG)	PP_1487 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,691,381	0	A	C	7.1%	157.3 / 2.3	54	E110A (GAG→GCG)	PP_1487	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/22);  new base C (2/2);  total (30/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

GCAAGTACGGTGGCAAGGCCCTGACTTACGGCGGCCTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAGGGTGAATTCACGCGCCTGGACAGCGACCGCGAATTGCAGCACTGGCTGCATGCCGAGCT  >  NC_002947/1691234‑1691528                                                                                                                                                     |                                                                                                                                                    gCAAGTACGGTGGCAAGGCCCTGACTTACGGCGGCCTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGag                                                                                                                                                    <  1:170783/149‑1 (MQ=255) gCAAGTACGGTGGCAAGGCCCTGACTTACGGCGGCCTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGag                                                                                                                                                    >  7:144281/1‑149 (MQ=255)     gTACGGTGGCAAGGCCCTGACTTACGGCGGCCTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGagcagc                                                                                                                                                <  6:326379/149‑1 (MQ=255)     gTACGGTGGCAAGGCCCTGACTTACGGCGGCCTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGagcag                                                                                                                                                 >  4:195240/1‑148 (MQ=255)     gTACGGTGGCAAGGCCCTGACTTACGGCGGCCTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGagcag                                                                                                                                                 >  3:151603/1‑148 (MQ=255)     gTACGGTGGCAAGGCCCTGACTTACGGCGGCCTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGCCACCAGGGCGcgcag                                                                                                                                                 >  2:250586/1‑149 (MQ=255)             gCAAGGCCCTGACTTACGGCGGCCTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAg                                                                                                                                        >  6:57112/1‑149 (MQ=255)               aaGGCCCTGACTTACGGCGGCCTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGc                                                                                                                                      >  6:232415/1‑149 (MQ=255)               aaGGCCCTGACTTACGGCGGCCTCGCCGCACTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGc                                                                                                                                      >  6:318184/1‑149 (MQ=255)                        aCTTACGGCGGCCTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTgc                                                                                                                              <  2:333314/148‑1 (MQ=255)                        aCTTACGGCGGCCTCGCCGCACTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTgc                                                                                                                              <  5:318184/148‑1 (MQ=255)                          ttACGGCGGCCTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTgcgcg                                                                                                                           >  4:245699/1‑149 (MQ=255)                          ttACGGCGGCCTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTgcgcg                                                                                                                           >  8:154171/1‑149 (MQ=255)                              ggcggcCTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCt                                                                                                                       <  7:281303/149‑1 (MQ=255)                               gcggcCTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCtg                                                                                                                      <  8:166402/149‑1 (MQ=255)                                   ccTCGCCGCGCTGGGGGACCTGGCCTATAAACCTTACGGCAACTGCCCGGCAAACCCGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCtggtgg                                                                                                                  <  7:114834/149‑1 (MQ=255)                                       gccgcGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTggcgg                                                                                                               <  4:151603/148‑1 (MQ=255)                                       gccgcGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGCCACCAGGGCGCGCAGTACCGCCTGGCTGTGCTGCTCGGCCtgggccgg                                                                                                                <  1:250586/148‑6 (MQ=21)                                        ccgcGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTggcggcg                                                                                                             >  2:67178/1‑149 (MQ=255)                                        ccgcGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTggcggcg                                                                                                             >  3:313146/1‑149 (MQ=255)                                          gcgcTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTggcggcggc                                                                                                           >  8:217672/1‑149 (MQ=255)                                                       tGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATggg                                                                                              >  3:259661/1‑149 (MQ=255)                                                       tGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATggg                                                                                              >  1:84386/1‑149 (MQ=255)                                                          ccTATAAAGGTTTCGGGAACTTCCAGGCAAACCAGGGGGGGGCCGCAGCCGAGCCTTAAACCATCGACCGCCTGCCTCGGGC‑CCAGGCCGAGCAGCACAGCCAGGCGTTACTGCGCGCCCTGGTGGCGGCGTCCAAGTCAGATGGGCat                                                                                           <  1:77341/149‑1 (MQ=255)                                                             ataAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATc                                                                                        >  7:94031/1‑149 (MQ=255)                                                              taAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATc                                                                                        <  8:94031/148‑1 (MQ=255)                                                                aaGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATcgac                                                                                     <  1:67178/149‑1 (MQ=255)                                                                              tGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGgc                                                                       >  7:222687/1‑149 (MQ=255)                                                                               ggcaggcaAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGgcg                                                                      >  3:45684/1‑149 (MQ=255)                                                                                gcaggcaAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGgcgc                                                                     >  5:254034/1‑149 (MQ=255)                                                                                gcaggcaAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTCGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGgcgc                                                                     >  4:32144/1‑149 (MQ=255)                                                                                    gcaAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGAt                                                                 >  5:346508/1‑149 (MQ=255)                                                                                         ccAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAgg                                                            <  7:217672/149‑1 (MQ=255)                                                                                          cAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATg                                                                                                >  1:160471/1‑112 (MQ=255)                                                                                          cAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATg                                                                                                <  2:160471/112‑1 (MQ=255)                                                                                           aGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAGGGt                                                          <  4:313146/149‑1 (MQ=255)                                                                                            gCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAGGGt                                                          >  6:202665/1‑148 (MQ=255)                                                                                                   gCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCCGC‑CCCGGCCGCGCAGCCCAGCCCGCCGGTACTGCGCCCCCTGGGGGCGGCGGCCAAGACCCATGGGGATAACGACGACCCCCAGCGGGCCCTGATAGCGGGTGAATTCAc                                                  >  2:338140/1‑149 (MQ=255)                                                                                                          ccGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAGGGTGAATTCACGCGCCTg                                           <  6:346508/149‑1 (MQ=255)                                                                                                             aGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAGGGTGAATTCACGCGCCTGGa                                         >  2:195782/1‑148 (MQ=255)                                                                                                              gCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAGGGTGAATTCACGCGCCTGGACa                                       <  3:32144/149‑1 (MQ=255)                                                                                                               ccTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAGGGTGAATTCACGCGCCTGGACAg                                      <  4:45684/149‑1 (MQ=255)                                                                                                                cTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAGGGTGAATTCACGCGCCTGGACAGc                                     <  2:84386/149‑1 (MQ=255)                                                                                                                 tCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAGGGTGAATTCACGCGCCTGGACAGc                                     >  4:233839/1‑148 (MQ=255)                                                                                                                 tCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAGGGTGAATTCACGCGCCTGGACAGCg                                    >  7:125401/1‑149 (MQ=255)                                                                                                                   aaaCCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAGGGTGAATTCACGCGCCTGGACAGCGAc                                  >  6:189077/1‑149 (MQ=255)                                                                                                                           gACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAGGGTGAATTCACGCGCCTGGACAGCGACCGCGAAtt                          <  5:202665/149‑1 (MQ=255)                                                                                                                           gACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAGGGTGAATTCACGCGCCTGGACAGCGACCGCGAAt                           <  4:259661/148‑1 (MQ=255)                                                                                                                            aCCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAGGGTGAATTCACGCGCCTGGACAGCGACCGCGAATTg                         >  3:324745/1‑149 (MQ=255)                                                                                                                              cGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAGGGTGAATTCACGCGCCTGGACAGCGACCGCGAATTgca                       <  5:232415/149‑1 (MQ=255)                                                                                                                                    ccgccgGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAGGGTGAATTCACGCGCCTGGACAGCGACCGCGAATTGCAGCACTg                 >  6:175849/1‑149 (MQ=255)                                                                                                                                        cgGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAGGGTGAATTCACGCGCCTGGACAGCGACCGCGAATTGCAGCACTGGCTg             <  8:125401/149‑1 (MQ=255)                                                                                                                                          gC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAGGGTGAATTCACGCGCCTGGACAGCGACCGCGAAGTGCAGCACTGGCTGCa           >  7:246841/1‑149 (MQ=255)                                                                                                                                          ccaccAGGGCGCGCAGTACCGCCTGGCTGTGCTGCTCGGCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAGGGTGAATTCAc                                                  <  8:281206/108‑1 (MQ=21)                                                                                                                                             ccAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAGGGTGAATTCACGCGCCTGGACAGCGACCGCGAATTGCAGCAc                   <  3:169252/139‑1 (MQ=255)                                                                                                                                             ccAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAGGGTGAATTCACGCGCCTGGACAGCGACCGCGAATTGCAGCAc                   >  4:169252/1‑139 (MQ=255)                                                                                                                                               aGGGCGCGCAGTACCGCCTGGCTGTGCTGCTCGGCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAGGGTGAATTCACGCGCCTGGACAGCGACCGCGAATTGCAGCACTGGCTGCATGcc       >  3:31034/1‑149 (MQ=38)                                                                                                                                                ggCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAGGGTGAATTCACGCGCCTGGACAGCGACCGCGAATTGCAGCACTGGCTGCATGCCg      <  3:245699/149‑1 (MQ=255)                                                                                                                                                   cGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAGGGTGAATTCACGCGCCTGGACAGCGACCGCGAATTGCAGCACTGGCTGCATGCCGAGc   >  2:244137/1‑149 (MQ=255)                                                                                                                                                    gAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAGGGTGAATTCACGCGCCTGGACAGCGACCGCGAATTGCAGCACTGGCTGCATGCCGAGCt  >  3:147354/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GCAAGTACGGTGGCAAGGCCCTGACTTACGGCGGCCTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGC‑CCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAGGGTGAATTCACGCGCCTGGACAGCGACCGCGAATTGCAGCACTGGCTGCATGCCGAGCT  >  NC_002947/1691234‑1691528

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