BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,388,615	A→G	V158A (GTG→GCG)	b4161 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,388,615	0	A	G	98.6%	964.4 / 1.1	275	V158A (GTG→GCG)	b4161	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/3); new base G (148/123); total (149/126)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

Only 100 of 287 total aligned reads displayed.

AATACGCGATAACCGATATTGCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGCAAAATGGCGGAGACA  >  NC_000913/4388525‑4388704

aaTACGCGATAACCGATATTGCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGcgca                                                                                         >  1:2940471/1‑93 (MQ=255)

 aTACGCGATAACCGATATTGCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCaa                                                                                        <  1:682975/93‑1 (MQ=255)

 aTACGCGATAACCGATATTGCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCaa                                                                                        <  1:3241301/93‑1 (MQ=255)

 aTACGCGATAACCGATATTGCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCaa                                                                                        <  1:2054166/93‑1 (MQ=255)

  tACGCGATAACCGATATTGCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCaat                                                                                       <  1:767726/93‑1 (MQ=255)

    cgcgATAACCGATATTGCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCaataa                                                                                     >  1:2016834/1‑93 (MQ=255)

     gcgATAACCGATATTGCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCaataat                                                                                    <  2:1171837/93‑1 (MQ=255)

      cgATAACCGATATTGCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATc                                                                                   <  1:3769415/93‑1 (MQ=255)

      cgATAACCGATATTGCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCACCAGGTCTATCTTGTTGAGCGCAATAATc                                                                                   >  1:2146755/1‑93 (MQ=255)

       gATAACCGATATTGCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCg                                                                                  >  2:3921148/1‑93 (MQ=255)

        aTAACCGATATTGCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCgg                                                                                 <  2:2634721/93‑1 (MQ=255)

           aCCGATATTGCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTc                                                                              >  2:445614/1‑93 (MQ=255)

            ccGATATTGCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCa                                                                             <  1:2820710/93‑1 (MQ=255)

             cGATATTGCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCaa                                                                            >  1:1929371/1‑93 (MQ=255)

              gATATTGCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAAt                                                                           >  1:4131830/1‑93 (MQ=255)

              gATATTGCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAAt                                                                           >  1:5683484/1‑93 (MQ=255)

              gATATTGCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAAt                                                                           >  1:596034/1‑93 (MQ=255)

              gATATTGCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAAt                                                                           >  2:821640/1‑93 (MQ=255)

               atatTGCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATc                                                                          <  1:543402/93‑1 (MQ=255)

                tatTGCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCt                                                                         >  2:1669350/1‑93 (MQ=255)

                 atTGCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTg                                                                        <  1:3700743/93‑1 (MQ=255)

                  ttGCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGc                                                                       >  2:3792566/1‑93 (MQ=255)

                   tGCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCa                                                                      <  2:214676/93‑1 (MQ=255)

                    gCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCaa                                                                     <  1:139708/93‑1 (MQ=255)

                    gCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCaa                                                                     >  2:2387514/1‑93 (MQ=255)

                     cGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAg                                                                    <  1:882998/93‑1 (MQ=255)

                     cGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAg                                                                    >  1:3974552/1‑93 (MQ=255)

                      ggTAGATATCCATCTGCTCGTTGCAGAACGCCCTGCCTTCGTCGTCCAGCAGGCCTACCTTTTTTAGCGCAATACTCGGCCCAATCTGCAAgg                                                                   >  1:2949904/1‑93 (MQ=255)

                      ggTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCTGCTCAATGTGCAAgg                                                                   >  2:1782749/1‑93 (MQ=255)

                       gTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGt                                                                  <  2:3095058/93‑1 (MQ=255)

                        tAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGTGCGCAATAATCGGCTCAATCTGCAAGGtt                                                                 >  2:556750/1‑93 (MQ=255)

                        tAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGtt                                                                 <  2:1824062/93‑1 (MQ=255)

                           atatCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCg                                                              >  2:3362423/1‑93 (MQ=255)

                            tatCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGc                                                             >  2:1547314/1‑93 (MQ=255)

                            tatCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGc                                                             >  1:221831/1‑93 (MQ=255)

                             atCCATCTCCGCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCa                                                            <  1:4568486/93‑1 (MQ=255)

                                cATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGc                                                         <  2:1462362/93‑1 (MQ=255)

                                cATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGc                                                         <  2:1312730/93‑1 (MQ=255)

                                  tCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGcca                                                       >  2:2469002/1‑93 (MQ=255)

                                  tCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGcca                                                       >  1:5105535/1‑93 (MQ=255)

                                  tCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGcca                                                       >  2:64863/1‑93 (MQ=255)

                                  tCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGcca                                                       >  2:3222549/1‑93 (MQ=255)

                                   cTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGGCTATCTTGTTTTGCGCAATAATCGGCTCTATCTGCAAGGGTTCGCAGGccac                                                      >  1:1489708/1‑93 (MQ=255)

                                    tGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGccacc                                                     >  1:1364519/1‑93 (MQ=255)

                                    tGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGccacc                                                     >  1:2324655/1‑93 (MQ=255)

                                    tGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGccacc                                                     >  1:3007113/1‑93 (MQ=255)

                                    tGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGccacc                                                     >  2:4939661/1‑93 (MQ=255)

                                    tGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGccacc                                                     >  1:4125181/1‑93 (MQ=255)

                                    tGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGccacc                                                     <  1:3596320/93‑1 (MQ=255)

                                     gCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGcgacca                                                    >  1:4907122/1‑93 (MQ=255)

                                     gCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGccacca                                                    <  1:1183527/93‑1 (MQ=255)

                                     gCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGccacca                                                    >  1:4510736/1‑93 (MQ=255)

                                      cTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAg                                                   <  2:4197450/93‑1 (MQ=255)

                                        cGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGt                                                 <  1:5361702/93‑1 (MQ=255)

                                         gTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTa                                                <  2:2806124/93‑1 (MQ=255)

                                          ttGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTaa                                               <  2:3622578/93‑1 (MQ=255)

                                          ttGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTaa                                               >  1:2697284/1‑93 (MQ=255)

                                          ttGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTaa                                               >  2:5761458/1‑93 (MQ=255)

                                           tGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAAc                                              <  1:2503496/93‑1 (MQ=255)

                                           tGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAAc                                              <  1:4455553/93‑1 (MQ=255)

                                            gACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACg                                             <  1:367322/93‑1 (MQ=255)

                                              cGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGt                                           <  2:5367489/93‑1 (MQ=255)

                                               gAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTc                                          <  1:1342643/93‑1 (MQ=255)

                                                aaCGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCg                                         >  1:604834/1‑93 (MQ=255)

                                                    ccATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGataa                                     <  1:5200442/93‑1 (MQ=255)

                                                    ccATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGataa                                     <  2:4859301/93‑1 (MQ=255)

                                                     cATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGataat                                    <  2:5104475/93‑1 (MQ=255)

                                                        gCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATAtt                                 >  2:2246460/1‑93 (MQ=255)

                                                        gCCGTCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATAtt                                 <  1:1878091/93‑1 (MQ=255)

                                                           ttCGTCGTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAg                              <  1:4715854/93‑1 (MQ=255)

                                                            tcgtcgTCCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGc                             >  2:4577452/1‑93 (MQ=255)

                                                               tcgtcCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGAc                          <  2:5603168/93‑1 (MQ=255)

                                                               tcgtcCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGAc                          <  2:5806218/93‑1 (MQ=255)

                                                               tcgtcCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGAc                          >  1:4809459/1‑93 (MQ=255)

                                                                  tcCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGc                       >  2:2287918/1‑93 (MQ=255)

                                                                  tcCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGc                       <  2:1309817/93‑1 (MQ=255)

                                                                  tcCAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGc                       >  1:2421725/1‑93 (MQ=255)

                                                                   ccAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCt                      >  2:4011689/1‑93 (MQ=255)

                                                                   ccAGCAGGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCt                      >  1:1009816/1‑93 (MQ=255)

                                                                    cagcagGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTc                     <  1:167919/93‑1 (MQ=255)

                                                                    cagcagGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTc                     <  2:4290061/93‑1 (MQ=255)

                                                                        agGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGc                 >  2:6060929/1‑93 (MQ=255)

                                                                        agGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGc                 >  1:5819101/1‑93 (MQ=255)

                                                                        agGTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGc                 >  1:3417597/1‑93 (MQ=255)

                                                                        aagTCTATCTTGTTGAGCGCAATAATCGGGTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGc                 >  1:4357963/3‑93 (MQ=255)

                                                                          gTCTATCTTGTTGAGCGCAATATTCGGCTCAATCTGCAAGGATTCGCAGGCCTCCAGGTAACGGTTGTTAATATTGAGCGACTGCTCCGGCaa               >  1:5463190/1‑93 (MQ=255)

                                                                          gTCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGCaa               >  1:170435/1‑93 (MQ=255)

                                                                           tCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGCaaa              >  2:4327012/1‑93 (MQ=255)

                                                                           tCTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGCaaa              >  1:802415/1‑93 (MQ=255)

                                                                            cTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGCaaaa             >  1:2628115/1‑93 (MQ=255)

                                                                            cTATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGCaaaa             <  1:1829299/93‑1 (MQ=255)

                                                                             tATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGCAAAAt            >  2:5183349/1‑93 (MQ=255)

                                                                             tATCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGCAAAAt            >  2:5239788/1‑93 (MQ=255)

                                                                             tATCTTGTTGAGCGAAATAATCGGCTCAATCTGCAAGGTTTCGCAGTCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGCAAAAt            <  1:100090/93‑1 (MQ=255)

                                                                              aTCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGCAAAATg           <  1:3192521/93‑1 (MQ=255)

                                                                               tCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGCAAAATgg          <  1:5500784/93‑1 (MQ=255)

                                                                               tCTTGTTGAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGCAAAATgg          <  1:4920059/93‑1 (MQ=255)

                                                                                 ttgttgAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGCAAAATGGCg        >  2:1272052/1‑93 (MQ=255)

                                                                                  tgttgAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGCAAAATGGCgg       >  2:1095174/1‑93 (MQ=255)

                                                                                   gttgAGCGCAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGCAAAATGGCGga      >  1:2167722/1‑93 (MQ=255)

AATACGCGATAACCGATATTGCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGCAAAATGGCGGAGACA  >  NC_000913/4388525‑4388704

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: `‑`