BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,568,197	A→G	M106T (ATG→ACG)	xasA ←	acid sensitivity protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,568,197	0	A	G	99.1%	692.6 / ‑2.9	219	M106T (ATG→ACG)	xasA	acid sensitivity protein
Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (99/118); total (101/118)
Fisher's exact test for biased strand distribution p-value = 2.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-02

Only 100 of 227 total aligned reads displayed.

ATGAGTGCTGCAATAGTTTTGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCG  >  NC_000913/1568107‑1568286

 tgaGTGCTGCAATAGTTTTGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCg                                                                                        >  2:5551644/1‑93 (MQ=255)

  gaGTGCTGCAATAGTTTTGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCgg                                                                                       <  2:6737037/93‑1 (MQ=255)

     tgctgcAATAGTTTTGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAAt                                                                                    >  2:3610140/1‑93 (MQ=255)

     tgctgcAATAGTTTTGGTAATGGGGTCTTCATTCAGCGCTCGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAAt                                                                                    <  1:5792239/93‑1 (MQ=255)

       ctgcAATAGTTTTGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATaa                                                                                  >  2:5971825/1‑93 (MQ=255)

        tgcAATAGTTTTGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATaaa                                                                                 >  2:7311344/1‑93 (MQ=255)

         gcAATAGTTTTGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATaaaa                                                                                <  1:1880203/93‑1 (MQ=255)

         gcAATAGTTTTGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGCGCCCCTAACACGAAATAGAGCGTCGGAATaaaa                                                                                >  2:9752187/1‑93 (MQ=255)

          cAATAGTTTTGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAAc                                                                               >  1:4114609/1‑93 (MQ=255)

          cAATAGTTTTGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAAc                                                                               >  1:9390637/1‑93 (MQ=255)

           aaTAGTTTTGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAAcc                                                                              <  2:9693851/93‑1 (MQ=255)

           aaTAGTCTTGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAAcc                                                                              <  1:7580675/93‑1 (MQ=255)

             tAGTTTTGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCaa                                                                            <  1:498634/93‑1 (MQ=255)

             aaGTTGTGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCaa                                                                            <  1:3225460/92‑1 (MQ=255)

              aGTTTTGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAAt                                                                           <  2:7540361/93‑1 (MQ=255)

                ttttGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATgg                                                                         >  1:6936418/1‑93 (MQ=255)

                   tGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGa                                                                      <  1:3211682/93‑1 (MQ=255)

                   tGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGa                                                                      <  2:6678517/93‑1 (MQ=255)

                     gTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACCCGAAATAAAGCGTCGGAATAAAACCAATGGCGAtt                                                                    <  1:678482/93‑1 (MQ=255)

                     gTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGAtt                                                                    >  2:7510333/1‑93 (MQ=255)

                     gTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGAtt                                                                    <  2:2045668/93‑1 (MQ=255)

                      tAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGAttt                                                                   >  1:1173942/1‑93 (MQ=255)

                      tAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGAttt                                                                   >  2:4429714/1‑93 (MQ=255)

                       aaTGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTg                                                                  >  1:3299505/1‑93 (MQ=255)

                        aTGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGc                                                                 >  2:3782552/1‑93 (MQ=255)

                         tGGGGTGTTCATGCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCa                                                                <  1:7507211/93‑1 (MQ=255)

                          ggggTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAg                                                               >  1:4819981/1‑93 (MQ=255)

                          ggggTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAg                                                               >  2:6076410/1‑93 (MQ=255)

                          ggggTCTGCATTCAGCGCTGGCCCTTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAg                                                               <  1:1042371/93‑1 (MQ=255)

                            ggTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGAt                                                             <  1:5890551/93‑1 (MQ=255)

                            ggTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGAt                                                             <  2:8804377/93‑1 (MQ=255)

                             gTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATa                                                            <  2:6914015/93‑1 (MQ=255)

                               cTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGc                                                          <  2:5018402/93‑1 (MQ=255)

                                ttcattcaGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGcc                                                         >  2:5049451/1‑93 (MQ=255)

                                 tcattcaGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCa                                                        <  2:909593/93‑1 (MQ=255)

                                 tcattcaGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCa                                                        >  1:573242/1‑93 (MQ=255)

                                   attcaGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCaaa                                                      >  1:2661593/1‑93 (MQ=255)

                                   attcaGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCaaa                                                      <  2:9633750/93‑1 (MQ=255)

                                     tcaGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATg                                                    <  1:5682327/93‑1 (MQ=255)

                                     tcaGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATg                                                    >  1:3585266/1‑93 (MQ=255)

                                       aGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATgag                                                  >  1:4212822/1‑93 (MQ=255)

                                         cgcTGTCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGAAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGAt                                                <  2:9044697/93‑1 (MQ=255)

                                         cgcTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGAt                                                <  2:387685/93‑1 (MQ=255)

                                         cgcTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCAGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGAt                                                <  1:4932410/93‑1 (MQ=255)

                                          gcTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATc                                               >  1:5641730/1‑93 (MQ=255)

                                          gcTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATc                                               >  2:8792776/1‑93 (MQ=255)

                                           cTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCg                                              >  2:4937713/1‑93 (MQ=255)

                                           cTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCg                                              >  2:309638/1‑93 (MQ=255)

                                            tGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGCGCTCGc                                             <  2:8920223/93‑1 (MQ=255)

                                            tGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGc                                             <  1:8093689/93‑1 (MQ=255)

                                            tGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGc                                             <  1:7769534/93‑1 (MQ=255)

                                            tGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGc                                             <  1:7274627/93‑1 (MQ=255)

                                             ggCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCt                                            <  2:7813061/93‑1 (MQ=255)

                                              gCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTg                                           >  2:9128720/1‑93 (MQ=255)

                                               ccATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCANATAGCCAAATGAGATCGCTGc                                          >  1:376232/1‑93 (MQ=255)

                                                 aTTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCaa                                        >  1:1351648/1‑93 (MQ=255)

                                                 aTTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCaa                                        >  1:820133/1‑93 (MQ=255)

                                                  tttCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCaaa                                       >  1:4820978/1‑93 (MQ=255)

                                                  tttCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCaaa                                       >  2:6791381/1‑93 (MQ=255)

                                                    tCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATc                                     <  1:9944478/93‑1 (MQ=255)

                                                       ggATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATcccc                                  >  2:3314352/1‑93 (MQ=255)

                                                        gATGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCa                                 <  2:230145/93‑1 (MQ=255)

                                                         aTGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCAt                                <  1:8885185/93‑1 (MQ=255)

                                                         aTGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGCCGATTTGAAGAAACCCAATTGAGATCGCTGCAAATCCCCac                                >  2:2431109/1‑92 (MQ=255)

                                                         aTGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCAt                                >  1:5971973/1‑93 (MQ=255)

                                                          tGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCAtc                               <  1:7169972/93‑1 (MQ=255)

                                                          tGTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCAtc                               <  1:2749598/93‑1 (MQ=255)

                                                           gTAGGAGAGTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCAtct                              <  1:4370516/93‑1 (MQ=255)

                                                           gTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCAtct                              >  1:5324089/1‑93 (MQ=255)

                                                                   gTGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCa                      >  1:2390737/1‑93 (MQ=255)

                                                                    tGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCag                     >  2:5615816/1‑93 (MQ=255)

                                                                    tGCCCCTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCag                     <  2:5733971/93‑1 (MQ=255)

                                                                      ccccTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCagag                   >  2:8564448/1‑93 (MQ=255)

                                                                      ccccTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCagag                   >  2:7634454/1‑93 (MQ=255)

                                                                      ccccTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCagag                   <  2:7260874/93‑1 (MQ=255)

                                                                        ccTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTa                 >  2:4966104/1‑93 (MQ=255)

                                                                         cTAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTAt                >  2:465071/1‑93 (MQ=255)

                                                                          tAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTAtt               >  2:2134604/1‑93 (MQ=255)

                                                                          tAACACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTAtt               >  2:849154/1‑93 (MQ=255)

                                                                           aaCACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTAttt              >  1:219017/1‑93 (MQ=255)

                                                                           aaCACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTAttt              <  1:8592434/93‑1 (MQ=255)

                                                                           aaCACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTAttt              <  2:8092928/93‑1 (MQ=255)

                                                                           aaCACGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCAAAATGAGATCGCGGCAAATCCCCATCTCGGCCCCAGAGTAttt              <  2:6351045/93‑1 (MQ=255)

                                                                            acacGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTg             <  2:5862701/93‑1 (MQ=255)

                                                                             cacGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGa            <  1:7965337/93‑1 (MQ=255)

                                                                             cacGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGa            <  1:9119667/93‑1 (MQ=255)

                                                                             cacGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGa            <  2:2234886/93‑1 (MQ=255)

                                                                             cacGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGa            <  2:279710/93‑1 (MQ=255)

                                                                              acGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGAt           <  1:8246693/93‑1 (MQ=255)

                                                                              acGAAATAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGAt           >  1:3050083/1‑93 (MQ=255)

                                                                                 aaaaAGAGCGTAGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCCGCAAATCCCCATCTCGGCCCCAGAGTATTTGATAcc        <  2:4203577/93‑1 (MQ=255)

                                                                                 aaaTAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATAcc        >  2:3493979/1‑93 (MQ=255)

                                                                                 aaaTAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATAcc        <  2:3833639/93‑1 (MQ=255)

                                                                                  aaTAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATAccc       >  2:1557346/1‑93 (MQ=255)

                                                                                   aTAGAGCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCa      >  1:6963558/1‑93 (MQ=255)

                                                                                     agagCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTTCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAgg    <  2:4243425/93‑1 (MQ=255)

                                                                                     agagCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAgg    <  1:1433391/93‑1 (MQ=255)

                                                                                     agagCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAgg    <  1:8438678/93‑1 (MQ=255)

                                                                                     agagCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAgg    <  1:3846619/93‑1 (MQ=255)

                                                                                       agCGTCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGca  >  1:8311822/1‑92 (MQ=255)

ATGAGTGCTGCAATAGTTTTGGTAATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCG  >  NC_000913/1568107‑1568286

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