BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	31,611	G→A	6.5%	V70I (GTC→ATC) ‡	PP_0027 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	31,611	0	G	A	6.5%	313.2 / 7.1	107	V70I (GTC→ATC) ‡	PP_0027	membrane protein
Reads supporting (aligned to +/- strand): ref base G (49/50); new base A (3/4); total (52/55)
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.79e-01

Only 100 of 110 total aligned reads displayed.

GCCCTCGTCAACAACCTCCAGAAGGCTCAGCGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGATGTTTATTGCTTTTCCATTGCTGCTACCTAAGCTCGGCTTTTG  >  NC_002947/31462‑31753

gCCCTCGTCAACAACCTCCAGAAGGCTCAGCGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑Cgg                                                                                                                                                <  1:1937555/150‑1 (MQ=255)

 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                                                                                                                               <  1:885067/149‑1 (MQ=255)

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

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

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

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

                          tCAGCGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGACCGTGAccaa                                                                                                                                      <  1:961245/134‑4 (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)

                                          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)

                                              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                                                                                                       >  1:433625/1‑129 (MQ=255)

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

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

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

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

                                                                  ccGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTg                                                                                <  1:529284/148‑1 (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                                                                           <  2:1901519/146‑1 (MQ=255)

                                                                          gTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGACGATCTCCAACCATGCTTGGTACACGTTc                                                                      >  2:495547/1‑150 (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)

                                                                                  cGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGt                                                              <  2:1819409/150‑1 (MQ=255)

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

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

                                                                                        agagCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCa                                                                                                              <  1:106195/96‑1 (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)

                                                                                                  aaaCGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCg                                              >  1:2056472/1‑150 (MQ=255)

                                                                                                    aCGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGAt                                            <  2:449909/150‑1 (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)

                                                                                                              aaGCTGGGCGTCCTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGTTGTTTATTGCt                                  >  2:2342363/1‑150 (MQ=255)

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

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

                                                                                                                gCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGATGTTTATTGCttt                                <  2:95569/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)

                                                                                                                              tCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGATGTTTATTGCTTTTCCATTGCTGCTAc                  >  1:936634/1‑150 (MQ=255)

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

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

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

                                                                                                                                       tACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGATGTTTATTGCTTTTCCATTGCTGCTACCTAAGCTCg         <  2:68681/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                                     >  2:1220968/1‑113 (MQ=255)

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

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

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

GCCCTCGTCAACAACCTCCAGAAGGCTCAGCGATGGCCTGGATCATCACGAAGTACCTACTCACTGCCGGTATAGTGATTCTCGTGTCAGAGCTGGCGAAACGCAGCGACAAGCTGGGCGGCTTGGTCGCCGCCCTACCCTTGGTCA‑CGGTCGTGACGTTGGTCTGGCTCTATTTAGAGAACCAGAGCATGGAGAAGATCTCCAACCATGCTTGGTACACGTTCTGGTACGTGCTTCCGACCCTGCCGATGTTTATTGCTTTTCCATTGCTGCTACCTAAGCTCGGCTTTTG  >  NC_002947/31462‑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: `‑`