BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,679,869	A→G	S62S (AGT→AGC)	ydgC ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,679,869	0	A	G	86.3%	689.2 / 92.1	284	S62S (AGT→AGC)	ydgC	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (19/20); new base G (121/124); total (140/144)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

Only 100 of 287 total aligned reads displayed.

ACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAAT  >  NC_000913/1679777‑1679961

aCAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATa                                                                                              >  2:3278977/1‑93 (MQ=255)

    aaCCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTaa                                                                                          <  1:2324917/93‑1 (MQ=255)

       caacaaAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTaaaaa                                                                                       <  2:2270692/93‑1 (MQ=255)

        aacaaAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAat                                                                                      >  2:4025161/1‑93 (MQ=255)

          caaAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAatga                                                                                    <  1:6881971/93‑1 (MQ=255)

           aaaaGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAatgat                                                                                   >  1:6627638/1‑93 (MQ=255)

           aaaaGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAatgat                                                                                   >  2:3170771/1‑93 (MQ=255)

            aaaGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAatgatg                                                                                  <  1:1764532/93‑1 (MQ=255)

             aaGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATgg                                                                                 <  2:2869752/93‑1 (MQ=255)

                cggcggGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTg                                                                              >  1:3476608/1‑93 (MQ=255)

                cggcggGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTg                                                                              <  1:1402184/93‑1 (MQ=255)

                 ggcggGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTgc                                                                             <  2:2006803/93‑1 (MQ=255)

                 ggcggGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTgc                                                                             <  2:3106400/93‑1 (MQ=255)

                    gggCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGcgcg                                                                          >  1:2342476/1‑93 (MQ=255)

                    gggCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGcgcg                                                                          >  1:501296/1‑93 (MQ=255)

                    gggCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGcgcg                                                                          >  2:5813856/1‑93 (MQ=255)

                    gggCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGcgcg                                                                          <  1:3229516/93‑1 (MQ=255)

                     ggCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGt                                                                         >  1:1321408/1‑93 (MQ=255)

                        aaGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAg                                                                      <  1:4825226/93‑1 (MQ=255)

                        aaGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAg                                                                      <  1:255663/93‑1 (MQ=255)

                         aGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAgg                                                                     >  2:7035600/1‑93 (MQ=255)

                           cgcATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCt                                                                   <  2:686668/93‑1 (MQ=255)

                             catcatTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTc                                                                 >  2:2755588/1‑93 (MQ=255)

                              atcatTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCa                                                                >  2:5448862/1‑93 (MQ=255)

                               tcatTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCaa                                                               >  1:1036404/1‑93 (MQ=255)

                               tcatTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCaa                                                               <  2:1951548/93‑1 (MQ=255)

                               tcatTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCaa                                                               <  1:5399023/93‑1 (MQ=255)

                               tcatTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCaa                                                               >  2:6231314/1‑93 (MQ=255)

                                catTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAAt                                                              <  1:6120679/93‑1 (MQ=255)

                                 atTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATg                                                             <  2:5908316/93‑1 (MQ=255)

                                  ttCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGc                                                            >  2:164084/1‑93 (MQ=255)

                                     cGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCcgc                                                         >  1:6104073/1‑93 (MQ=255)

                                       gTGAAATACCACAGCGACACCAGGTAAAAAAAATAGGGAATAATCGACCACATACTAAAAATGATGGTTGCGCGAAAGGCTTCAATGCCGCGt                                                       <  1:3110874/93‑1 (MQ=255)

                                        tGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGtt                                                      <  1:26688/93‑1 (MQ=255)

                                        tGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGtt                                                      <  2:5601739/93‑1 (MQ=255)

                                        tGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGtt                                                      <  1:7120252/93‑1 (MQ=255)

                                        tGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGtt                                                      <  1:2155791/93‑1 (MQ=255)

                                           aaTACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGc                                                   <  2:178678/93‑1 (MQ=255)

                                           aaTACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGc                                                   >  1:1527017/1‑93 (MQ=255)

                                            aTACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCt                                                  <  1:3126306/93‑1 (MQ=255)

                                             tACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTg                                                 <  1:4056709/93‑1 (MQ=255)

                                             tACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTg                                                 >  1:6997139/1‑93 (MQ=255)

                                              aCCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTgg                                                <  1:7154812/93‑1 (MQ=255)

                                              aCCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTgg                                                >  2:3693907/1‑93 (MQ=255)

                                                 acaGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGcaa                                             <  2:6672552/93‑1 (MQ=255)

                                                  caGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGcaac                                            >  2:2389477/1‑93 (MQ=255)

                                                   aGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGcaaca                                           >  1:749575/1‑93 (MQ=255)

                                                    gCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGcaacaa                                          >  1:2139697/1‑93 (MQ=255)

                                                      gACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAata                                        <  1:2743759/93‑1 (MQ=255)

                                                        cacCAGGTAAACAAAATACGGAATAATCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAAtata                                      <  2:2653838/93‑1 (MQ=255)

                                                         acCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATaa                                     <  2:5856478/93‑1 (MQ=255)

                                                         acCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATaa                                     >  2:3829148/1‑93 (MQ=255)

                                                          ccAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAAt                                    <  1:7476008/93‑1 (MQ=255)

                                                          ccAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAAt                                    <  1:3533493/93‑1 (MQ=255)

                                                          ccAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAAt                                    >  2:5309820/1‑93 (MQ=255)

                                                          ccAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAAt                                    >  2:1426718/1‑93 (MQ=255)

                                                          ccAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAAt                                    >  2:4712328/1‑93 (MQ=255)

                                                           cAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATg                                   <  2:5125327/93‑1 (MQ=255)

                                                           cAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATg                                   >  1:5179651/1‑93 (MQ=255)

                                                            aGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATgc                                  <  1:3845505/93‑1 (MQ=255)

                                                               tAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGcgcg                               <  2:2542267/93‑1 (MQ=255)

                                                                aaaCAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGa                              <  2:5483220/93‑1 (MQ=255)

                                                                   cAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATaa                           <  2:4963380/93‑1 (MQ=255)

                                                                   cAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATaa                           <  1:895869/93‑1 (MQ=255)

                                                                   cAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATaa                           >  1:5928917/1‑93 (MQ=255)

                                                                   cAAAATAGGGAATAATCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATaa                           <  2:4251915/93‑1 (MQ=255)

                                                                    aaaaTAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAg                          >  1:1920228/1‑93 (MQ=255)

                                                                    aaaaTAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAg                          <  2:3556620/93‑1 (MQ=255)

                                                                       aTAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAgcgc                       >  2:3625407/1‑93 (MQ=255)

                                                                        tAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCa                      >  2:484975/1‑93 (MQ=255)

                                                                        tAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCa                      >  1:1280229/1‑93 (MQ=255)

                                                                          gggAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCaaa                    >  2:6360845/1‑93 (MQ=255)

                                                                          gggAATAATCGACCACATGCAAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCGGGCACCAATATAATGCGCGATAAGCGCaaa                    <  1:3670548/93‑1 (MQ=255)

                                                                           ggAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAg                   >  1:5719420/1‑93 (MQ=255)

                                                                           ggAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAg                   >  2:6495927/1‑93 (MQ=255)

                                                                           ggAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAg                   <  1:4220661/93‑1 (MQ=255)

                                                                           ggAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAg                   <  1:4412579/93‑1 (MQ=255)

                                                                            gAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAgg                  <  1:3337706/93‑1 (MQ=255)

                                                                            gAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAgg                  <  2:4647449/93‑1 (MQ=255)

                                                                            gAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAgg                  <  2:7283985/93‑1 (MQ=255)

                                                                            gAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAgg                  <  1:5053025/93‑1 (MQ=255)

                                                                             aataatCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGt                 <  1:2520134/93‑1 (MQ=255)

                                                                               taatCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCg               <  1:564562/93‑1 (MQ=255)

                                                                               taatCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCg               <  2:104480/93‑1 (MQ=255)

                                                                                  tCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGa            >  2:1533080/1‑93 (MQ=255)

                                                                                  tCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGa            >  2:3441024/1‑93 (MQ=255)

                                                                                    gACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGaaa          >  2:533394/1‑93 (MQ=255)

                                                                                     aCCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGaaaa         >  2:7390678/1‑93 (MQ=255)

                                                                                     aCCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGaaaa         >  2:1845870/1‑93 (MQ=255)

                                                                                     aCCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGaaaa         >  2:4497346/1‑93 (MQ=255)

                                                                                      ccACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAg        >  1:6015559/1‑93 (MQ=255)

                                                                                      ccACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAg        >  2:3495682/1‑93 (MQ=255)

                                                                                       cacaTGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGt       >  1:5080991/1‑93 (MQ=255)

                                                                                       cacaTACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGt       >  1:7233529/1‑93 (MQ=255)

                                                                                        acaTGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTg      <  2:6005704/93‑1 (MQ=255)

                                                                                        acaTGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTg      <  2:2939098/93‑1 (MQ=255)

                                                                                        acaTGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTg      <  1:1003082/93‑1 (MQ=255)

                                                                                        acaTGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATACCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTg      <  2:2165966/93‑1 (MQ=255)

                                                                                         caTGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTgg     >  1:2727859/1‑93 (MQ=255)

                                                                                            aCTAAAAATGATGGGTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAAt  >  2:7485584/1‑93 (MQ=255)

ACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAAT  >  NC_000913/1679777‑1679961

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