BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,446,589	G→A	5.8%	G314S (GGC→AGC)	gpD →	tail formation

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,446,589	0	G	A	5.8%	306.4 / 5.7	103	G314S (GGC→AGC)	gpD	tail formation
Reads supporting (aligned to +/- strand): ref base G (35/62); new base A (4/2); total (39/64)
Fisher's exact test for biased strand distribution p-value = 1.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 103 total aligned reads displayed.

ATCGCAGCACCGCCGGCGTGCGGCTGGAAATGCCCGGCCGCACCGACCTGTTTGCCGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCT  >  NC_002947/3446446‑3446735

aTCGCAGCACCGCCGGCGTGCGGCTGGAAATGCTCTGTTTCTCCGACCTGTTTGCCGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGa                                                                                                                                              <  1:1865781/150‑1 (MQ=255)

 tCGCAGCACCGCCGGCGTGCGGCTGGAAATGCCCGGCCGCACCGACCTGTTTGCCGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGAc                                                                                                                                             >  1:1583706/1‑150 (MQ=255)

   gcagcaCCGCCGGCGTGCGGCTGGAAATGCCCGGCCGCACCGACCTGTTTGCCGAACGGTCGCTCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCGGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGATGGACTa                                                                                                                                           >  2:85097/1‑150 (MQ=255)

    cagcaCCGCCGGCGTGCGGCTGGAAATGCCCGGCCGCACCGACCTGTTTGCCGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTAc                                                                                                                                          >  1:1730953/1‑150 (MQ=255)

      gcaCCGCCGGCGTGCGGCTGGAAATGCCCGGCCGCACCGACCTGTTTGCCGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGTTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGa                                                                                                                                        <  1:815905/150‑1 (MQ=255)

        aCCGCCGGCGTGCGGCTGGAAATGCCCGGCCGCACCGACCTGTTTGCCGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGAcc                                                                                                                                      <  2:564507/150‑1 (MQ=255)

             cgGCGTGCGGCTGGAAATGCCCGGCCGCACCGACCTGTTTGCCGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGa                                                                                                                                 <  2:920622/150‑1 (MQ=255)

             cgGCGTGCGGCTGGAAATGCCCGGCCGCACCGACCTGTTTGCCGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGa                                                                                                                                 <  1:39898/150‑1 (MQ=255)

             cgGCGTGCGGCTGGAAATGCCCGGCCGCACCGACCTGTTTGCCGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGa                                                                                                                                 <  2:1208841/150‑1 (MQ=255)

               gCGTGCGGCCGGAAATGCCCGGCCGCACCGACCTGTTTGCCGAACGCTCGAGCAACGCCCAAGGCTTCAAGCCGGGACTGGCTGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGt                                                                                                                               <  2:1325339/150‑1 (MQ=255)

                    cGGCTGGAAATGCCCGGCCGCACCGACCTGTTTGCCGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAAc                                                                                                                          <  1:614560/150‑1 (MQ=255)

                      gCTGGAAATGCCCGGCCGCACCGACCTGTTTGCCGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAAcgg                                                                                                                        >  1:620354/1‑150 (MQ=255)

                       cTGGAAATGCCCGGCCGCACCGACCTGTTTGCCGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAAcggc                                                                                                                       <  1:1779671/150‑1 (MQ=255)

                       cTGGAAATGCCCGGCCGCACCGACCTGTTTGCCGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAAcggc                                                                                                                       >  2:1422416/1‑150 (MQ=255)

                         ggAAATGCCCGGCCGCACCGACCTGTTTGCCGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAAcggcgg                                                                                                                     <  1:326194/150‑1 (MQ=255)

                          gAAATGCCCGGCCGCACCGACCTGTTTGCCGAACGCTCGATCAACGCCCCAGGCTTCCAGCCGGGGTTCCATGTTGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACggcggc                                                                                                                    <  1:1772325/150‑1 (MQ=255)

                          gAAATGCCCGGCCGCACCGACCTGTTTGCCGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACggcggc                                                                                                                    >  2:1379692/1‑150 (MQ=255)

                            aaTGCCCGGCCGCACCGACCTGTTTGCCGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCaa                                                                                                                  <  1:722755/150‑1 (MQ=255)

                              tGCCCGGCCGCACCGACCTGTTTGCCGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCaaga                                                                                                                <  1:970375/150‑1 (MQ=255)

                               gCCCGGCCGCACCGACCTGTTTGCCGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCaagaa                                                                                                               <  2:1080162/150‑1 (MQ=255)

                               gCCCGGCCGCACCGACCTGTTTGCCGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCaagaa                                                                                                               <  1:219995/150‑1 (MQ=255)

                                  cGGCCGCACCGACCTGTTTGCCGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAggg                                                                                                            <  1:1658771/150‑1 (MQ=255)

                                           cGACCTGTTTGCCGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGcaaggccaa                                                                                                   <  2:144465/150‑1 (MQ=255)

                                            gACCTGTTTGCCGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGcaaggccaag                                                                                                  <  1:471178/150‑1 (MQ=255)

                                             aCCTGTTTGCCGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGcaaggccaagg                                                                                                 <  2:1583706/150‑1 (MQ=255)

                                               cTGTTTGCCGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCaaggccaaggcc                                                                                               >  1:355645/1‑150 (MQ=255)

                                                   ttGCCGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCgg                                                                                           <  2:91483/150‑1 (MQ=255)

                                                   ttGCCGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCgg                                                                                           <  1:536950/150‑1 (MQ=255)

                                                      ccGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCgg                                                                                           >  2:58640/1‑147 (MQ=255)

                                                      ccGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCgg                                                                                           <  1:58640/147‑1 (MQ=255)

                                                      ccGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCa                                                                                        >  2:1648244/1‑150 (MQ=255)

                                                        gAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCaag                                                                                      >  2:1614089/1‑150 (MQ=255)

                                                         aaCGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGcaaggccaagg                                                                                                 <  1:344026/138‑1 (MQ=255)

                                                         aaCGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGcaaggccaagg                                                                                                 >  2:344026/1‑138 (MQ=255)

                                                            gCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCaagaaaa                                                                                  <  1:1422416/150‑1 (MQ=255)

                                                            gCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGGATTTACCCAGTCCGGTTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCaagaaaa                                                                                  >  2:460888/1‑150 (MQ=255)

                                                              tCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCaagaaaaag                                                                                >  1:1518446/1‑150 (MQ=255)

                                                              tCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCaagaaaaag                                                                                >  2:1718538/1‑150 (MQ=255)

                                                               cGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCaagaaaaaga                                                                               <  1:1312624/150‑1 (MQ=255)

                                                                gATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCaagaaaaagga                                                                              >  1:753367/1‑148 (MQ=255)

                                                                 aTCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCaagaaaaagaaa                                                                             <  1:158192/150‑1 (MQ=255)

                                                                   cAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGTTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAaagaaaga                                                                           <  1:460888/150‑1 (MQ=255)

                                                                    aaCGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTc                                                                                                                                   <  2:497015/93‑1 (MQ=255)

                                                                    aaCGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTc                                                                                                                                   >  1:497015/1‑93 (MQ=255)

                                                                     aCGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGa                                                                                                                                 <  1:668710/94‑1 (MQ=255)

                                                                     aCGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGa                                                                                                                                 >  2:668710/1‑94 (MQ=255)

                                                                     aCGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGaca                                                                         <  1:211343/150‑1 (MQ=255)

                                                                     aCGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCAGCCGGACtgg                                                                                                                                          <  2:1245776/85‑3 (MQ=255)

                                                                     aCGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCAGCCGGACtgg                                                                                                                                          >  1:1245776/1‑83 (MQ=255)

                                                                      cGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGacac                                                                        <  2:1162092/150‑1 (MQ=255)

                                                                         ccAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCaa                                                                     <  1:522421/150‑1 (MQ=255)

                                                                         ccAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCaa                                                                     >  1:932190/1‑150 (MQ=255)

                                                                         ccAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCaa                                                                     <  2:355645/150‑1 (MQ=255)

                                                                           aaGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGCCTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAACAGAAAGACACCAAGc                                                                   >  2:811311/1‑150 (MQ=255)

                                                                             ggCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAACAAAAAACAGCaccgag                                                                 >  2:1355518/1‑145 (MQ=255)

                                                                                ttCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCAGCTGGACTACGACCGTCGAGTGCAACGGCGGCaagaag                                                                                                              >  2:679309/1‑102 (MQ=255)

                                                                                ttCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCAGCTGGACTACGACCGTCGAGTGCAACGGCGGCaagaag                                                                                                              <  1:679309/102‑1 (MQ=255)

                                                                                ttCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGCATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTc                                                              >  2:1173549/1‑150 (MQ=255)

                                                                                    aGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAgg                                                          >  1:1258474/1‑150 (MQ=255)

                                                                                        ggggCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGa                                                      <  2:1825313/150‑1 (MQ=255)

                                                                                             tCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCaagaaaaagaaa                                                                             <  2:1666708/122‑1 (MQ=255)

                                                                                             tCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCaagaaaaagaaa                                                                             >  1:1666708/1‑122 (MQ=255)

                                                                                              cGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCt                                                <  2:260821/150‑1 (MQ=255)

                                                                                              cGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCt                                                <  2:1258474/150‑1 (MQ=255)

                                                                                                  ggCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTaa                                            <  1:304824/150‑1 (MQ=255)

                                                                                                    cGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAAcc                                          <  2:1690209/150‑1 (MQ=255)

                                                                                                     gAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCTGCTCAACCTTGAGCGGCttg                                              <  1:328348/145‑3 (MQ=255)

                                                                                                     gAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCTGCTCAACCTTGAGCGGCttg                                              >  2:328348/1‑143 (MQ=255)

                                                                                                       gTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCa                                       <  2:1382388/150‑1 (MQ=255)

                                                                                                           caggtggACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGc                                   <  2:1111986/148‑1 (MQ=255)

                                                                                                           cTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAATGTTGAGCAGCTGTAACCCCACGGc                                   >  2:522583/1‑150 (MQ=255)

                                                                                                            tggtggACGGGTTCGAGCAGTTATTTACGCAGTCCGGCTGGACTACTTCTGTTGAGTGCAACGGCGGCTAGAAGGGCAAGGCCAAGGCCTCGTGCAAGAAAACGTTAGACATCAAGCCGCTTTATGTTGATCAGCTGTATCTCCACTGcc                                  <  2:623818/150‑1 (MQ=255)

                                                                                                            tggtggACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGcc                                  >  1:852747/1‑150 (MQ=255)

                                                                                                            tggtggACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCAAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGcc                                  >  1:1898030/1‑150 (MQ=255)

                                                                                                                  aCGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTacac                            <  1:1355518/150‑1 (MQ=255)

                                                                                                                  aCGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTacac                            >  2:178441/1‑150 (MQ=255)

                                                                                                                      ggTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGcc                        >  2:533226/1‑150 (MQ=255)

                                                                                                                        tCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCAt                      >  2:1256929/1‑150 (MQ=255)

                                                                                                                         cGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAAccc                                         >  2:769187/1‑130 (MQ=255)

                                                                                                                         cGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAAccc                                         <  1:769187/130‑1 (MQ=255)

                                                                                                                          gAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTGTACGGccgtgg                    <  2:1730953/150‑5 (MQ=255)

                                                                                                                             cAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTg                 <  1:345531/150‑1 (MQ=255)

                                                                                                                               ggTATTTACCCATACCGGATGGACTACGACCGTAGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAACAAAAAGAAAGACACCAAGCCGCTAAAGGTTCACAAGCTGTAACCCCACGGCCGTACTCGGCCATTACTTga               <  2:1058338/150‑1 (MQ=255)

                                                                                                                               ggTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGga               <  2:295194/150‑1 (MQ=255)

                                                                                                                               ggTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGga               <  1:784738/150‑1 (MQ=255)

                                                                                                                               ggTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGga               <  1:766919/150‑1 (MQ=255)

                                                                                                                                gTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAgg                                                          <  1:475245/106‑1 (MQ=255)

                                                                                                                                gTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAgg                                                          >  2:475245/1‑106 (MQ=255)

                                                                                                                                      tcgtaGTCCAGCCGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCa        >  2:859793/5‑150 (MQ=255)

                                                                                                                                      aCCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGtt                                                        >  1:219933/1‑102 (MQ=255)

                                                                                                                                      aCCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGtt                                                        <  2:219933/102‑1 (MQ=255)

                                                                                                                                      aCCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGACATCACTGGAGAAACCa        <  1:1499069/150‑1 (MQ=255)

                                                                                                                                       cccAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTACCCACACGGACGTACACGGCCATCACTGGAGAAACCaa       <  2:117810/150‑1 (MQ=255)

                                                                                                                                       cccAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTg                 <  1:1136276/140‑1 (MQ=255)

                                                                                                                                       cccAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTg                 >  2:1136276/1‑140 (MQ=255)

                                                                                                                                       cccAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCaa       <  1:1806057/150‑1 (MQ=255)

                                                                                                                                       cccAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCaa       <  2:1518446/150‑1 (MQ=255)

                                                                                                                                       cccAGGCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCaa       <  1:400778/150‑1 (MQ=255)

                                                                                                                                        gtaGTCCAGCCGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAAt      >  2:267637/3‑150 (MQ=255)

                                                                                                                                        ccAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAAt      <  2:620354/150‑1 (MQ=255)

ATCGCAGCACCGCCGGCGTGCGGCTGGAAATGCCCGGCCGCACCGACCTGTTTGCCGAACGCTCGATCAACGCCCAAGGCTTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCT  >  NC_002947/3446446‑3446735

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