BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,506,198	G→A	G479E (GGG→GAG)	ydcP →	putative collagenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,506,198	0	G	A	99.6%	963.5 / ‑5.0	250	G479E (GGG→GAG)	ydcP	putative collagenase
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (130/119); total (130/120)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

Only 100 of 252 total aligned reads displayed.

GCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTG  >  NC_000913/1506108‑1506287

gCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGc                                                                                         <  2:6372011/93‑1 (MQ=255)

 cATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCa                                                                                        >  1:5901511/1‑93 (MQ=255)

  ctCACACATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCaa                                                                                       <  2:2566408/92‑1 (MQ=255)

  aTCACGCATCCGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGACCAATCTGAAGGATAGTCTGGCAAAACTGGAGCaa                                                                                       <  1:2485565/93‑1 (MQ=255)

        cATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCctc                                                                                 >  2:476289/1‑93 (MQ=255)

        cATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCctc                                                                                 >  2:4545475/1‑93 (MQ=255)

        cATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATCAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCctc                                                                                 >  1:4539958/1‑93 (MQ=255)

         aTACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCctct                                                                                <  1:6440160/93‑1 (MQ=255)

          tACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTa                                                                               <  1:2282745/93‑1 (MQ=255)

          tACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTa                                                                               <  2:2114576/93‑1 (MQ=255)

           aCGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTAt                                                                              <  2:1043422/93‑1 (MQ=255)

           aCGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTAt                                                                              <  2:4611768/93‑1 (MQ=255)

           aCGCTGGACAGGCAGTACGACGAAGCCAATAACGCCGAAAAAGCAATGAGCAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTAt                                                                              <  2:5597231/93‑1 (MQ=255)

            cGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTAtt                                                                             <  2:3215382/93‑1 (MQ=255)

             cggggACGGGCAGTTCGAGGAAGCCAATAACCCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTa                                                                            <  1:2079414/90‑1 (MQ=255)

               tGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGGAAAACTGGAGCAAACCCTCTATTACg                                                                          >  2:2779854/1‑93 (MQ=255)

               tGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACg                                                                          >  2:4388240/1‑93 (MQ=255)

               tGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACg                                                                          >  1:1167117/1‑93 (MQ=255)

                 gACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGcc                                                                        <  1:2150499/93‑1 (MQ=255)

                 gACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGcc                                                                        >  2:5504510/1‑93 (MQ=255)

                  aCGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGccc                                                                       <  2:1197646/93‑1 (MQ=255)

                  aCGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGccc                                                                       <  1:1703933/93‑1 (MQ=255)

                   cGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCcg                                                                      >  1:1001167/1‑93 (MQ=255)

                    gggCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCcgc                                                                     >  1:4992448/1‑93 (MQ=255)

                    gggCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCcgc                                                                     <  2:802995/93‑1 (MQ=255)

                     ggCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCcgcg                                                                    >  2:3334779/1‑93 (MQ=255)

                      gCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGa                                                                   >  1:5843280/1‑93 (MQ=255)

                         gTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGAtgt                                                                >  2:3453039/1‑93 (MQ=255)

                           tCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGc                                                              <  2:754936/93‑1 (MQ=255)

                            cgacgaAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCa                                                             <  2:5702723/93‑1 (MQ=255)

                              acgaAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCaaa                                                           >  2:1005084/1‑93 (MQ=255)

                               cgaAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAaat                                                          >  2:4073004/1‑93 (MQ=255)

                                gaAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAaatt                                                         >  1:3506475/1‑93 (MQ=255)

                                gaAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGGGCAaatt                                                         >  2:4121873/1‑93 (MQ=255)

                                  aGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAaattaa                                                       >  1:4831371/1‑93 (MQ=255)

                                  aGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAaattaa                                                       >  2:5653608/1‑93 (MQ=255)

                                   gCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAaattaat                                                      >  1:3075934/1‑93 (MQ=255)

                                    ccAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAaattaatt                                                     <  1:6247257/93‑1 (MQ=255)

                                    ccAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAaattaatt                                                     >  1:2638917/1‑93 (MQ=255)

                                     cAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAAttt                                                    <  1:2987414/93‑1 (MQ=255)

                                     cAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAAttt                                                    >  2:2466225/1‑93 (MQ=255)

                                     cAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAAttt                                                    >  2:1890853/1‑93 (MQ=255)

                                     cAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAAttt                                                    >  1:3124177/1‑93 (MQ=255)

                                       aTAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGc                                                  >  1:1425277/1‑93 (MQ=255)

                                       aTAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGc                                                  >  2:3970098/1‑93 (MQ=255)

                                       aTAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGc                                                  <  2:4605161/93‑1 (MQ=255)

                                         aaCGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCg                                                <  2:1518825/93‑1 (MQ=255)

                                         aaCGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCg                                                <  1:1937845/93‑1 (MQ=255)

                                         aaCGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTCGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCg                                                <  2:4462143/93‑1 (MQ=255)

                                          aCGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCgg                                               >  2:3414720/1‑93 (MQ=255)

                                           cGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCggg                                              <  2:1405071/93‑1 (MQ=255)

                                           cGCCGAAAAAGCAATGAACAATCTGAAGCATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCggg                                              <  1:5668158/93‑1 (MQ=255)

                                            gCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCgggg                                             >  1:3894884/1‑93 (MQ=255)

                                              cGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGgc                                           >  1:6342029/1‑93 (MQ=255)

                                              cGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGgc                                           <  1:5289667/93‑1 (MQ=255)

                                              cGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGgc                                           <  2:3361196/93‑1 (MQ=255)

                                              cGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGgc                                           >  1:1386763/1‑93 (MQ=255)

                                              cGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGgc                                           <  1:1536382/93‑1 (MQ=255)

                                                aaaaaGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGgcgc                                         <  1:5476216/93‑1 (MQ=255)

                                                aaaaaGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGgcgc                                         <  2:3494550/93‑1 (MQ=255)

                                                aaaaaGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGgcgc                                         >  2:3012114/1‑93 (MQ=255)

                                                aaaaaGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGgcgc                                         <  2:52723/93‑1 (MQ=255)

                                                 aaaaGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCt                                        >  2:1900219/1‑93 (MQ=255)

                                                 aaaaGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCt                                        <  1:2878979/93‑1 (MQ=255)

                                                  aaaGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTg                                       <  1:4457849/93‑1 (MQ=255)

                                                    aGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGtt                                     <  1:2287363/93‑1 (MQ=255)

                                                     gCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGttt                                    >  1:3720952/1‑93 (MQ=255)

                                                       aaTGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGt                                  >  2:2192495/1‑93 (MQ=255)

                                                        aTGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTa                                 >  2:1751689/1‑93 (MQ=255)

                                                        aTGAACAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTa                                 >  1:1863699/1‑93 (MQ=255)

                                                           aaCAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCa                              >  2:6440274/1‑93 (MQ=255)

                                                           aaCAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCa                              >  1:3374573/1‑93 (MQ=255)

                                                           aaCAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAAATTGCCGGGGGCGCTGTTTGTACCa                              >  2:1284363/1‑93 (MQ=255)

                                                            aCAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTTCCGGGGGCGCTGTTTGTACCaa                             >  1:611161/1‑93 (MQ=255)

                                                            aCAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCaa                             >  2:3175046/1‑93 (MQ=255)

                                                             cAATCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCaaa                            <  1:2413267/93‑1 (MQ=255)

                                                                tCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCAAACAg                         >  2:901890/1‑93 (MQ=255)

                                                                tCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCAAACAg                         >  1:4125599/1‑93 (MQ=255)

                                                                tCTGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCAAACAg                         <  2:3784667/93‑1 (MQ=255)

                                                                  tGAAGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCAAACAGTc                       <  1:4008235/93‑1 (MQ=255)

                                                                   gAAGGATGGCCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCAAACAGTCt                      >  1:4294868/1‑93 (MQ=255)

                                                                     aGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCAAACAGTCTGt                    >  1:760758/1‑93 (MQ=255)

                                                                     aGGATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCAAACAGTCTGt                    <  1:282546/93‑1 (MQ=255)

                                                                      ggATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCAAACAGTCTGtt                   <  2:4323710/93‑1 (MQ=255)

                                                                      ggATGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCAAACAGTCTGtt                   <  1:1661258/93‑1 (MQ=255)

                                                                        aTGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCAAACAGTCTGTTaa                 <  1:470894/93‑1 (MQ=255)

                                                                         tagTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCAAACAGTCTGTTaaa                <  1:3160976/91‑1 (MQ=255)

                                                                         tGGTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGCTTGTACCAAACAGTCTGTTaaa                >  1:5609905/1‑93 (MQ=255)

                                                                           gTCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCAAACAGTCTGTTAAAcc              >  1:1595104/1‑93 (MQ=255)

                                                                            tCTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCa             <  1:6397358/93‑1 (MQ=255)

                                                                             cTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAg            <  2:513746/93‑1 (MQ=255)

                                                                             cTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAg            >  1:4925059/1‑93 (MQ=255)

                                                                             cTGGCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAg            >  1:5654128/1‑93 (MQ=255)

                                                                                gCAAAACTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTc         >  1:2679767/1‑93 (MQ=255)

                                                                                  aaaaCTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTccg       <  1:4767934/93‑1 (MQ=255)

                                                                                  aaaaCTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTccg       >  1:436034/1‑93 (MQ=255)

                                                                                  aaaaCTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTccg       >  2:6168607/1‑93 (MQ=255)

                                                                                   aaaCTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTccgc      >  2:4929241/1‑93 (MQ=255)

                                                                                    aaCTGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTccgcc     >  1:5373399/1‑93 (MQ=255)

                                                                                       tGGAGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTg  >  2:6058659/1‑93 (MQ=255)

GCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTGCCGGGGGCGCTGTTTGTACCAAACAGTCTGTTAAACCAGTTCCGCCGTG  >  NC_000913/1506108‑1506287

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