BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,390,147	G→A	61.1%	A151V (GCC→GTC)	mpaA ←	murein peptide amidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,390,147	0	G	A	61.1%	77.2 / 132.5	125	A151V (GCC→GTC)	mpaA	murein peptide amidase A
Reads supporting (aligned to +/- strand): ref base G (21/27); new base A (33/43); total (54/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.29e-01

Only 100 of 131 total aligned reads displayed.

AGAGGCAATTCAAACGCCTGAGCCAGCCATTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATACGACATCGCGTTCTTCAGCG  >  NC_000913/1390051‑1390245

agagGCAATTCAAACGCCTGAGCCAGCCATTCGCCTAATTCTCTGTGTCTTGGATCTTCAATACATGTCATTTGATCGTGGAATGTGTCCACTtaggcag                                                                                                 >  2:3646741/1‑100 (MQ=255)

agagGCAATTCAAACGCCTGAGCCAGCCATTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCcaggcag                                                                                                 >  1:3490128/1‑100 (MQ=255)

agagGCAATTCAAACGCCTGAGCCAGCCATTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCcagacag                                                                                                 >  2:463921/1‑100 (MQ=255)

agagGCAATTCAAACGCCTGAGCCAGCCATGCGCCTAATTCGCTGTGTCTGGGATCTTCAAGACAGGCAAGTGGATCGTGGAATGGGACCACCcaggca                                                                                                  <  2:654213/99‑1 (MQ=255)

  agGCAATTCAAACGCCTGAGCCAGCCATTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGt                                                                                               <  1:3184731/100‑1 (MQ=255)

    gCAATTCAAACGCCTGAGCCAGCCATTCGCCTAATTCGCTGTGTCTGGGATCTGCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTg                                                                                             <  2:2305859/100‑1 (MQ=255)

    cgatctCAAACGCCTGAGCCAGCCATTCGCCTAATTCGCTGTGTCTGTGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTg                                                                                             <  1:763702/95‑1 (MQ=255)

     cAATTCAAACGCCTGAGCCAGCCATTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGa                                                                                            <  2:3913159/100‑1 (MQ=255)

       aTTCAAACGCCTGAGCCAGCCATTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAAt                                                                                          <  1:1394730/100‑1 (MQ=255)

           aaaCGCCTGAGCCAGCCATTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCgg                                                                                      >  1:521740/1‑100 (MQ=255)

           aaaCGCCTGAGCCAGCCATTCGCCTAATTCGCTATATCTGGGATCTTCAATACAGGCCAGTGAATCGTGGAATGAGACAACCCAGGCAAGTTGAATGCgg                                                                                      >  2:1687686/1‑100 (MQ=255)

            aaCGCCTGAGCCAGCCATTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGgt                                                                                     <  1:54377/100‑1 (MQ=255)

            aaCGCCTGAGCCAGCCATTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGgt                                                                                     <  1:3771896/100‑1 (MQ=255)

              cGCCTGAGCCAGCCATTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGgtgt                                                                                   >  1:1388977/1‑100 (MQ=255)

                ccTGAGCCAGCCATTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGtat                                                                                 >  2:2882437/1‑100 (MQ=255)

                ccTGAGCCAGCCATTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGtat                                                                                 >  1:145524/1‑100 (MQ=255)

                ccTGAGCAAGACATTCGCCTAATAAGAAGTGTCCGGGAGCTCCAATAGAGGTCAGTGGGTTGTCGTACAAGACCACGGCGACAGGTTGAATGCGgtgt                                                                                   <  2:2406983/98‑1 (MQ=255)

                 cTGAGCCAGCCATTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGtata                                                                                <  1:1358091/100‑1 (MQ=255)

                     gccagccaTTCGCCTAATTCGCTGTGTCTGGGATCTCCAATACAGGCCAGTGGGTCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCt                                                                            <  2:2670599/100‑1 (MQ=255)

                     gccagccaTTCGCCTAATGCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCt                                                                            <  1:2548635/100‑1 (MQ=255)

                       cagccaTTCGCCAAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTgg                                                                          >  1:2715304/1‑100 (MQ=255)

                         gccaTTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGATTGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTggta                                                                        <  1:2794604/100‑3 (MQ=255)

                         gccaTTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGca                                                                        <  1:80781/100‑1 (MQ=255)

                          ccaTTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGcac                                                                       >  1:1000916/1‑100 (MQ=255)

                             ttCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAgc                                                                    <  1:2517057/100‑1 (MQ=255)

                               cGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAgcgc                                                                  <  1:595282/100‑1 (MQ=255)

                               cGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTCGAATGCGGTGTATAAGCAGGCACAgcg                                                                   >  2:794490/1‑99 (MQ=255)

                                nccTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGcc                                                                 >  2:3677043/2‑100 (MQ=255)

                                gCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGcc                                                                 >  2:3897081/1‑100 (MQ=255)

                                gCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGcc                                                                 >  1:396098/1‑100 (MQ=255)

                                 ccTAATTCGCTGTGTCTGGGATCTTCAATGCAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCt                                                                <  1:2560014/100‑1 (MQ=255)

                                 ccTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCt                                                                >  1:2756314/1‑100 (MQ=255)

                                  cTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTg                                                               <  1:1454767/100‑1 (MQ=255)

                                   ttaTTGGCTGTGTGTGGGTTCTTCTATGCAGGTGAGGGGATCGTGGAATGAGACTTCTCATTCTGGTTGATTGCGGTGTATGAGCTGGGACAgcgcg                                                                 <  2:128518/95‑2 (MQ=255)

                                    aaTTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTTAATGCGGTGTATAAGCTGGCACAGCGCCTggg                                                             <  1:3922355/100‑1 (MQ=255)

                                    aaTTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTggg                                                             <  1:3665418/100‑1 (MQ=255)

                                    aaTTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGATGAATGCGGTGTATAAGCTGGCACAGCGCCTggg                                                             >  2:2179908/1‑100 (MQ=255)

                                     aTTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGTGCCTGGGt                                                            >  1:657916/1‑100 (MQ=255)

                                     aTTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGt                                                            <  1:766136/100‑1 (MQ=255)

                                      ttCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGtt                                                           >  2:1669552/1‑100 (MQ=255)

                                         gCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCa                                                        <  1:3501313/100‑1 (MQ=255)

                                         gCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTACAAGCTGGCACCGCGCCTGGGTTTCa                                                        >  2:1882435/1‑100 (MQ=255)

                                           tgtgTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAgg                                                      <  1:3001887/100‑1 (MQ=255)

                                            gtgtCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGt                                                     <  1:2064665/100‑1 (MQ=255)

                                              gtCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTc                                                   >  2:646293/1‑100 (MQ=255)

                                               tCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCg                                                  <  1:2394268/100‑1 (MQ=255)

                                                cTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCgg                                                 >  1:1972463/1‑100 (MQ=255)

                                                cTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACATGTTGAATGCGGTGTATAAACTGACACAGCGCCtccgcttca                                                        >  2:400862/1‑85 (MQ=255)

                                                cTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCgg                                                 >  1:1468837/1‑100 (MQ=255)

                                                cTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCg                                                  >  2:1164813/1‑99 (MQ=255)

                                                cTGGCATCTGCAAAACAGGCCAGTCGATCGTGCAATGAGCCCACCCAGGCATGTTGAATGCGGTGTATAAGCTTGCACAGCGCCCCGGTTTCAGGTTCg                                                  <  2:3645494/99‑1 (MQ=255)

                                                  gggATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACATCGCCTGGGTTTCAGGTTCGGa                                                <  1:3982112/99‑1 (MQ=255)

                                                  gggATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAg                                               <  1:1673260/100‑1 (MQ=255)

                                                  gggATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAg                                               <  1:1293586/100‑1 (MQ=255)

                                                  gggATCTTCAATACAGGCCAGTGGAGCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAg                                               <  1:3149837/100‑1 (MQ=255)

                                                   ggATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGc                                              <  2:1398319/100‑1 (MQ=255)

                                                    gATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGcc                                             >  1:2291011/1‑100 (MQ=255)

                                                     aTCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCa                                            <  1:3714215/100‑1 (MQ=255)

                                                       cTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAgg                                          <  1:2839293/100‑1 (MQ=255)

                                                       cTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAgg                                          >  1:2749363/1‑100 (MQ=255)

                                                           aaTACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTg                                      >  1:294239/1‑100 (MQ=255)

                                                           aaTACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTg                                      <  1:510120/100‑1 (MQ=255)

                                                            aTACAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGt                                     >  2:2233728/1‑100 (MQ=255)

                                                               cAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGCTTGTCAc                                  >  2:525111/1‑100 (MQ=255)

                                                               cAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCAc                                  <  1:2582532/100‑1 (MQ=255)

                                                               cAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCAc                                  >  1:1116957/1‑100 (MQ=255)

                                                               cAGGCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCAc                                  <  1:2580792/100‑1 (MQ=255)

                                                                aGGGCAGTGGAGCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCAcc                                 <  2:3813933/100‑1 (MQ=255)

                                                                 ggCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGGTTCAGGTTCGGAGCCAGGTTTGTCACCt                                >  2:3148959/1‑100 (MQ=255)

                                                                 ggCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCt                                >  1:928205/1‑100 (MQ=255)

                                                                 ggCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCt                                >  1:1908399/1‑100 (MQ=255)

                                                                  gCCAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTg                               <  2:3589692/100‑1 (MQ=255)

                                                                    cAGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTc                             <  2:1468826/100‑1 (MQ=255)

                                                                     aGTGGATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCa                            <  2:8316/100‑1 (MQ=255)

                                                                        ggATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCCGTa                         >  2:1968806/1‑100 (MQ=255)

                                                                        ggATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTa                         <  2:4002574/100‑1 (MQ=255)

                                                                        ggATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCGGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTa                         >  2:2437772/1‑100 (MQ=255)

                                                                         gATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTaa                        <  1:1926683/100‑1 (MQ=255)

                                                                         gATCGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGAGTCAGGCTCGGAGCCAGGTTTGTCACCCGTCAGt                          >  2:134511/1‑98 (MQ=255)

                                                                          aTCGTGGAATGAGACCCCCCAGGCACGTTGAATGCGGTGTATAAGCTGGCACCGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCCGTCAt                       >  2:1384811/1‑100 (MQ=255)

                                                                           tCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATa                      <  2:2377276/100‑1 (MQ=255)

                                                                            cGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATAc                     >  1:250516/1‑100 (MQ=255)

                                                                            cGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATAc                     <  2:2953682/100‑1 (MQ=255)

                                                                            cGTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATa                      <  2:304246/99‑1 (MQ=255)

                                                                             gTGGAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATACg                    >  1:3265123/1‑100 (MQ=255)

                                                                               ggAATGAGACCACCCAGACAGGTTGACTGCGGTGTATAAGCTGGCACAGCGCTTGGGTTTGAGGTTCGGAGCCAGGTTTGTCACCTGTCCGTAATACGAc                  >  2:1540721/1‑100 (MQ=255)

                                                                                gAATGAGACCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATACGACa                 >  1:1001223/1‑100 (MQ=255)

                                                                                     agaCCGCCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATACGACATcgcg            <  2:1864252/100‑1 (MQ=255)

                                                                                      gaCCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATACGACATCGCGt           <  2:1722559/100‑1 (MQ=255)

                                                                                      gaCCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATACGACATCGCGt           <  1:2179921/100‑1 (MQ=255)

                                                                                      gaCCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATACGACATCGCGt           <  2:222657/100‑1 (MQ=255)

                                                                                      gaCCACCCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATACGACATCGCGt           >  1:1707171/1‑100 (MQ=255)

                                                                                       accaccCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATACGACATCGCGtt          >  1:2077976/1‑100 (MQ=255)

                                                                                       accaccCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATACGACATCGCGtt          <  1:3938823/100‑1 (MQ=255)

                                                                                       accaccCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATACGACATCGCGtt          <  1:3605483/100‑1 (MQ=255)

                                                                                        ccaccCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATACGACATCGCGttc         <  2:2853924/100‑1 (MQ=255)

                                                                                         caccCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATACGACATCGCGtt          <  2:1388968/98‑1 (MQ=255)

                                                                                         caccCAGACAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATACGACATCGCGttc         >  2:886740/1‑99 (MQ=255)

                                                                                             caggcaggTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATACGACATCGCGTTCTTCAg    >  1:936156/1‑100 (MQ=255)

                                                                                               ggcaggTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATACGACATCGCGTTCTTCAgcg  <  1:646299/100‑1 (MQ=255)

AGAGGCAATTCAAACGCCTGAGCCAGCCATTCGCCTAATTCGCTGTGTCTGGGATCTTCAATACAGGCCAGTGGATCGTGGAATGAGACCACCCAGGCAGGTTGAATGCGGTGTATAAGCTGGCACAGCGCCTGGGTTTCAGGTTCGGAGCCAGGTTTGTCACCTGTCAGTAATACGACATCGCGTTCTTCAGCG  >  NC_000913/1390051‑1390245

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