BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,840,766	G→A	C243C (TGC→TGT)	yicO ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,840,766	0	G	A	99.7%	1306.1 / NA	365	C243C (TGC→TGT)	yicO	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (197/167); minor base C (0/1); total (197/168)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 397 total aligned reads displayed.

GCTCAAATCTACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAGCAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAACCCTAAAATACCTAA  >  NC_000913/3840676‑3840855

 cTCAAATCTACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCAcaacaa                                                                                        >  2:4279484/1‑93 (MQ=255)

 cTCAAATCTACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCAcaacaa                                                                                        <  2:7604343/93‑1 (MQ=255)

  tCAAATCTACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCAcaacaaa                                                                                       >  2:158809/1‑92 (MQ=255)

   cAAATCTACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGa                                                                                      <  2:3150230/93‑1 (MQ=255)

     aaTCTACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACg                                                                                    >  2:3792171/1‑93 (MQ=255)

       tCTACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTc                                                                                  <  1:4310797/93‑1 (MQ=255)

       tCTACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTc                                                                                  <  2:5311588/93‑1 (MQ=255)

       tCTACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTc                                                                                  >  2:6993537/1‑93 (MQ=255)

       tCTACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTc                                                                                  >  1:6396111/1‑93 (MQ=255)

       tCTACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTc                                                                                  <  2:2957315/93‑1 (MQ=255)

         tACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTcac                                                                                >  1:5672417/1‑93 (MQ=255)

         tACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTcac                                                                                >  2:8579378/1‑93 (MQ=255)

         tACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAACGAACATCACCGAAAAATAATCCACAACAAGACGTcac                                                                                <  1:96250/93‑1 (MQ=255)

         tACGTCACCAATGACGCCGCTAATATCAGGAGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTcac                                                                                <  2:8084925/93‑1 (MQ=255)

          aCTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTcacc                                                                               >  1:5121226/1‑93 (MQ=255)

           cTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTcacca                                                                              <  1:5741636/93‑1 (MQ=255)

            ttCACCAATGACGCCGCTAATTTCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTcaccaa                                                                             >  2:7139301/1‑92 (MQ=255)

             tCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAATAAGACGTCACcaca                                                                            >  2:945148/1‑93 (MQ=255)

              cACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACaa                                                                           >  1:2506205/1‑93 (MQ=255)

              cACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACaa                                                                           <  1:1985509/93‑1 (MQ=255)

               aCCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAAt                                                                          >  2:8534094/1‑93 (MQ=255)

                ccAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATa                                                                         >  2:6972305/1‑93 (MQ=255)

                 cAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAg                                                                        >  2:5106118/1‑93 (MQ=255)

                  aaTGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAGCATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGa                                                                       >  2:7494503/1‑93 (MQ=255)

                  aaTGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGa                                                                       <  2:6500503/93‑1 (MQ=255)

                  aaTGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGa                                                                       <  2:333715/93‑1 (MQ=255)

                    tGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGaaa                                                                     <  2:5948751/93‑1 (MQ=255)

                    tGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGaaa                                                                     >  2:558914/1‑93 (MQ=255)

                    tGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGaaa                                                                     >  2:388880/1‑93 (MQ=255)

                     gACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAAc                                                                    >  1:1741195/1‑93 (MQ=255)

                           gcTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCAc                                                              >  2:1141254/1‑93 (MQ=255)

                           gcTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCAc                                                              <  1:6550853/93‑1 (MQ=255)

                            cTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCAcc                                                             <  2:6653505/93‑1 (MQ=255)

                             tAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACcg                                                            <  2:7734477/93‑1 (MQ=255)

                             tAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACcg                                                            <  2:2802794/93‑1 (MQ=255)

                              aaTATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACcgc                                                           >  1:7291577/1‑93 (MQ=255)

                               atatCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACcgcg                                                          <  1:7174154/93‑1 (MQ=255)

                                 atCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATATAAACCAGCACCGCGGc                                                        >  1:7787268/1‑93 (MQ=255)

                                   cAGGCGGAATGGAATAGACGCCGCTAAAACGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCAt                                                      >  2:17398/1‑93 (MQ=255)

                                       cGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGaaa                                                  <  2:3442481/93‑1 (MQ=255)

                                       cGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGaaa                                                  >  1:8761826/1‑93 (MQ=255)

                                            tggaatAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGAc                                             <  1:5269783/93‑1 (MQ=255)

                                             ggaatAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACg                                            <  2:8188398/93‑1 (MQ=255)

                                             ggaatAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACg                                            <  2:1186820/93‑1 (MQ=255)

                                             ggaatAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACg                                            <  1:1187357/93‑1 (MQ=255)

                                             ggaatAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACg                                            <  1:7054795/93‑1 (MQ=255)

                                             ggaatAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACg                                            <  1:474050/93‑1 (MQ=255)

                                              gaatAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGt                                           >  2:1023656/1‑93 (MQ=255)

                                              aaatAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCGCAACAAGACGTCGCCACAATAGAAACCAGCACCGCGGCATGAAAATGACGt                                           <  1:2604926/92‑1 (MQ=255)

                                               aatAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGtg                                          >  1:3067885/1‑93 (MQ=255)

                                                atAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGtga                                         >  1:7604872/1‑93 (MQ=255)

                                                atAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGtga                                         >  2:6195607/1‑93 (MQ=255)

                                                atAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGtga                                         <  1:2241338/93‑1 (MQ=255)

                                                 tAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGtgat                                        >  2:4540899/1‑93 (MQ=255)

                                                 tAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGtgat                                        <  2:195004/93‑1 (MQ=255)

                                                 tAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGtgat                                        >  2:1318102/1‑93 (MQ=255)

                                                    aCGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGac                                     >  2:1375707/1‑93 (MQ=255)

                                                    aCGCCGCTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGac                                     <  1:130719/93‑1 (MQ=255)

                                                     cgccgcTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGaca                                    >  1:3710943/1‑93 (MQ=255)

                                                      gccgcTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGacaa                                   <  1:3915907/93‑1 (MQ=255)

                                                       ccgcTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGacaac                                  <  1:8489492/93‑1 (MQ=255)

                                                          cTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACg                               >  2:5812872/1‑93 (MQ=255)

                                                          cTAAAATGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACg                               >  1:1042116/1‑93 (MQ=255)

                                                            aaaaTGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGGAACGCGGCATGAAAATGAGGTAATGACAACACGGt                             >  1:5949075/1‑93 (MQ=255)

                                                             aaaTGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGtt                            >  1:325067/1‑93 (MQ=255)

                                                              aaTGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTta                           <  2:3416312/93‑1 (MQ=255)

                                                              aaTGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTta                           <  2:3316827/93‑1 (MQ=255)

                                                              aaTGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTta                           >  2:2357730/1‑93 (MQ=255)

                                                               aTGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTtat                          <  1:6675626/93‑1 (MQ=255)

                                                                tGAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTtata                         <  1:726690/93‑1 (MQ=255)

                                                                 gAACATCACCGAAAAATAATCCACAACAAGACGTCCCCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGCCAACACGGTTataa                        >  2:1792449/1‑93 (MQ=255)

                                                                 gAACATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTataa                        <  1:4529218/93‑1 (MQ=255)

                                                                    cATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAACATGACGTGATGACAACACGGTTAtaataa                     >  1:161119/1‑93 (MQ=255)

                                                                    cATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTAtaataa                     >  1:3391581/1‑93 (MQ=255)

                                                                    cATCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTAtaataa                     >  2:4864126/1‑93 (MQ=255)

                                                                      tCACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATaaaa                   >  2:6066541/1‑93 (MQ=255)

                                                                       cACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATaaaaa                  <  2:1240519/93‑1 (MQ=255)

                                                                       cACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATaaaaa                  >  2:2628941/1‑93 (MQ=255)

                                                                       cACCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATaaaaa                  <  2:7883300/93‑1 (MQ=255)

                                                                        aCCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATaaaaaa                 <  2:4713541/93‑1 (MQ=255)

                                                                        aCCGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATaaaaaa                 >  1:2369293/1‑93 (MQ=255)

                                                                         ccGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAAc                <  2:3698583/93‑1 (MQ=255)

                                                                          cGAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAAcc               <  1:1845316/93‑1 (MQ=255)

                                                                           gAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAAccc              >  2:5471105/1‑93 (MQ=255)

                                                                           gAAAAATAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAAccc              >  2:7526606/1‑93 (MQ=255)

                                                                             aaaaTAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAACCCTa            >  2:1090510/1‑93 (MQ=255)

                                                                              aaaTAATCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAACCCTaa           <  2:2168617/93‑1 (MQ=255)

                                                                               aataatCCAGAACAAGACGTCACCACAGTAGAAACCAGCACCGCGGCATGAAAATGCCGTGATGCCAACACGGTTATAATAAAAAACCCTaaa          <  1:3158140/93‑1 (MQ=255)

                                                                               aataatCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAACCCTaaa          >  1:80682/1‑93 (MQ=255)

                                                                               aataatCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAACCCTaaa          >  1:3966841/1‑93 (MQ=255)

                                                                               aataatCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAACCCTaaa          <  2:5191415/93‑1 (MQ=255)

                                                                                ataatCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAACCCTaaaa         >  1:7851985/1‑93 (MQ=255)

                                                                                ataatCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAACCCTaaaa         <  1:6979209/93‑1 (MQ=255)

                                                                                 taatCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGCCAACACGGTTATAATAAAAAACCCTAAAAt        >  2:7266097/1‑93 (MQ=255)

                                                                                 taatCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAACCCTAAAAt        >  2:246859/1‑93 (MQ=255)

                                                                                  aagccACAACAAGACGTCACCATAATAGAAACCAGCCCCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAACCCTAAAATa       <  1:1903536/90‑1 (MQ=255)

                                                                                   atCCACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAACCCTAAAATAc      <  2:4018472/93‑1 (MQ=255)

                                                                                     ccACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAACCCTAAAATACCt    <  1:6155507/93‑1 (MQ=255)

                                                                                     ccACAACAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAACCCTAAAATACCt    >  1:5460929/1‑93 (MQ=255)

                                                                                       acaacaAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAACCCTAAAATACCTaa  <  2:6894107/93‑1 (MQ=255)

GCTCAAATCTACTTCACCAATGACGCCGCTAATATCAGGCGGAATGGAATAGACGCCGCTAAAATGAACATCACCGAAAAATAATCCACAGCAAGACGTCACCACAATAGAAACCAGCACCGCGGCATGAAAATGACGTGATGACAACACGGTTATAATAAAAAACCCTAAAATACCTAA  >  NC_000913/3840676‑3840855

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