Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,098,376	A→C	5.9%	K247T (AAG→ACG)	PP_2707 →	exodeoxyribonuclease III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,098,376	0	A	C	5.9%	218.1 / 3.2	68	K247T (AAG→ACG)	PP_2707	exodeoxyribonuclease III
Reads supporting (aligned to +/- strand):  ref base A (34/30);  new base C (2/2);  total (36/32)
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.93e-01

GCGTGTCTACACCTTCTGGGACTATTTTCGTAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGACTTCAGATTGCATTCAGATTCGGCCCCGATAATTCCCT  >  NC_002947/3098233‑3098524                                                                                                                                                |                                                                                                                                                      gCGTGTCTACACCTTCTGGGACTATTTTCGTAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCaggccaaggcca                                                                                                                                                  <  3:284625/149‑1 (MQ=255) gCGTGTCTACACCTTCTGGGACTATTTTCGTAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGAGTGTCGATGCATGGGTGCGCAACCaggccaaggc                                                                                                                                                    >  8:67312/1‑147 (MQ=255)           acCTTCTGGGACTATTTTCGTAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCACGGCCAGCGACCATgc                                                                                                                                        >  6:78721/1‑149 (MQ=255)           acCTTCTGGGACTATTTTCGTAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATgc                                                                                                                                        >  8:35665/1‑149 (MQ=255)                     aCTATTTTCGTAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTgg                                                                                                                             <  6:9165/149‑1 (MQ=255)                     aCTATTTTCGTAAACACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACAATGCCAAGGCAAGCGACCATGCG‑CCCGTCTg                                                                                                                              >  8:31451/1‑148 (MQ=255)                       tATTTTCGTAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGAt                                                                                                                           >  8:131733/1‑149 (MQ=255)                         ttttCGTAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATc                                                                                                                          >  4:113325/1‑148 (MQ=255)                              gTAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTg                                                                                                                    >  3:286284/1‑149 (MQ=255)                               tAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTgg                                                                                                                   >  4:249187/1‑149 (MQ=255)                               tAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTgg                                                                                                                   <  3:113325/149‑1 (MQ=255)                                     cTGGGCCCGCAATGCCTGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCt                                                                                                                     <  8:263756/141‑1 (MQ=255)                                     cTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCt                                                                                                                     >  7:263756/1‑141 (MQ=255)                                     cTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTAcgc                                                                                                             <  7:131733/149‑1 (MQ=255)                                            cGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGc                                                                                                      <  4:286284/149‑1 (MQ=255)                                               aaTGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCg                                                                                                   <  3:249187/149‑1 (MQ=255)                                                 tGCCGGGTTGCGAATCGATCACGTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGt                                                                                                 <  4:39559/149‑1 (MQ=255)                                                    cGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGTCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCAt                                                                                              <  5:149162/149‑1 (MQ=255)                                                     gggTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATa                                                                                             <  8:256590/149‑1 (MQ=255)                                                            gAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCACGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGtt                                                                                      <  5:78721/149‑1 (MQ=255)                                                               tCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAg                                                                                   <  7:142305/149‑1 (MQ=255)                                                                 gATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGt                                                                                       >  3:176922/1‑143 (MQ=255)                                                                 gATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGt                                                                                       <  4:176922/143‑1 (MQ=255)                                                                  aTCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTtcatca                                                                                 <  3:232612/148‑1 (MQ=255)                                                                   tCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGcaca                                                                               >  2:159587/1‑149 (MQ=255)                                                                   tCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGcaca                                                                               >  6:114083/1‑149 (MQ=255)                                                                   tCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGcaca                                                                               >  5:158730/1‑149 (MQ=255)                                                                     aCCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACAt                                                                              <  6:158730/148‑1 (MQ=255)                                                                       cTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGcgtc                                                                                                          <  1:112789/118‑1 (MQ=255)                                                                       cTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGcgtc                                                                                                          >  2:112789/1‑118 (MQ=255)                                                                       cTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGaa                                                                           >  6:287601/1‑149 (MQ=255)                                                                        tGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCa                                                                                         <  6:259998/134‑1 (MQ=255)                                                                        tGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCa                                                                                         >  5:259998/1‑134 (MQ=255)                                                                        tGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAAc                                                                          >  1:14540/1‑149 (MQ=255)                                                                              tCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCg                                                                                                   <  4:230997/118‑1 (MQ=255)                                                                              tCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCg                                                                                                   >  3:230997/1‑118 (MQ=255)                                                                              tCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAgcgc                                                                    >  4:165343/1‑149 (MQ=255)                                                                               cAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGGATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGcc                                                                   >  1:64899/1‑149 (MQ=255)                                                                                aaTCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGcc                                                                   <  1:159587/148‑1 (MQ=255)                                                                                    gcgAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCt                                                              <  2:64899/149‑1 (MQ=255)                                                                                    gcgAGCTTGCGCCTTATCTGAAACGGGCGAGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGc                                                               <  7:67312/148‑1 (MQ=255)                                                                                     cgAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCACGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCt                                                              >  6:239537/1‑148 (MQ=255)                                                                                          ttGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCgg                                                                 <  3:133772/140‑1 (MQ=255)                                                                                          ttGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCgg                                                                 >  4:133772/1‑140 (MQ=255)                                                                                            gcgcCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCt                                                      >  2:105259/1‑149 (MQ=255)                                                                                                attATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCATGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGa                                                  >  6:163526/2‑149 (MQ=255)                                                                                                   aTCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCACGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACa                                                <  5:239537/148‑1 (MQ=255)                                                                                                       gAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCggg                                                                >  6:132209/1‑128 (MQ=255)                                                                                                       gAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCggg                                                                <  5:132209/128‑1 (MQ=255)                                                                                                        aaaCGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGcc                                          <  3:154071/149‑1 (MQ=255)                                                                                                          aCGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGa                                        >  6:325438/1‑149 (MQ=255)                                                                                                              gcgggTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTc                                                     >  3:23075/1‑132 (MQ=255)                                                                                                              gcgggTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTc                                                     <  4:23075/132‑1 (MQ=255)                                                                                                                  gtgtCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGACTTCAGAt                                >  1:270712/1‑149 (MQ=255)                                                                                                                  gtgtCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGAACATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGTGCTGAGGCGCTCCGACAGTAGCCGACTTCAGAt                                <  7:31451/149‑1 (MQ=255)                                                                                                                        aTGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGACTTCAGATTGCAt                           >  1:65656/1‑149 (MQ=255)                                                                                                                         tGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGACTTCAGATTGCATTc                         <  7:275031/149‑1 (MQ=255)                                                                                                                             tGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGACTTCAGATTGCATTCAGAt                     >  1:111018/1‑149 (MQ=255)                                                                                                                               ggTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGACTTCAGATTGCATTCAGATTc                   <  3:165343/149‑1 (MQ=255)                                                                                                                                gTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGACTTCAGATTGCATTCAGATTCg                  >  2:158690/1‑149 (MQ=255)                                                                                                                                 tGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGa                                                                        >  1:26862/1‑94 (MQ=255)                                                                                                                                 tGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGa                                                                        <  2:26862/94‑1 (MQ=255)                                                                                                                                  gcgcAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGACTTCAGATTGCATTCAGATTCGGc                >  3:72252/1‑149 (MQ=255)                                                                                                                                   cgcAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGACTTCa                                   >  5:130579/1‑129 (MQ=255)                                                                                                                                   cgcAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGACTTCa                                   <  6:130579/129‑1 (MQ=255)                                                                                                                                    gcAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGACTTCAGATTGCATTCAGATTCGGccc              <  4:72252/149‑1 (MQ=255)                                                                                                                                    gcAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGACTTCAGATTGCATTCAGATTCGGccc              >  2:231252/1‑149 (MQ=255)                                                                                                                                       aCCAGGCCAAGGCCAGCGCCCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGACTTCAGATTGCATTCAGATTCGGCCCCGa           >  4:178282/1‑149 (MQ=255)                                                                                                                                        ccAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGtt                                                                                      <  6:238790/73‑1 (MQ=255)                                                                                                                                        ccAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGtt                                                                                      >  5:238790/1‑73 (MQ=255)                                                                                                                                              caaggccaGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCATGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGACTTCAGATTGCATTCAGATTCGGCCCCGATAATTcc    <  5:163526/149‑1 (MQ=255)                                                                                                                                                aggccaGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGACTTCAGATTGCATTCAGATTCGGCCCCGATAATTCCCt  <  5:287601/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                      GCGTGTCTACACCTTCTGGGACTATTTTCGTAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCG‑CCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGACTTCAGATTGCATTCAGATTCGGCCCCGATAATTCCCT  >  NC_002947/3098233‑3098524

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