BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,923,762	A→C	6.1%	L242R (CTG→CGG) ‡	gcvP‑II ←	glycine dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,923,762	0	A	C	6.1%	378.0 / 14.5	131	L242R (CTG→CGG) ‡	gcvP‑II	glycine dehydrogenase
Reads supporting (aligned to +/- strand): ref base A (65/58); new base C (4/4); total (69/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.21e-01

Only 100 of 131 total aligned reads displayed.

GTCGCGACAGGCGAAATAGGCCGCGTGAGGCCCGCCGTAGCCCATCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAGGCGCCGAACACCGAGTGCTTGGCAAGGTTATCCACAGAATCGACGATCAGCTCAAAGCCGAA  >  NC_002947/5923614‑5923904

gtcgCGACAGGCGAAATAGGCCGCGTGAGGCCCGCCGTAGCCCATCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAg                                                                                                                                                  >  2:1041454/1‑150 (MQ=255)

 tcgCGACAGGCGAAATAGGCCGCGTGAGGCCCGCCGTAGCCCATCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAgg                                                                                                                                                 <  1:272998/150‑1 (MQ=255)

 tcgCGACAGGCGAAATAGGCCGCGTGAGGCCCGCCGTAGCCCATCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAgg                                                                                                                                                 <  1:1238619/150‑1 (MQ=255)

 tcgCGACAGGCGAAATAGGCCGCGTGAGGCCCGCCGTAGCCCATCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAgg                                                                                                                                                 >  1:1476178/1‑150 (MQ=255)

 tcgCGACAGGCGAAATAGGCCGCGTGAGGCCCGCCGTAGCCCATCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAgg                                                                                                                                                 <  2:1091251/150‑1 (MQ=255)

   gcgACAGGCGAAATAGGCCGCGTGAGGCCCGCCGTAGCCCATCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGc                                                                                                                                               >  2:445238/1‑150 (MQ=255)

      aCAGGCGAAATAGGCCGCGTGAGGCCCGCCGTAGCCCATCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCctg                                                                                                                                            <  1:1426785/150‑1 (MQ=255)

        aGGCGAAATAGGCCGCGTGAGGCCCGCCGTAGCCCATCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCctgct                                                                                                                                          >  2:155215/1‑150 (MQ=255)

        aGGCGAAATAGGCCGCGTGAGGCCCGCCGTAGCCCATCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCctgct                                                                                                                                          <  2:2112846/150‑1 (MQ=255)

         ggCGAAATAGGCCGCGTGAGGCCCGCCGTAGCCCATCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCGGCAGGTCAGCGGCTACGCAGGCCAGGACctgctg                                                                                                                                         >  2:764452/1‑150 (MQ=255)

         ggCGAAATAGGCCGCGTGAGGCCCGCCGTAGCCCATCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCctgctg                                                                                                                                         >  1:275720/1‑150 (MQ=255)

          acGAAATAGGCCGCGTGAGGCCAGCCGTAGCCCATCGGCACACCAAAGCGCTGGGTCGAGCTGAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCCCAGCAGGTCGGCGGCTACGCAGGCCAGGGCCTgctgg                                                                                                                                        >  2:899122/2‑150 (MQ=255)

             aaaTAGGCCGCGTGAGGCCCGCCGTAGCCCATCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTgctggctg                                                                                                                                     >  1:2025064/1‑150 (MQ=255)

                tAGGCCGCGTGAGGCCCGCCGTAGCCCATCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑c                                                                                                                                 <  2:664402/150‑1 (MQ=255)

                   gCCGCGTGAGGCCCGCCGTAGCCCATCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAAc                                                                                                                              >  1:812421/1‑150 (MQ=255)

                    ccGCGTGAGGCCCGCCGTAGCCCATCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACt                                                                                                                             >  1:361524/1‑150 (MQ=255)

                     cgcgTGAGGCCCGCCGTAGCCCATCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTg                                                                                                                            <  2:438377/150‑1 (MQ=255)

                             gCCCGCCGTAGCCCATCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCa                                                                                                                  >  2:795333/1‑150 (MQ=255)

                             gCCCGCCGTAGCCCATCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCa                                                                                                                  >  2:13475/1‑150 (MQ=255)

                               ccgccgTAGCCCATCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGc                                                                                                                >  2:1333823/1‑150 (MQ=255)

                               ccgccgTAGCCCATCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGc                                                                                                                >  2:904444/1‑150 (MQ=255)

                               ccgccgTAGCCCATCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGc                                                                                                                <  1:1334374/150‑1 (MQ=255)

                                          cATCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑c                                                                                                                                 <  1:1388492/124‑1 (MQ=255)

                                          cATCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑c                                                                                                                                 >  2:1388492/1‑124 (MQ=255)

                                          cATCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATc                                                                                                     <  1:1560309/150‑1 (MQ=255)

                                          cATCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATc                                                                                                     >  2:2180125/1‑150 (MQ=255)

                                           aTCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAAc                                                                                                                              >  2:2083933/1‑126 (MQ=255)

                                           aTCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAAc                                                                                                                              <  1:2083933/126‑1 (MQ=255)

                                           aTCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATcg                                                                                                    >  2:1019895/1‑150 (MQ=255)

                                             cGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATcgcg                                                                                                  <  1:1675321/150‑1 (MQ=255)

                                               gcccacCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCa                                                                                                >  2:601790/4‑150 (MQ=255)

                                               gCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGt                                                                                                                          >  2:1548348/1‑126 (MQ=255)

                                                cacaCCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCAc                                                                                               >  1:1262926/1‑150 (MQ=255)

                                                   gccAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCctgc                                                                                                                                           <  1:714218/105‑1 (MQ=255)

                                                   gccAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCctgc                                                                                                                                           >  2:714218/2‑106 (MQ=255)

                                                     cAACGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTgctgg                                                                                                                                        <  1:618464/107‑1 (MQ=255)

                                                     cAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTgctgg                                                                                                                                        >  2:618464/1‑107 (MQ=255)

                                                     cAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGc                                                                                          >  1:1056156/1‑150 (MQ=255)

                                                      aaaGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGcc                                                                                         >  1:1040852/1‑150 (MQ=255)

                                                      aaaGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGcc                                                                                         <  1:137031/150‑1 (MQ=255)

                                                        aGCGCTGGGTCGAGCCCAGCACCACCTTGCCGCCGAGCTCTCCGGGCGGGCCCAGCACCACCAGGCTCCCCAGGTTGGGGGCTACGCAGGGCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGt                                                                                       <  1:899122/150‑1 (MQ=255)

                                                          cgcTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTgg                                                                                     <  1:155215/150‑1 (MQ=255)

                                                            cTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGt                                                                                       >  2:2351955/1‑146 (MQ=255)

                                                            cTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGt                                                                                       <  1:2351955/146‑1 (MQ=255)

                                                               ggTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTc                                                                                >  2:2003683/1‑150 (MQ=255)

                                                                gTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTcc                                                                               <  2:361524/150‑1 (MQ=255)

                                                                   gAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCggg                                                                            <  1:445238/150‑1 (MQ=255)

                                                                           caccacATCGGCGCCGAGCTCGCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGCGCGCAAAGGGCGCACTTCGCCGTGGATGTCCGGCTACTgcag                                                                    <  2:357033/150‑1 (MQ=255)

                                                                             ccacATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGcc                                                                                         <  2:1325163/127‑1 (MQ=255)

                                                                             ccacATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGcc                                                                                         >  1:1325163/1‑127 (MQ=255)

                                                                              cacaTCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTgc                                                                      >  2:139581/1‑145 (MQ=255)

                                                                              cacaTCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCAG‑TAAAAGAGATAGCTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTACGGGTACTGCAGCaa                                                                 >  1:1353890/1‑150 (MQ=255)

                                                                              atcaTCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCaa                                                                 <  1:1512736/148‑1 (MQ=255)

                                                                                caTCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAgg                                                               <  2:659024/150‑1 (MQ=255)

                                                                                 aTCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTgg                                                                                     >  2:1305839/1‑127 (MQ=255)

                                                                                 aTCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAGgc                                                              <  1:1037677/150‑1 (MQ=255)

                                                                                  tCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTgcag                                                                    >  1:1303998/1‑143 (MQ=255)

                                                                                  tCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAGgcg                                                             >  2:614005/1‑150 (MQ=255)

                                                                                       gcCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGc                                                                                                                                               >  1:730783/1‑66 (MQ=255)

                                                                                          gAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTgg                                                                                     <  2:423065/118‑1 (MQ=255)

                                                                                          gAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGACG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAGGCGCCGAacac                                                     >  1:1850860/1‑150 (MQ=255)

                                                                                             ctctCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGTGCCTAAATCGCGCACTTCGCCGTGGGTGTCCGGTTTCTGTAATGATGCGGCGAACACCGa                                                  <  2:962447/150‑1 (MQ=255)

                                                                                                tcCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCAACGGGTGTCCGGGTACTGCAGCAAGGCGCCGAACACCGAGTg                                               >  1:1600568/1‑150 (MQ=255)

                                                                                                 ccGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGATGTG‑CAACTGGTCG‑‑ATCAGCGGGCTCAAATCGCGCACTTCGCCGTGGGTTTCCGGGTACTGCAGCAAGGCGCCGAACACCGAGTGc                                              <  2:275720/150‑1 (MQ=255)

                                                                                                  cGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTgg                                                                                     >  2:1619902/1‑110 (MQ=255)

                                                                                                  cGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTgg                                                                                     <  1:1619902/110‑1 (MQ=255)

                                                                                                   gggCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAGGCGCCGAACACCGAGTGCtt                                            <  1:1333823/150‑1 (MQ=255)

                                                                                                       ggTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAGGCGCCGAACACCGAGTGCTTGGCa                                        >  1:334593/1‑150 (MQ=255)

                                                                                                        gTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAGGCGCCGAACACCGAGTGCTTGGCaa                                       <  1:1959545/150‑1 (MQ=255)

                                                                                                           ccAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAGGCGCCGAACACCGAGTGCTTGGCAAGGt                                    >  2:2115104/1‑150 (MQ=255)

                                                                                                           ccAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGAGTCCGGGTGCTGCAGCTATTCGCCGAACACCGATGGCTTGGCAAGGt                                    >  2:1932406/1‑150 (MQ=255)

                                                                                                               caccacCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAGGCGCCGAACACCGAGTGCTTGGCAAGGTTATc                                <  1:1041454/150‑1 (MQ=255)

                                                                                                                 ccaccaGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAGGCGCCGAACACCGAGTGCTTGGCAAGGTTATCca                              <  2:1040852/150‑1 (MQ=255)

                                                                                                                  caccaGGCTCAGCAGGTCGGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAGGCGCCGAACACCGAGTGCTTGGCAAGGTTATCcaa                             >  2:249771/1‑149 (MQ=255)

                                                                                                                  caccaGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑At                                                                                                                    >  2:1448712/1‑63 (MQ=255)

                                                                                                                  caccaGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑At                                                                                                                    <  1:1448712/63‑1 (MQ=255)

                                                                                                                  caccaGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAGGCGCCGAACACCGAGTGCTTGGCAAGGTTATCcac                             <  1:795333/150‑1 (MQ=255)

                                                                                                                   accaGGCTCAGCAGGTCAGCGGCTAAGGAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAGGCGCCGAACACCGAGTGCTTGGCAAGGTTATCcaca                            >  1:64569/1‑150 (MQ=255)

                                                                                                                    ccaGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAGGCGCCGAACACCGAGTGCTTGGCAAGGTTATCCACAg                           >  1:403576/1‑150 (MQ=255)

                                                                                                                     caGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAGGCGCCGAACACCGAGTGCTTGGCAAGGTTATCCACAGa                          >  1:1784593/1‑150 (MQ=255)

                                                                                                                      aGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAGGCGCCGAACACCGAGTGCTTGGCAAGGTTATCCACAGaa                         >  2:1790476/1‑150 (MQ=255)

                                                                                                                          tCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAGGCGCCGAACACCGAGTGCTTGGCAAGGTTATCCACAGAATcga                     >  1:920899/1‑150 (MQ=255)

                                                                                                                           cagcagGTCAGCGGCTACGCAGACAAGGGACTACTGGCTGTC‑CATATGCGCG‑‑ATTATCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAGGCGCCGAACACCGAGTGCTTGGCAAGGTTATCCACAGAATcgac                    <  2:1600568/150‑1 (MQ=255)

                                                                                                                             gcagGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAGGCGCCGAACACCGAGTGCTTGGCAAGGTTATCCACAGAATcgacga                  >  2:2062444/1‑150 (MQ=255)

                                                                                                                              cagGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCggg                                                                            >  1:589279/1‑91 (MQ=255)

                                                                                                                              cagGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAGGCGCCGAACACCGAGTGCTTGGCAAGGTTATCCACAGAATCGACGAt                 <  2:403576/150‑1 (MQ=255)

                                                                                                                                  atctcGGGTACGCAGGCGAGGTCCTGCTGGCTGTG‑CAACTGGTCGCTACCCACGCCAGCACCAGGCGCAGCTCG‑‑GTTGGGGTCCGGGGACAACAGCAAGGAGCAGCACAACGAGGCCTTGGAAAGGTGATCCACAGAATCGACGATCAGc             <  1:551884/146‑1 (MQ=255)

                                                                                                                                   cAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGGGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAGGCGCNGAANACNGAGTGCTTGGCANGGTTATCCACAGAATCGACGATCAGCt            <  1:2305133/150‑1 (MQ=255)

                                                                                                                                     gCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAGGCGCCGAACACCGAGTGCTTGGCAAGGTTATCCACAGAATCGACGATCAGCTCa          >  1:2384141/1‑150 (MQ=255)

                                                                                                                                       ggCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTgg                                                                                     >  2:2249428/1‑73 (MQ=255)

                                                                                                                                       ggCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTgg                                                                                     <  1:2249428/73‑1 (MQ=255)

                                                                                                                                       ggCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAGGCGCCGAACACCGAGTGCTTGGCAAGGTTATCCACAGAATCGACGATCAGCtagat        >  2:551884/1‑146 (MQ=255)

                                                                                                                                        agccagCAGGCCCTGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAGGCGCCGAACACCGAGTGCtt                                            >  2:1785338/6‑113 (MQ=255)

                                                                                                                                         cTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAGGCGCCGAACACCGAGTGCTTGGCAAGGTTATCCACAGAATCGACGATCAGCTCAAAGc      >  1:2094123/1‑150 (MQ=255)

                                                                                                                                          ccagCAGGCCCTGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTgg                                                                                     <  1:2172068/67‑1 (MQ=255)

                                                                                                                                           cagCAGGCCCTGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAGGCGCCGAACACCGAGTGCTTGGCAAGGTTATcc                               <  2:74419/121‑1 (MQ=255)

                                                                                                                                             gCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGGGTCCGGGTACTGCAGCAAGGCGCCGAACACCGAGTGCTTGGCAAGGTTATCCACAGAATCGACGATCAGCTCAAAGCCGaa  >  2:1981265/1‑150 (MQ=255)

                                                                                                                                                 gCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGgt                                                                                   >  1:964331/1‑65 (MQ=255)

                                                                                                                                                   cAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAGGCGCCGAACACCGAGTGCTTGGCAAGGTTATCCACAGAATCGACGAt                 >  2:249259/1‑129 (MQ=255)

                                                                                                                                                   cAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAGGCGCCGAACACCGAGTGCTTGGCAAGGTTATCCACAGAATCGACGAt                 <  1:249259/129‑1 (MQ=255)

GTCGCGACAGGCGAAATAGGCCGCGTGAGGCCCGCCGTAGCCCATCGGCACACCAAAGCGCTGGGTCGAGCCCAGCACCACATCGGCGCCGAGCTCTCCGGGCGGTGCCAGCACCACCAGGCTCAGCAGGTCAGCGGCTACGCAGGCCAGGGCCTGCTGGCTGTG‑CAACTGGTCG‑‑ATCAGCGGGCGCAAATCGCGCACTTCGCCGTGGGTGTCCGGGTACTGCAGCAAGGCGCCGAACACCGAGTGCTTGGCAAGGTTATCCACAGAATCGACGATCAGCTCAAAGCCGAA  >  NC_002947/5923614‑5923904

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