Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,989,017	A→T	7.4%	D512V (GAC→GTC)	oplA →	5‑oxoprolinase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,989,017	0	A	T	7.4%	189.8 / 4.4	67	D512V (GAC→GTC)	oplA	5‑oxoprolinase A
Reads supporting (aligned to +/- strand):  ref base A (29/32);  new base T (2/3);  total (32/35)
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.40e-01

TGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGATTCATGAGTGCAGTGTGGA  >  NC_002947/3988870‑3989164                                                                                                                                                     |                                                                                                                                                    tggtgCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGAc                                                                                                                                                    >  2:23827/1‑149 (MQ=255) tggtgCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGAc                                                                                                                                                    <  4:171420/149‑1 (MQ=255)      cgcCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACCTACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTgg                                                                                                                                                <  1:84771/149‑1 (MQ=255)            gCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCAGGTCGACGTGCCCGGCGTCCAGGCGCat                                                                                                                                           <  5:222669/147‑2 (MQ=38)            gCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCAGGTCGACGTGCCCGGCGTCCAGGCGCat                                                                                                                                           <  5:121480/147‑2 (MQ=38)             ccGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACt                                                                                                                                         >  7:143665/1‑148 (MQ=255)               gaagaaATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCTCCGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCATGGACGCCGGGCACGTCGCCCTGGCGCACCTCaa                                                                                                                                      >  2:233932/1‑149 (MQ=255)               gaagaaATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCaa                                                                                                                                      >  4:81414/1‑149 (MQ=255)               gaagaaATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCACAACTCaa                                                                                                                                      >  5:139440/1‑149 (MQ=255)                    aaTTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAgg                                                                                                                                 <  1:319/149‑1 (MQ=255)                     aTTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGc                                                                                                                                <  4:279567/149‑1 (MQ=255)                                cacCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGaa                                                                                                                      >  6:102450/1‑148 (MQ=255)                                cacCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAAc                                                                                                                     <  8:143665/149‑1 (MQ=255)                                   cGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCATTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTAc                                                                                                                  <  5:7142/149‑1 (MQ=255)                                       cTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTc                                                                                                               <  1:23827/148‑1 (MQ=255)                                          aTCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGc                                                                                                            <  6:143230/148‑1 (MQ=255)                                           tCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCaa                                                                                                          <  3:81414/149‑1 (MQ=255)                                            ccccAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCAGGTCGACGTGCCCGGCGTCCAGgc                                                                                                                                               <  8:286200/111‑1 (MQ=37)                                            ccccAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCAGGTCGACGTGCCCGGCGTCCAGgc                                                                                                                                               >  7:286200/1‑111 (MQ=37)                                             cccAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAgg                                                                                                        >  8:188450/1‑149 (MQ=255)                                             cccAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAg                                                                                                         <  3:64888/148‑1 (MQ=255)                                              ccAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGc                                                                                                       >  6:94143/1‑149 (MQ=255)                                                      gCCTCGGGCTTGTGCGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAAGCCTTGCCaa                                                                                               >  1:85783/1‑149 (MQ=255)                                                         tCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCtt                                                                                                    <  3:221771/141‑1 (MQ=255)                                                         tCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCtt                                                                                                    >  4:221771/1‑141 (MQ=255)                                                             gCTTGTGTGCGTTAGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCTCCGAACTACGTCAGCAAGGCCTTGCCAACCTgcgc                                                                                        <  7:189516/149‑1 (MQ=255)                                                               ttGTGTGCGTTCGGCCAGATCCTGTCCGACGAGCGCTA‑CGACCAGCTCACCAGCCTGTCCATCCGCCTGGACGCCGGGCACTACGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACCTCAGCAAGGCCTTGCCAACCTccgcc                                                                                       >  2:79556/1‑147 (MQ=255)                                                                    tgCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTg                                                                                                   <  8:79565/131‑1 (MQ=255)                                                                    tgCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTg                                                                                                   >  7:79565/1‑131 (MQ=255)                                                                      cGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGc                                                                               >  1:270341/1‑149 (MQ=255)                                                                       gTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCt                                                                              >  8:37334/1‑149 (MQ=255)                                                                                  tCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAg                                                                   <  7:37334/149‑1 (MQ=255)                                                                                  tCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAg                                                                   <  5:260114/149‑1 (MQ=255)                                                                                   ccTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAgg                                                                  >  6:109090/1‑149 (MQ=255)                                                                                    cTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGc                                                                 >  1:104974/1‑149 (MQ=255)                                                                                         cGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCcag                                                              >  5:153520/1‑147 (MQ=255)                                                                                         cGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCcag                                                              >  5:162554/1‑147 (MQ=255)                                                                                                 gcTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCAc                                                    <  2:270341/149‑1 (MQ=255)                                                                                                   tA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTa                                                  >  7:219168/1‑149 (MQ=255)                                                                                                         ccAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTAAGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACAcc                                              >  8:128388/1‑148 (MQ=255)                                                                                                           aGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTg                                           <  5:94143/149‑1 (MQ=255)                                                                                                             cTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGga                                         <  5:109090/149‑1 (MQ=255)                                                                                                              tCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTg                                           >  5:187300/1‑146 (MQ=255)                                                                                                                aCCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGAt                                      >  2:55412/1‑149 (MQ=255)                                                                                                                  cAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGaa                                                                                     <  5:253101/100‑1 (MQ=255)                                                                                                                  cAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGaa                                                                                     >  6:253101/1‑100 (MQ=255)                                                                                                                    gCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATgcgc                                  >  2:315013/1‑149 (MQ=255)                                                                                                                     ccTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAAGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCt                                 <  2:85783/149‑1 (MQ=255)                                                                                                                         tCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCt                             <  1:49085/149‑1 (MQ=255)                                                                                                                           gATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGATGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTa                                                  <  1:51191/126‑1 (MQ=255)                                                                                                                           gATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTa                                                  >  2:51191/1‑126 (MQ=255)                                                                                                                            aTGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTggg                          >  4:292329/1‑149 (MQ=255)                                                                                                                            aTGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTggg                          <  6:162554/149‑1 (MQ=255)                                                                                                                              gcgcCAGGTCGACGTGCCCGGCGTCCAGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTgggg                         >  5:162412/1‑148 (MQ=255)                                                                                                                               cgcCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGc                        <  1:233932/148‑1 (MQ=255)                                                                                                                                gcCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTgggg                         >  2:227057/1‑146 (MQ=255)                                                                                                                                 ccTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGa                     <  2:42621/149‑1 (MQ=255)                                                                                                                                      aCGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGATTCAt                <  6:153520/149‑1 (MQ=255)                                                                                                                                       cGCCGGGCACGTCTACCTGGCGCAACTCAACCAGGCACTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGCCCAAAAGTTGCCACTACACACTGGAGATGCGCTACCTGGGGCAGATTCAt                >  8:249860/1‑148 (MQ=255)                                                                                                                                        gCCGGGCACGTCGACCTGGCACAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGATTCATGa              <  6:139440/149‑1 (MQ=255)                                                                                                                                             gCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGATTCATGAGTg           >  3:19433/1‑147 (MQ=255)                                                                                                                                             gCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGATTCATGAGTGCa         <  1:55412/149‑1 (MQ=255)                                                                                                                                               aCGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGATTCATGAGTGCAgt       <  7:188450/149‑1 (MQ=255)                                                                                                                                                cGTCGACCTGGTGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGATTCATGAGTGCAgtg      >  1:76471/1‑149 (MQ=255)                                                                                                                                                  tCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGATTCATGAGTGCAGtgtg    <  3:292329/149‑1 (MQ=255)                                                                                                                                                    gACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTTCACCCTGGAGATGCGCTACCTGGGGCAGATTCATGAGTGCAGTGTGGa  <  1:79556/149‑1 (MQ=255)                                                                                                                                                    gACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGATTCATGAGTGCAGTGTGGa  <  5:102450/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGATTCATGAGTGCAGTGTGGA  >  NC_002947/3988870‑3989164

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