BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,861,203	A→G	A282A (GCA→GCG)	b2737 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,861,203	0	A	G	99.4%	499.7 / ‑5.1	160	A282A (GCA→GCG)	b2737	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (90/69); total (91/69)
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.95e-01

Only 100 of 169 total aligned reads displayed.

CGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCACTG  >  NC_000913/2861113‑2861290

cgTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGcgcg                                                                                       >  1:27053/1‑93 (MQ=255)

 gTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGgcgcgc                                                                                      <  1:2127282/93‑1 (MQ=255)

 gTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGgcgcgc                                                                                      <  2:4255022/93‑1 (MQ=255)

 gTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGgcgcgc                                                                                      >  1:3870254/1‑93 (MQ=255)

 gTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGgcgcgc                                                                                      >  1:1549556/1‑93 (MQ=255)

 gTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGgcgcgc                                                                                      <  2:1656990/93‑1 (MQ=255)

  tAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCt                                                                                     <  2:4664082/93‑1 (MQ=255)

    gTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGc                                                                                   <  1:4389583/93‑1 (MQ=255)

    gTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGc                                                                                   >  1:1279651/1‑93 (MQ=255)

    gTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGc                                                                                   >  1:1828642/1‑93 (MQ=255)

      gCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCct                                                                                 <  1:910481/93‑1 (MQ=255)

       ctctCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCctc                                                                                <  1:486481/93‑1 (MQ=255)

        tctcCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCctct                                                                               >  1:1374999/1‑93 (MQ=255)

         ctcCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCtctc                                                                              >  1:3873947/1‑93 (MQ=255)

          tcCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCa                                                                             >  1:3200417/1‑93 (MQ=255)

          tcCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCa                                                                             >  1:2745557/1‑93 (MQ=255)

             ggTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAAtt                                                                          <  2:4067822/93‑1 (MQ=255)

             ggTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAAtt                                                                          <  1:3162681/93‑1 (MQ=255)

               tttaTGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGa                                                                        >  1:3236218/4‑93 (MQ=255)

                  aTGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAAc                                                                     >  2:214333/1‑93 (MQ=255)

                  aTGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAAc                                                                     <  2:4183865/93‑1 (MQ=255)

                  aTGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAAc                                                                     <  1:1751988/93‑1 (MQ=255)

                     ctctCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAActct                                                                  <  1:1898759/93‑1 (MQ=255)

                      tctcAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTg                                                                 >  1:2921384/1‑93 (MQ=255)

                      tctcAAATGACCAACCGCCAAGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTg                                                                 >  2:1783251/1‑93 (MQ=255)

                        tcAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGc                                                               <  2:3631113/93‑1 (MQ=255)

                           aaTGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGc                                                            >  2:1468553/1‑93 (MQ=255)

                           aaTGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGc                                                            >  1:2075921/1‑93 (MQ=255)

                            aTGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCt                                                           >  2:4437319/1‑93 (MQ=255)

                            aTGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCt                                                           <  2:4041432/93‑1 (MQ=255)

                             tGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCtt                                                          >  1:2431305/1‑93 (MQ=255)

                             tGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCtt                                                          >  1:2242404/1‑93 (MQ=255)

                              gACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTa                                                         <  1:2591645/93‑1 (MQ=255)

                              gACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTa                                                         >  1:2464639/1‑93 (MQ=255)

                              gACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTa                                                         <  1:2283425/93‑1 (MQ=255)

                               aCCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTAt                                                        <  1:2645840/93‑1 (MQ=255)

                               aCCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGAGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTAt                                                        <  1:2573277/93‑1 (MQ=255)

                                ccAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATg                                                       <  1:1104064/93‑1 (MQ=255)

                                 cAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGc                                                      <  1:840151/93‑1 (MQ=255)

                                  aaCCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGca                                                     >  2:2217714/1‑93 (MQ=255)

                                  aaCCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGca                                                     <  1:3932451/93‑1 (MQ=255)

                                     cGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCacac                                                  >  1:2354614/1‑93 (MQ=255)

                                      gCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACg                                                 <  2:276014/93‑1 (MQ=255)

                                      gCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACg                                                 >  1:4431096/1‑93 (MQ=255)

                                       ccAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGa                                                >  1:4365291/1‑93 (MQ=255)

                                        cAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGac                                               >  1:3759525/1‑92 (MQ=255)

                                         aGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAAc                                              >  2:1455744/1‑93 (MQ=255)

                                         aGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAAc                                              <  2:4122082/93‑1 (MQ=255)

                                         aGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAAc                                              <  1:3639379/93‑1 (MQ=255)

                                          ggTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACt                                             >  1:3701776/1‑93 (MQ=255)

                                          ggTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACt                                             >  1:1963632/1‑93 (MQ=255)

                                          ggTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACt                                             >  2:1949836/1‑93 (MQ=255)

                                              gCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTATGGCCGCTTATGCACACGAACTGGCa                                         >  1:137561/1‑93 (MQ=255)

                                                acaTTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAAATCTGGCCGCTTATGCACACGAACTGGCaga                                       >  1:1572375/1‑93 (MQ=255)

                                                 caTTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCagag                                      >  2:1388562/1‑93 (MQ=255)

                                                  aTTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCagagg                                     >  1:3387292/1‑92 (MQ=255)

                                                  aTTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGt                                     <  1:4315711/93‑1 (MQ=255)

                                                  aTTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGt                                     >  1:2642554/1‑93 (MQ=255)

                                                   ttACCGTCGAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTg                                    >  2:4580655/1‑93 (MQ=255)

                                                   ttACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGTACACGAACTGGCAGAGTg                                    <  2:4202656/93‑1 (MQ=255)

                                                   ttACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTg                                    >  1:961002/1‑93 (MQ=255)

                                                    tACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTgg                                   <  1:103558/93‑1 (MQ=255)

                                                     aCCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTggg                                  <  1:1011258/93‑1 (MQ=255)

                                                     aCCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTggg                                  <  1:1400439/93‑1 (MQ=255)

                                                     aCCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTggg                                  <  1:1311523/93‑1 (MQ=255)

                                                       cGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGtt                                >  2:340729/1‑93 (MQ=255)

                                                        gTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGttt                               >  2:1000375/1‑92 (MQ=255)

                                                        gTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTc                               >  2:2601880/1‑93 (MQ=255)

                                                          cAAATTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTg                             >  1:3980980/1‑93 (MQ=255)

                                                           aaaTTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTgg                            <  2:1073764/93‑1 (MQ=255)

                                                           aaaTTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTgg                            >  2:1800509/1‑93 (MQ=255)

                                                           aaaTTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTgg                            >  1:3137440/1‑93 (MQ=255)

                                                            aaTTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTggg                           >  1:465277/1‑93 (MQ=255)

                                                             aTTGCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGc                          >  2:654492/1‑93 (MQ=255)

                                                                gCACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAg                       <  2:3747047/93‑1 (MQ=255)

                                                                 cACCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAgg                      >  2:3804073/1‑93 (MQ=255)

                                                                  aCCAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGa                     >  2:3037958/1‑93 (MQ=255)

                                                                   ccAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGaa                    >  2:576249/1‑93 (MQ=255)

                                                                   ccAGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGaa                    <  1:1611741/93‑1 (MQ=255)

                                                                     aGCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGAAGAGTGGGTTCTGGGCCAGGAAAg                  >  2:305644/1‑93 (MQ=255)

                                                                      gCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAgt                 <  1:253200/93‑1 (MQ=255)

                                                                      gCCCGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAgt                 <  2:3005899/93‑1 (MQ=255)

                                                                        ccGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAgtgt               >  2:3439323/1‑93 (MQ=255)

                                                                        ccGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAgtgt               >  1:706288/1‑93 (MQ=255)

                                                                        ccGTGAAGTTGATGTGGCGCGCTCCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAgtgt               <  2:1024240/93‑1 (MQ=255)

                                                                         cGTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTa              <  1:1328014/93‑1 (MQ=255)

                                                                          gTGAAGTTGATGTGGCGTTCTGCCTCTCAATTGAAAATCTGGCCGCTTATGCACACGAACTGGCAGAGTGCGTTGTGGGGCAGGAAAGAGTAc             >  2:2726686/1‑93 (MQ=255)

                                                                          gTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGTTTCTGGGCCAGGAAAGTGTAc             >  2:3585239/1‑93 (MQ=255)

                                                                          gTGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTAc             <  1:4228842/93‑1 (MQ=255)

                                                                           tGAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACt            <  2:3196059/93‑1 (MQ=255)

                                                                            gAAGTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACtt           >  2:4198945/1‑93 (MQ=255)

                                                                               gTTGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCt        >  2:3473903/1‑93 (MQ=255)

                                                                                ttGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTc       >  2:828731/1‑93 (MQ=255)

                                                                                ttGATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTc       >  2:879198/1‑93 (MQ=255)

                                                                                  gATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCa     >  2:3798701/1‑93 (MQ=255)

                                                                                  gATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCa     >  2:2781600/1‑93 (MQ=255)

                                                                                  gATGTGGCGCGCTGCCTCTCAATTGAAACTCTGGACGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTccc     >  1:6005/1‑92 (MQ=255)

                                                                                   aTGTGGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCAc    <  1:3276192/93‑1 (MQ=255)

                                                                                    tgtgGCGCGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCACt   >  2:3864713/1‑93 (MQ=255)

                                                                                     gtgGCGCGCTGCCTCTCAATTGAAACTCTGGCCGATTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCACTg  <  1:2176958/93‑1 (MQ=255)

CGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAAAGTGTACTTGCTCCACTG  >  NC_000913/2861113‑2861290

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