BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,127,469	G→A	G254G (GGC→GGT)	b2981 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,127,469	0	G	A	99.4%	617.3 / ‑5.3	181	G254G (GGC→GGT)	b2981	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (88/91); total (89/92)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.15e-01

Only 100 of 183 total aligned reads displayed.

GTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGCATCT  >  NC_000913/3127379‑3127560

gTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCAcc                                                                                           >  2:4675410/1‑93 (MQ=255)

gTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCAcc                                                                                           >  1:360104/1‑93 (MQ=255)

 ttCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCAccc                                                                                          >  1:873996/1‑93 (MQ=255)

  tCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCa                                                                                         <  2:1386352/93‑1 (MQ=255)

  tCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCa                                                                                         >  1:1869919/1‑93 (MQ=255)

    ggCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATg                                                                                       >  1:2541852/1‑93 (MQ=255)

    ggCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATg                                                                                       >  2:3755239/1‑93 (MQ=255)

       aaGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTc                                                                                    <  2:2299777/93‑1 (MQ=255)

       aaGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTc                                                                                    <  2:2794034/93‑1 (MQ=255)

         gCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTcgc                                                                                  <  2:404471/93‑1 (MQ=255)

           cGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGa                                                                                <  2:2078340/93‑1 (MQ=255)

            gACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGaa                                                                               >  1:1771109/1‑93 (MQ=255)

             aCAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAAt                                                                              <  1:1992523/93‑1 (MQ=255)

              cAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAAtt                                                                             <  1:3613075/93‑1 (MQ=255)

               aaTAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTc                                                                            >  1:4707508/1‑93 (MQ=255)

                aTAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCa                                                                           <  2:1596042/93‑1 (MQ=255)

                aTAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCa                                                                           <  1:2522685/93‑1 (MQ=255)

                 tAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCAt                                                                          >  1:80766/1‑93 (MQ=255)

                 tAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCAt                                                                          >  1:4539893/1‑93 (MQ=255)

                 tAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCAt                                                                          >  2:2697891/1‑93 (MQ=255)

                        gcAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGt                                                                   >  1:2976185/1‑93 (MQ=255)

                        gcAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGt                                                                   >  1:1526800/1‑93 (MQ=255)

                          aGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTaa                                                                 <  1:2788562/93‑1 (MQ=255)

                          aGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTaa                                                                 >  2:443426/1‑93 (MQ=255)

                           gggTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTaaa                                                                <  2:1374022/93‑1 (MQ=255)

                           gggTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTaaa                                                                <  1:2956338/93‑1 (MQ=255)

                           gggTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTaaa                                                                >  2:2659405/1‑93 (MQ=255)

                             gTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACg                                                              <  2:2043943/93‑1 (MQ=255)

                             gTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACg                                                              >  2:2201799/1‑93 (MQ=255)

                                 aGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGcac                                                          <  2:3736544/93‑1 (MQ=255)

                                 aGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGcac                                                          >  2:4579194/1‑93 (MQ=255)

                                 aGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGcac                                                          >  2:4239190/1‑93 (MQ=255)

                                  gAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGcaca                                                         >  2:4658371/1‑93 (MQ=255)

                                   aatcaatcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACaa                                                        >  2:2754026/1‑93 (MQ=255)

                                   aatcaatcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTAATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACaa                                                        >  2:4632959/1‑93 (MQ=255)

                                    atcaatcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACaaa                                                       <  2:3710645/93‑1 (MQ=255)

                                    atcaatcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACaaa                                                       >  1:1069244/1‑93 (MQ=255)

                                    atcaatcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACaaa                                                       >  1:4351739/1‑93 (MQ=255)

                                     tcaatcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAAt                                                      >  1:951559/1‑93 (MQ=255)

                                     tcaatcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCACGTTCGGTAAACGGCACAAAt                                                      >  2:3206381/1‑93 (MQ=255)

                                     acaatcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCTCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAAt                                                      <  1:2641312/92‑1 (MQ=255)

                                      caatcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATc                                                     >  2:1317902/1‑93 (MQ=255)

                                       aatcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATcc                                                    <  1:208793/93‑1 (MQ=255)

                                       aatcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATcc                                                    >  1:1589031/1‑93 (MQ=255)

                                       aatcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATcc                                                    >  1:3060795/1‑93 (MQ=255)

                                       aatcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATcc                                                    >  2:1697242/1‑93 (MQ=255)

                                        atcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATccc                                                   >  1:3573906/1‑93 (MQ=255)

                                         tctcAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATcccc                                                  >  1:2804737/1‑93 (MQ=255)

                                          ctcAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCa                                                 <  1:1578168/93‑1 (MQ=255)

                                          ctcAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCa                                                 >  1:645728/1‑93 (MQ=255)

                                           tcAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAg                                                <  2:3408577/93‑1 (MQ=255)

                                           tcAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAg                                                >  1:1806905/1‑93 (MQ=255)

                                            cAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGt                                               >  2:2852406/1‑93 (MQ=255)

                                            cAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGt                                               <  1:2248319/93‑1 (MQ=255)

                                            cAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGt                                               <  2:562497/93‑1 (MQ=255)

                                              gCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTtg                                             <  2:3292579/93‑1 (MQ=255)

                                                  tataGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGtgct                                         >  1:903584/1‑93 (MQ=255)

                                                  tataGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGtgcg                                         >  2:3608111/1‑92 (MQ=255)

                                                   ataGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGtgctg                                        <  1:1106026/93‑1 (MQ=255)

                                                   ataGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGtgctg                                        >  1:40267/1‑93 (MQ=255)

                                                     aGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTCAACGGCACAAATCCCCAGTTGTGCTgcc                                      >  1:2010220/1‑93 (MQ=255)

                                                     aGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgcc                                      >  1:814420/1‑93 (MQ=255)

                                                     aGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgcc                                      <  2:734743/93‑1 (MQ=255)

                                                      gATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgcca                                     >  2:3801958/1‑93 (MQ=255)

                                                          atatCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgccagcca                                 >  2:3020131/1‑93 (MQ=255)

                                                          atatCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgccagcca                                 <  2:1778972/93‑1 (MQ=255)

                                                          atatCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgccagcca                                 <  1:3259306/93‑1 (MQ=255)

                                                          atatCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgccagcca                                 >  1:206498/1‑93 (MQ=255)

                                                           tatCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGccagccag                                >  1:3982421/1‑93 (MQ=255)

                                                            atCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCaga                               <  2:132547/93‑1 (MQ=255)

                                                              cGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGt                             <  1:1075195/93‑1 (MQ=255)

                                                              cGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGt                             <  1:1325236/93‑1 (MQ=255)

                                                              cGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGt                             <  1:2420397/93‑1 (MQ=255)

                                                              cGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGt                             <  1:3352565/93‑1 (MQ=255)

                                                               gTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGtt                            >  1:4407234/1‑93 (MQ=255)

                                                               gTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGtt                            <  2:81421/93‑1 (MQ=255)

                                                               gTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGtt                            <  2:1189427/93‑1 (MQ=255)

                                                               gTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGtt                            <  2:1995550/93‑1 (MQ=255)

                                                                tCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGtta                           >  1:3796988/1‑92 (MQ=255)

                                                                 ccGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGa                          <  2:3251932/93‑1 (MQ=255)

                                                                  cGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGaa                         >  2:4534490/1‑93 (MQ=255)

                                                                    gCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGaaaa                       <  1:2104250/93‑1 (MQ=255)

                                                                    gCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGaaaa                       <  1:3528823/93‑1 (MQ=255)

                                                                    gCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGaaaa                       <  1:3403173/93‑1 (MQ=255)

                                                                      aCCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAAtc                     <  2:2988473/93‑1 (MQ=255)

                                                                      aCCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGACAGCCAGAGTTGAAAAtc                     >  2:4739529/1‑93 (MQ=255)

                                                                        cAGATATTTAAGTTGATCACCCATTGTCGCGAATTCATGCTCGGTAAACGGCACAAATCCCCAGCTGTGCTGCCAGCCAGAGATGAAAAtctc                   >  1:618922/1‑93 (MQ=255)

                                                                         aGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCa                  >  1:884552/1‑93 (MQ=255)

                                                                         aGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCa                  >  1:2707860/1‑93 (MQ=255)

                                                                          gATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCAc                 <  2:1492153/93‑1 (MQ=255)

                                                                           atatTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACg                <  2:1614095/93‑1 (MQ=255)

                                                                            tatTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGc               <  1:336774/93‑1 (MQ=255)

                                                                             atTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCa              <  2:3637850/93‑1 (MQ=255)

                                                                             atTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCa              <  2:1889296/93‑1 (MQ=255)

                                                                              tttAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAg             <  1:551851/93‑1 (MQ=255)

                                                                               ttAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAgg            <  2:247897/93‑1 (MQ=255)

                                                                                tAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGa           <  2:1108189/93‑1 (MQ=255)

                                                                                   gTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCt        >  1:503443/1‑93 (MQ=255)

                                                                                     tGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGc      <  2:4396478/93‑1 (MQ=255)

                                                                                      gATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGCa     <  2:4552238/93‑1 (MQ=255)

GTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGCATCT  >  NC_000913/3127379‑3127560

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