BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,405,348	G→C	5.1%	W115C (TGG→TGC)	PP_2108 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,405,348	0	G	C	5.1%	379.0 / 7.7	117	W115C (TGG→TGC)	PP_2108	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (52/59); new base C (4/2); total (56/61)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.12e-01

Only 100 of 121 total aligned reads displayed.

GTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTC‑ACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGGCAGCAAGGGGATGTTCGTATGCTGTATATAAATACAGTGTCGAGGT  >  NC_002947/2405201‑2405494

gTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCGACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATgg                                                                                                                                                     >  1:91127/1‑149 (MQ=255)

 ttGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCTGCGGACTTC‑ACGGTGATCGACGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGt                                                                                                                                                   <  5:248772/149‑1 (MQ=255)

 ttGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTC‑ACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGt                                                                                                                                                   >  6:277805/1‑149 (MQ=255)

 ttGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTC‑ACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGt                                                                                                                                                   <  6:141787/149‑1 (MQ=255)

 ttGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTC‑ACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGt                                                                                                                                                   <  5:253268/149‑1 (MQ=255)

 ttGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTC‑ACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGt                                                                                                                                                   <  5:69498/149‑1 (MQ=255)

 gctaTCGGAGGCGCCACCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTC‑ACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGt                                                                                                                                                   >  8:229918/4‑149 (MQ=255)

   gATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTC‑ACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTg                                                                                                                                                  >  2:218306/1‑148 (MQ=255)

         tGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTC‑ACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGt                                                                                                                                           <  5:259931/149‑1 (MQ=255)

         tGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTC‑ACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGt                                                                                                                                           >  1:192815/1‑149 (MQ=255)

           gcgcCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTC‑ACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCg                                                                                                                                         <  5:277805/149‑1 (MQ=255)

                tCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTC‑ACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATgg                                                                                                                                                     >  8:62630/1‑132 (MQ=255)

                  cGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTC‑ACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGCACCCATCTGTCgaattcc                                                                                                                                   <  1:294365/148‑7 (MQ=255)

                       gCCAGCACGTGCAAGGGCTGGCGGCGGACTTC‑ACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGc                                                                                                                             <  1:200046/149‑1 (MQ=255)

                        ccAGCACGTGCAAGGGCTGGCGGCGGACTTC‑ACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCa                                                                                                                            >  2:88968/1‑149 (MQ=255)

                         cAGCACGTGCAAGGGCTGGCGGCGGACTTC‑ACGGTGATCGAGGTCAGCCCGCGCGAGCCCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCAc                                                                                                                           >  8:42168/1‑149 (MQ=255)

                          aGCACGTGCAAGGGCTGGCGGCGGACTTC‑ACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCAcg                                                                                                                          >  5:96048/1‑149 (MQ=255)

                             aCGTGCAAGGGCTGGCGGCGGACTTC‑ACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCAcgccg                                                                                                                       <  4:115473/149‑1 (MQ=255)

                             aCGTGCAAGGGCTGGCGGCGGACTTC‑ACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCAcgccg                                                                                                                       >  3:11878/1‑149 (MQ=255)

                                 gCAAGGGCTGGCGGCGGACTTC‑ACGGTGATCGAGGTCAGCCCGGGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCgca                                                                                                                   <  4:81392/149‑1 (MQ=255)

                                 gCAAGGGCTGGCGGCGGACTTC‑ACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCcgcg                                                                                                                     >  4:95/1‑147 (MQ=255)

                                  cAAGGGCTGGCGGCGGACTTC‑ACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCgcag                                                                                                                  <  4:6237/149‑1 (MQ=255)

                                   aaGGGCTGTCGGCGTAATTC‑ACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCgcagg                                                                                                                 <  8:257155/149‑1 (MQ=255)

                                         tGGCGGCGGACTTC‑ACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTg                                                                                                           >  1:172647/1‑149 (MQ=255)

                                           gcggcgGACTTC‑ACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGATTCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCt                                                                                                         >  8:140519/1‑149 (MQ=255)

                                           gcggcgGACTTC‑ACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTg                                                                                                           <  2:253324/147‑1 (MQ=255)

                                                 gACTTC‑ACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCAt                                                                                                   <  2:172647/149‑1 (MQ=255)

                                                    ttC‑ACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATgcg                                                                                                >  3:86116/1‑149 (MQ=255)

                                                    ttC‑ACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATgcg                                                                                                <  8:242756/149‑1 (MQ=255)

                                                         cGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCaaa                                                                                            <  3:221719/149‑1 (MQ=255)

                                                          ggTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGTTGCCCTTCGATCAGCGGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGAGCAAAg                                                                                           >  5:108891/1‑149 (MQ=255)

                                                           gTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAg‑g                                                                                         <  4:11878/149‑1 (MQ=255)

                                                           gTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAg‑g                                                                                         <  7:229918/149‑1 (MQ=255)

                                                              aTCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAg                                                                                      >  3:147142/1‑149 (MQ=255)

                                                                cGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCgcaggcag                                                                                                              <  6:135128/124‑1 (MQ=255)

                                                                 gAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGCCCATGCGCAAAG‑GCAGCg                                                                                    >  6:230991/1‑148 (MQ=255)

                                                                 gAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGccc                                                                                   >  7:56536/1‑147 (MQ=255)

                                                                 gAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCgg                                                                                   >  3:116709/1‑149 (MQ=255)

                                                                      cctatcgcgcgAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTAcc                                                                              <  1:281989/144‑1 (MQ=255)

                                                                          ccGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGc                                                                          >  6:145274/1‑149 (MQ=255)

                                                                           cgcgcgAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGc                                                                                     <  7:60871/137‑1 (MQ=255)

                                                                           cgcgcgAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGc                                                                                     >  8:60871/1‑137 (MQ=255)

                                                                            gcgcgAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCgg                                                                        <  1:218306/149‑1 (MQ=255)

                                                                              gcgAGACCGTCCGGCGAATCAGCAAAAGTTCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAg                                                                                           >  7:270009/1‑129 (MQ=255)

                                                                              gcgAGACCGTCCGGCGAATCAGCAAAAGTTCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAg                                                                                           <  8:270009/129‑1 (MQ=255)

                                                                              gcgAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCgggg                                                                      <  4:147142/149‑1 (MQ=255)

                                                                                         cGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAg                                                           <  5:145274/149‑1 (MQ=255)

                                                                                          ggCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCg                                                                         <  6:246956/134‑1 (MQ=255)

                                                                                           gCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAgcg                                                         <  1:33763/149‑1 (MQ=255)

                                                                                           gCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAgcg                                                         >  3:28239/1‑149 (MQ=255)

                                                                                           gCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAgcg                                                         >  3:24099/1‑149 (MQ=255)

                                                                                             gAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCa                                                       >  4:258396/1‑149 (MQ=255)

                                                                                                 cAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGt                                                                                                                                       <  1:188835/66‑1 (MQ=255)

                                                                                                 cAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGt                                                                                                                                       >  2:188835/1‑66 (MQ=255)

                                                                                                 cAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGa                                                   >  4:282991/1‑149 (MQ=255)

                                                                                                 cAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGa                                                   >  6:163024/1‑149 (MQ=255)

                                                                                                   gCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACg                                                 >  8:169726/1‑149 (MQ=255)

                                                                                                   gCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACg                                                 >  7:181173/1‑149 (MQ=255)

                                                                                                    cAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACgg                                                <  5:299706/149‑1 (MQ=255)

                                                                                                     aaaaGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACgg                                                >  5:253988/1‑148 (MQ=255)

                                                                                                      aaaGTGCGGTGCCCTTCGATCAGCTGATGCTGTAATTCGACAGATGGCTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGGTGACCATGCGCAAAG‑GGAGCGGCTACCTGGCGTGGCTGCAATGAAGCTCATTGAc                                                  <  7:84801/145‑1 (MQ=255)

                                                                                                      aaaGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGgca                                              <  5:163024/149‑1 (MQ=255)

                                                                                                      aaaGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGgca                                              <  3:282991/149‑1 (MQ=255)

                                                                                                      aaaGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGgca                                              <  4:116709/149‑1 (MQ=255)

                                                                                                      aaaGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGAATCGGTCGGTGACGCGCGGCACGCCCCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGATGCAATGAAGCGAATTGACgag                                               >  8:84801/1‑146 (MQ=255)

                                                                                                       aaGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGgcag                                             <  7:42168/149‑1 (MQ=255)

                                                                                                       aaGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGgc                                               <  1:316988/147‑1 (MQ=255)

                                                                                                       aaGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGgc                                               >  2:316988/1‑147 (MQ=255)

                                                                                                           gCGGAGCCCTTCGATCAGCTGAACCTGGAATTCGACAGATCGGTGCATCTGTCGGAGACGCGCGGCACGCCGCGCACGCAGGTGCTGACCATGCGCAAAG‑GCAGAGGCTACCTCGCGCGGCTGCAATGAAACGCATTGACGGCAGCAAg                                         >  6:270134/1‑149 (MQ=255)

                                                                                                            cGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGGTAGCAAgg                                        >  6:253139/1‑149 (MQ=255)

                                                                                                            cGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGGCAGCAAgg                                        <  6:96048/149‑1 (MQ=255)

                                                                                                            cGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGGCAGCAAg                                         <  7:169726/148‑1 (MQ=255)

                                                                                                              gTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGGCAGCAAgggg                                      <  4:55561/149‑1 (MQ=255)

                                                                                                                 cccTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGGCAGCAAGGGGAt                                    <  4:24099/148‑1 (MQ=255)

                                                                                                                  ccTTCGATCAGCTGATCCTGGAAATCGACAGCTGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAATCGCATTGACGGCAGCAAGGGGATGt                                  <  4:55596/149‑1 (MQ=255)

                                                                                                                   cTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGGCAGCAAGGGGATGtt                                 >  7:207830/1‑149 (MQ=255)

                                                                                                                   cTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGGCAGCAAGGGGATGtt                                 >  8:264961/1‑149 (MQ=255)

                                                                                                                   cTTCGATCAGCTGATCCTGCAATTCGTCAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGTGGCTGCAATGAAGCGCATTGACGGCAGCAAGGGGATGtt                                 <  2:239812/149‑1 (MQ=255)

                                                                                                                    ttCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCgggg                                                                      <  7:30322/111‑1 (MQ=255)

                                                                                                                      cGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGGCAGCAAGGGGATGTTc                                <  2:242869/147‑1 (MQ=255)

                                                                                                                      cGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGGCAGCAAGGGGATGTTc                                <  2:91127/147‑1 (MQ=255)

                                                                                                                       gATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGGCAGCAAGGGGATGTTCGTa                             <  3:258396/149‑1 (MQ=255)

                                                                                                                         tCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGGCAGCAAGGGGATGTTCGTATg                           >  2:58518/1‑149 (MQ=255)

                                                                                                                            gCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGGCAGCAAGGGGATGTTCGTATGCTg                        >  8:144323/1‑149 (MQ=255)

                                                                                                                                aTCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGGCAGCAAGGGGATGTTCGTATGCTGtata                    >  5:62079/1‑149 (MQ=255)

                                                                                                                                   cTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGGCAGCAAGGGGATGTTCGTATGCTGtatata                  >  5:183629/1‑148 (MQ=255)

                                                                                                                                     ggAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGGCAGCAAGGGGATGTTCGTATGCTg                        >  3:136046/1‑140 (MQ=255)

                                                                                                                                     ggAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGGCAGCAAGGGGATGTTCGTATGCTg                        <  4:136046/140‑1 (MQ=255)

                                                                                                                                         ttCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGGCAGCAAGGGGATGTTCGTATGCTGTATATAAATACAg           >  3:290073/1‑149 (MQ=255)

                                                                                                                                          ccGACAGATGCACCCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGGCAGCAAGGGGATGTTCGTATGCTGTATATAAATACAgt          >  3:104898/2‑149 (MQ=255)

                                                                                                                                           cGACAGATGCACCCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCt                                                                    >  8:169269/1‑90 (MQ=255)

                                                                                                                                           cGACAGATGCACCCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCt                                                                    <  7:169269/90‑1 (MQ=255)

                                                                                                                                           cGACAGATGCACCCATATGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGGCAGCAAGGGGATGTTCGTATGCTGTATATAAATACAgtg         >  1:305360/1‑149 (MQ=255)

                                                                                                                                            gACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGGCAGCAAGGGGATGTTCGTATGCTGTATATAAATACAgtgt        >  2:149796/1‑149 (MQ=255)

                                                                                                                                             aCAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGc                                                                                                          >  4:2345/1‑51 (MQ=255)

                                                                                                                                                gATGGGTTCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGGCAGCAAGGGGATGTTCGTATGCTGTATATAAATACAGTGTCGAg    <  7:140519/149‑1 (MQ=255)

                                                                                                                                                gATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGGCAGCAAGGGGATGTTCGTATGCTGTATATAAATACAGTGTCGAg    <  6:62079/149‑1 (MQ=255)

                                                                                                                                                gATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGGCAGCAAGGGGATGTTCGTATGCTGTATATAAATACAGTGTCGAg    <  6:183629/149‑1 (MQ=255)

                                                                                                                                                  tGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTTCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGGCAGCAAGGGGATGTTCGTATGCTGTATATAAATACAGTGTCGAGGt  <  7:227078/149‑1 (MQ=255)

                                                                                                                                                    ggTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAGAGAAGCGGATACATGGCGGGCCTGCAATGAAGCGCATTGACGGAACCAAGGGCATTTTCGTATGCTTTATATAAAAAAAGTTACGAGGt  >  1:95254/1‑148 (MQ=255)

GTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTC‑ACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCACGCCGCGCAGGCAGGTGCTGACCATGCGCAAAG‑GCAGCGGCTACCTGGCGGGGCTGCAATGAAGCGCATTGACGGCAGCAAGGGGATGTTCGTATGCTGTATATAAATACAGTGTCGAGGT  >  NC_002947/2405201‑2405494

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