Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,989,002	A→T	8.1%	D507V (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,002	0	A	T	8.1%	176.3 / 4.3	61	D507V (GAC→GTC)	oplA	5‑oxoprolinase A
Reads supporting (aligned to +/- strand):  ref base A (29/27);  new base T (2/3);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

CGGCGGGCATGCACGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGATTCAT  >  NC_002947/3988855‑3989150                                                                                                                                                     |                                                                                                                                                     cggcggGCATGCACGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGa                                                                                                                                                      <  4:328362/148‑1 (MQ=255)    cggGCATGCAAGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGATGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTAGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGcc                                                                                                                                                  >  8:189516/1‑149 (MQ=255)       gCATGCCCGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACCCCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCCGCCTGTCGATGCGCCTGGACGCCggg                                                                                                                                               >  2:49085/1‑149 (MQ=255)       gCATGCACGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCggg                                                                                                                                               <  8:269268/149‑1 (MQ=255)            cACGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACg                                                                                                                                           <  5:178320/148‑1 (MQ=255)                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:121480/147‑2 (MQ=38)                           gCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCAGGTCGACGTGCCCGGCGTCCAGGCGCat                                                                                                                             <  5:222669/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                                                                                      >  4:221771/1‑141 (MQ=255)                                                                        tCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCtt                                                                                      <  3:221771/141‑1 (MQ=255)                                                                            gCTTGTGTGCGTTAGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCTCCGAACTACGTCAGCAAGGCCTTGCCAACCTgcgc                                                                          <  7:189516/149‑1 (MQ=255)                                                                              ttGTGTGCGTTCGGCCAGATCCTGTCCGACGAGCGCTA‑CGACCAGCTCACCAGCCTGTCCATCCGCCTGGACGCCGGGCACTACGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACCTCAGCAAGGCCTTGCCAACCTccgcc                                                                         >  2:79556/1‑147 (MQ=255)                                                                                   tgCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTg                                                                                     >  7:79565/1‑131 (MQ=255)                                                                                   tgCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTg                                                                                     <  8:79565/131‑1 (MQ=255)                                                                                     cGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGc                                                                 >  1:270341/1‑149 (MQ=255)                                                                                      gTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCt                                                                >  8:37334/1‑149 (MQ=255)                                                                                                 tCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAg                                                     <  5:260114/149‑1 (MQ=255)                                                                                                 tCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAg                                                     <  7:37334/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                                                                       >  6:253101/1‑100 (MQ=255)                                                                                                                                 cAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGaa                                                                       <  5:253101/100‑1 (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            <  6:162554/149‑1 (MQ=255)                                                                                                                                           aTGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTggg            >  4:292329/1‑149 (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)                                                                                                                                                     |                                                                                                                                                     CGGCGGGCATGCACGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGATTCAT  >  NC_002947/3988855‑3989150

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