BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	31,612	T→C	7.3%	V70A (GTC→GCC) ‡	PP_0027 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	31,612	0	T	C	7.3%	314.4 / 7.9	108	V70A (GTC→GCC) ‡	PP_0027	membrane protein
Reads supporting (aligned to +/- strand): ref base T (50/50); new base C (3/5); total (53/55)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 109 total aligned reads displayed.

CCCTCGTCAACAACCTCCAGAAGGCTCAGCGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGATGTTTATTGCTTTTCCATTGCTGCTACCTAAGCTCGGCTTTTG  >  NC_002947/31463‑31753

cccTCGTCAACAACCTCCAGAAGGCTCAGCGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGt                                                                                                                                               <  1:1048975/150‑1 (MQ=255)

 ccTCGTCAACAACCTCCAGAAGGCTCAGCGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTc                                                                                                                                              >  1:1100602/1‑150 (MQ=255)

  cTCGTCAACAACCTCCAGAAGGCTCAGCGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCg                                                                                                                                             >  2:2095987/1‑150 (MQ=255)

   tCGTCAACAACCTCCAGAAGGCTCAGCGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGt                                                                                                                                            >  2:421456/1‑150 (MQ=255)

   tCGTCAACAACCTCCAGAAGGCTCAGCGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGt                                                                                                                                            >  1:449727/1‑150 (MQ=255)

     gTCAACATCTTCTAGAAGGCTCAGCTATAGCCTGGTTCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGa                                                                                                                                          <  1:100421/150‑1 (MQ=255)

        aacaacCTCCAGAAGGCTCAGCGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGt                                                                                                                                       >  2:105263/1‑150 (MQ=255)

         acaacCTCCAGAAGGCTCAGCGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGtt                                                                                                                                      >  1:1539769/1‑150 (MQ=255)

         acaacCTCCAGAAGGCTCAGCGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGACCGTGAccaa                                                                                                                                      <  1:1171763/150‑4 (MQ=255)

          caacCTCCAGAAGGCTCAGCGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTg                                                                                                                                     >  1:95569/1‑150 (MQ=255)

          caacCTCCAGAAGGCTCAGCGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTAGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTg                                                                                                                                     <  1:194677/150‑1 (MQ=255)

            acCTCCAGAAGGCTCAGCGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGt                                                                                                                                   <  2:2115895/150‑1 (MQ=255)

             ccTCCAGAAGGCTCAGCGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGt                                                                                                                                   <  1:1670666/149‑1 (MQ=255)

             ccTCCAGAAGGCTCAGCGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGt                                                                                                                                   >  2:1670666/1‑149 (MQ=255)

             ccTCCAGAAGGCTCAGCGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTc                                                                                                                                  <  2:449727/150‑1 (MQ=255)

               tCCAGAAGGCTCAGCGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTGGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTg                                                                                                                                >  1:38783/1‑150 (MQ=255)

               tCCAGAAGGCTCAGCGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTg                                                                                                                                >  1:1819409/1‑150 (MQ=255)

                 cAGAAGGCTCAGCGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTgg                                                                                                                               >  2:885067/1‑149 (MQ=255)

                 cAGAAGGCTCAGCGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTgg                                                                                                                               <  1:885067/149‑1 (MQ=255)

                 cAGAAGGCTCAGCGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGc                                                                                                                              <  1:2233025/150‑1 (MQ=255)

                         tCAGCGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTa                                                                                                                      >  1:936570/1‑150 (MQ=255)

                         tCAGCGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGACCGTGAccaa                                                                                                                                      <  1:961245/134‑4 (MQ=255)

                         tCAGCGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGACCGTGAccaa                                                                                                                                      >  2:961245/1‑131 (MQ=255)

                          cAGCGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTag                                                                                                                     <  2:322663/150‑1 (MQ=255)

                             cGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAgaga                                                                                                                  >  2:1271985/1‑150 (MQ=255)

                              gATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGaa                                                                                                                 >  2:2349534/1‑150 (MQ=255)

                               aTGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAAc                                                                                                                <  2:1653549/150‑1 (MQ=255)

                                     tGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGc                                                                                                          <  2:1100602/150‑1 (MQ=255)

                                        atcatcACGAAGTACCTACTCAATGCCGTTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATg                                                                                                       <  1:978478/150‑1 (MQ=255)

                                         tcatcaCGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATgg                                                                                                      <  1:105263/150‑1 (MQ=255)

                                            tcaCGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTGTGGCTCTATTTAGAGAACCAGAGCATGgaga                                                                                                   <  2:1486043/150‑1 (MQ=255)

                                            tcaCGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCt                                                                                                                                 <  2:1055801/120‑1 (MQ=255)

                                            tcaCGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCt                                                                                                                                 >  1:1055801/1‑120 (MQ=255)

                                            tcaCGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGgaga                                                                                                   >  1:1516952/1‑150 (MQ=255)

                                             caCGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTCTTTAGAGAACCAGAGCATGGagaa                                                                                                  >  2:316203/1‑150 (MQ=255)

                                              aCGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGagaag                                                                                                 >  1:1762766/1‑150 (MQ=255)

                                                   gTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCAGCGTCCTTGTTGTT‑GTGGGGCTCTATTTAGAGAACCAGAGCATGGAGAAGAtctc                                                                                            <  1:421456/150‑1 (MQ=255)

                                                    tACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTcc                                                                                           <  2:1516952/150‑1 (MQ=255)

                                                      ccTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCaa                                                                                         <  1:2095987/150‑1 (MQ=255)

                                                      ccTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGACGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTAGTGACGTTGGTCTGGCTCGATTTAGAGAACCAGAGCATGGAGAAGATCTCCaa                                                                                         >  1:2300737/1‑150 (MQ=255)

                                                       cTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAAc                                                                                        <  1:2173159/150‑1 (MQ=255)

                                                          cTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCAt                                                                                     <  1:1879547/150‑1 (MQ=255)

                                                            cACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGACCTCCAACCATGc                                                                                   <  2:4532/150‑1 (MQ=255)

                                                             aCTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATg                                                                                                       <  2:433625/129‑1 (MQ=255)

                                                             aCTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCt                                                                                  >  1:1928222/1‑150 (MQ=255)

                                                                gCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTgg                                                                               >  1:1779732/1‑150 (MQ=255)

                                                                gCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTgg                                                                               <  1:345709/150‑1 (MQ=255)

                                                                 ccGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTg                                                                                <  1:529284/148‑1 (MQ=255)

                                                                 ccGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATGTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTg                                                                                >  2:529284/1‑148 (MQ=255)

                                                                 ccGGTATAGTGAGTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAACCTGGGCTGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGGATGGAGAAGATCTCCAACCATGCTTGGt                                                                              <  2:2300737/150‑1 (MQ=255)

                                                                    gTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTaca                                                                           >  2:2233521/1‑150 (MQ=255)

                                                                       taGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGAtc                                                                                              <  2:92742/128‑1 (MQ=255)

                                                                       taGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGAtc                                                                                              >  1:92742/1‑128 (MQ=255)

                                                                        aGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTaca                                                                           >  1:1901519/1‑146 (MQ=255)

                                                                        aGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTaca                                                                           <  2:1901519/146‑1 (MQ=255)

                                                                          tGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCt                                                                     <  2:1539769/150‑1 (MQ=255)

                                                                          tGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCt                                                                     >  1:68681/1‑150 (MQ=255)

                                                                             ttCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGt                                                                  >  2:2274736/1‑150 (MQ=255)

                                                                                tcGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACg                                                               >  2:1814043/1‑150 (MQ=255)

                                                                                  gtgtCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGATCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTg                                                             >  1:1898513/1‑150 (MQ=255)

                                                                                       agagCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCa                                                                                                              <  1:106195/96‑1 (MQ=255)

                                                                                       agagCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCa                                                                                                              >  2:106195/1‑96 (MQ=255)

                                                                                        gagCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCg                                                       >  2:2361833/1‑150 (MQ=255)

                                                                                              gCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTgg                                                                               <  2:583659/120‑1 (MQ=255)

                                                                                              gCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTgg                                                                               >  1:583659/1‑120 (MQ=255)

                                                                                               cGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACAATGTTTGGTACACGTTCTGGTACGAGATTCCGACACTGc                                                >  2:797706/1‑150 (MQ=255)

                                                                                                   aCGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGAt                                            <  2:449909/150‑1 (MQ=255)

                                                                                                   aCGCAGCGACAAGCTGGGCGGCTTGGTCGCCCCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACAAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGAt                                            <  1:797706/150‑1 (MQ=255)

                                                                                                     gCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGATGt                                          >  2:2007162/1‑150 (MQ=255)

                                                                                                        gCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGtt                                                                       >  2:842454/1‑118 (MQ=255)

                                                                                                        gCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGtt                                                                       <  1:842454/118‑1 (MQ=255)

                                                                                                        gCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGATGTTTa                                       >  1:938465/1‑150 (MQ=255)

                                                                                                         cGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGATGTTTAt                                      <  2:936570/150‑1 (MQ=255)

                                                                                                          gACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTATGTGCTTCCGACCCTGCCGATGTTTAtt                                     <  1:1271985/150‑1 (MQ=255)

                                                                                                          gACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGATGTTTAtt                                     <  2:1653925/150‑1 (MQ=255)

                                                                                                             aaGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGATGTTTATTGCt                                  >  2:1545681/1‑150 (MQ=255)

                                                                                                               gCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGATGTTTATTGCttt                                <  2:95569/150‑1 (MQ=255)

                                                                                                               gCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGATGTTTATTGCttt                                <  1:2361833/150‑1 (MQ=255)

                                                                                                                cTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTACAACCATGCTTGCTACAAGTTCTGTTACGTGCTTCCGACCCTGCCGATGTTTATTGCtttt                               <  1:2349534/150‑1 (MQ=255)

                                                                                                                    gcggcTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGATGTTTATTGCTTTTCCAt                           <  2:2056472/150‑1 (MQ=255)

                                                                                                                     cggcTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGATGTTTATTGCTTTTCCAtt                          >  2:1789672/1‑150 (MQ=255)

                                                                                                                     cggcTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGATGTTTATTGCTTTTCCAtt                          <  2:1202988/150‑1 (MQ=255)

                                                                                                                         ttGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGtt                                                                       <  1:973047/101‑1 (MQ=255)

                                                                                                                         tgggTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGATGTTTATTGCTTTTCCATtgctg                      >  1:1157569/3‑150 (MQ=255)

                                                                                                                            gTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCagag                                                                                                           <  2:2323895/62‑1 (MQ=255)

                                                                                                                            gTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCagag                                                                                                           >  1:2323895/1‑62 (MQ=255)

                                                                                                                             tCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGATGTTTATTGCTTTTCCATTGCTGCTAc                  <  2:2164181/150‑1 (MQ=255)

                                                                                                                               gccgccCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGATGTTTATTGCTTTTCCATTGCTGCTACCt                >  2:1202869/1‑150 (MQ=255)

                                                                                                                               gccgccCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGATGTTTATTGCTTTTCCATTGCTGCTACCt                >  1:1743471/1‑150 (MQ=255)

                                                                                                                                      tACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGATGTTTATTGCTTTTCCATTGCTGCTACCTAAGCTCg         <  2:68681/150‑1 (MQ=255)

                                                                                                                                      tACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGATGTTTATTGCTTTTCCATTGCTGCTACCTAAGCTCg         <  2:1928222/150‑1 (MQ=255)

                                                                                                                                       aCCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGATGTTTATTGCTTTTCCATTGCTGCTACCTAAGCTCgg        <  2:581384/150‑1 (MQ=255)

                                                                                                                                          caacgTCA‑CGACCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTAc                                                                <  2:482506/87‑1 (MQ=255)

                                                                                                                                          caacgTCA‑CGACCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTAc                                                                >  1:482506/5‑91 (MQ=255)

                                                                                                                                           aacgTCA‑CGACCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTacac                                                                          <  1:900277/77‑1 (MQ=255)

                                                                                                                                           aacgTCA‑CGACCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGATGTTTATTGCTTTTCCATTGCTGCTACCTAAGCTCGGCttt    >  2:1565707/4‑150 (MQ=255)

                                                                                                                                             ggTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGATGTTTATTGCTTTTCCATTGCTGCTACCTAAGCTCGGCTTTTg  <  1:1814043/150‑1 (MQ=255)

                                                                                                                                               tCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGATGTTTAtt                                     <  1:1220968/113‑1 (MQ=255)

                                                                                                                                                    ggTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATgg                                                                                                      <  2:353331/44‑1 (MQ=255)

                                                                                                                                                    ggTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATgg                                                                                                      >  1:353331/1‑44 (MQ=255)

CCCTCGTCAACAACCTCCAGAAGGCTCAGCGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGATGTTTATTGCTTTTCCATTGCTGCTACCTAAGCTCGGCTTTTG  >  NC_002947/31463‑31753

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`