BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,729,429	G→A	S773S (TCG→TCA)	lhr →	Lhr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,729,429	0	G	A	98.7%	565.8 / ‑2.8	153	S773S (TCG→TCA)	lhr	Lhr
Reads supporting (aligned to +/- strand): ref base G (1/1); new base A (87/64); total (88/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

Only 100 of 153 total aligned reads displayed.

TGACACCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCACTCCGGGCCATCGCACCCCGCTTTGGCA  >  NC_000913/1729339‑1729518

  acacCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCAc                                                                                       >  1:3022612/1‑93 (MQ=255)

   cacCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACt                                                                                      <  2:2650662/93‑1 (MQ=255)

    acCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACtt                                                                                     <  2:423080/93‑1 (MQ=255)

    acCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACtt                                                                                     <  2:3921557/93‑1 (MQ=255)

     ccGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACttt                                                                                    <  2:5213147/93‑1 (MQ=255)

     ccGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACttt                                                                                    <  1:3700794/93‑1 (MQ=255)

      cGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACtttt                                                                                   >  2:290061/1‑93 (MQ=255)

      cGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACtttt                                                                                   <  2:136801/93‑1 (MQ=255)

      cGAGGGTAAATTGCCAGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACtttt                                                                                   <  1:2906505/93‑1 (MQ=255)

       gATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTc                                                                                  >  2:2871268/1‑93 (MQ=255)

       gATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTc                                                                                  <  1:5442497/93‑1 (MQ=255)

       gATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTc                                                                                  >  1:5078020/1‑93 (MQ=255)

       gATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTc                                                                                  <  1:429979/93‑1 (MQ=255)

        aTGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTcg                                                                                 >  1:5237225/1‑93 (MQ=255)

         tGGTAAATTGCCCGATGCGGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTcgc                                                                                <  1:3713986/93‑1 (MQ=255)

         tGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTcgc                                                                                <  1:3515989/93‑1 (MQ=255)

         tGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTcgc                                                                                <  1:269011/93‑1 (MQ=255)

          ggTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTcgcc                                                                               <  1:5705078/93‑1 (MQ=255)

          ggTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTcgcc                                                                               >  2:189993/1‑93 (MQ=255)

           gTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTcgccg                                                                              >  2:4683317/1‑93 (MQ=255)

           gTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTcgccg                                                                              <  2:3576190/93‑1 (MQ=255)

            tAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTcgccgc                                                                             <  1:4687600/93‑1 (MQ=255)

            tAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTcgccgc                                                                             >  1:770506/1‑93 (MQ=255)

             aaaTTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCgccgcc                                                                            <  1:4499053/93‑1 (MQ=255)

              aaTTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGccc                                                                           >  2:1827155/1‑93 (MQ=255)

                 tGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGtt                                                                        <  1:2267127/93‑1 (MQ=255)

                      gATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTcgc                                                                   >  1:2185587/1‑93 (MQ=255)

                      gATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTcgc                                                                   >  1:5361360/1‑93 (MQ=255)

                         gCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGaa                                                                <  1:723499/93‑1 (MQ=255)

                         gCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGaa                                                                <  2:1431618/93‑1 (MQ=255)

                         gCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGaa                                                                >  1:4073776/1‑93 (MQ=255)

                          ccGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGTCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAAt                                                               <  1:3179276/93‑1 (MQ=255)

                          ccGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAAt                                                               <  1:3637637/93‑1 (MQ=255)

                           cgcgATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCACTTTTCGCCGCCCGTTTTCGCGAATg                                                              <  1:4142071/93‑1 (MQ=255)

                           cgcgATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATg                                                              >  1:91390/1‑93 (MQ=255)

                           cgcgATTGTTTTGTGTGATCCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATg                                                              <  1:2929115/93‑1 (MQ=255)

                            gcgATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATgc                                                             <  1:3631638/93‑1 (MQ=255)

                               aTTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGcgcc                                                          >  1:1676374/1‑93 (MQ=255)

                                tttttttGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGcgccg                                                         >  2:5103715/1‑93 (MQ=255)

                                 ttttttGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGcgccgc                                                        >  2:4630434/1‑93 (MQ=255)

                                 ttttttGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGcgccgc                                                        >  1:4534524/1‑93 (MQ=255)

                                  tttttGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCcgcg                                                       >  2:5530559/1‑93 (MQ=255)

                                  tttttGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCCCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATgcgcgggg                                                       >  2:285193/1‑89 (MQ=255)

                                   ttttGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCcgcgc                                                      <  2:2167420/93‑1 (MQ=255)

                                    tttgtttgAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCcgcgcg                                                     >  2:3155008/1‑93 (MQ=255)

                                    tttgtttgAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCcgcgcg                                                     >  1:3164553/1‑93 (MQ=255)

                                    tttgtttgAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCcgcgcg                                                     >  2:2484032/1‑93 (MQ=255)

                                       gtttgAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGc                                                  >  1:2142070/1‑93 (MQ=255)

                                        tttgAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCa                                                 >  2:6064918/1‑93 (MQ=255)

                                           gAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCatta                                              <  2:529235/93‑1 (MQ=255)

                                           gAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCatta                                              >  2:5318103/1‑93 (MQ=255)

                                           gAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCatta                                              <  1:3765362/93‑1 (MQ=255)

                                           gAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCatta                                              <  2:1704865/93‑1 (MQ=255)

                                             aCCAGAAAAGTTGCTGCAAATTGTCGGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCattatt                                            <  1:3966750/93‑1 (MQ=255)

                                             aCCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCattatt                                            >  1:3850338/1‑93 (MQ=255)

                                             aCCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCattatt                                            >  2:2091954/1‑93 (MQ=255)

                                              ccAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCAttatta                                           <  2:4485048/93‑1 (MQ=255)

                                              ccAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCAttatta                                           >  2:1431684/1‑93 (MQ=255)

                                              ccAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCAttatta                                           >  1:1921460/1‑93 (MQ=255)

                                               cAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCCCTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTaa                                          <  1:1859142/93‑1 (MQ=255)

                                               cAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTaa                                          >  2:5728775/1‑93 (MQ=255)

                                               cAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTaa                                          <  2:4835062/93‑1 (MQ=255)

                                                aGAAAAGGTGCTGCAAATTGTCCGGGAGGCGGTAGGCAGCTCAGCACTTTGCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAAt                                         <  2:1549408/93‑1 (MQ=255)

                                                  aaaaGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGc                                       <  1:2444610/93‑1 (MQ=255)

                                                  aaaaGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGc                                       >  2:2003467/1‑93 (MQ=255)

                                                   aaaGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGcc                                      >  1:5635809/1‑93 (MQ=255)

                                                     aGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCgg                                    <  1:2954543/93‑1 (MQ=255)

                                                      gTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCgcg                                   >  2:2480776/1‑91 (MQ=255)

                                                       ttGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCgggg                                  >  1:4942456/1‑93 (MQ=255)

                                                        tgctgcAAATTGTCCTCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGgc                                 <  1:3140817/93‑1 (MQ=255)

                                                          ctgcAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGgcgc                               <  1:1359381/93‑1 (MQ=255)

                                                          ctgcAAATTGTCCGCGAGGCGGTAGGCAGCCCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGgcgc                               >  1:5923917/1‑93 (MQ=255)

                                                            gcAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCAc                             >  2:2351725/1‑93 (MQ=255)

                                                            gcAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCAc                             >  2:1260022/1‑93 (MQ=255)

                                                            gcAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCAc                             <  1:1530372/93‑1 (MQ=255)

                                                             gaaaTTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCACt                            <  1:1898349/92‑1 (MQ=255)

                                                             cAAATTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCACt                            >  1:3408927/1‑93 (MQ=255)

                                                              aaaTTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCACTc                           <  1:2250541/93‑1 (MQ=255)

                                                              aaaTTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCACTc                           >  2:3837371/1‑93 (MQ=255)

                                                               aaTTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCACTcc                          <  1:5910350/93‑1 (MQ=255)

                                                               aaTTGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCACTcc                          >  1:4691711/1‑93 (MQ=255)

                                                                 ttGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCACTCCgg                        >  2:3138427/1‑93 (MQ=255)

                                                                 tcgTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCACTCCgg                        <  1:1192746/91‑1 (MQ=255)

                                                                  tGTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCACTCCggg                       >  1:4987819/1‑93 (MQ=255)

                                                                   gTCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCACTCCGGGc                      <  2:4805500/93‑1 (MQ=255)

                                                                    tCCGCGAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCACTCCGGGcc                     >  2:6117223/1‑93 (MQ=255)

                                                                      cgcgAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCACTCCGGGCCAt                   >  2:291181/1‑93 (MQ=255)

                                                                      cgcgAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCACTCCGGGCCAt                   <  2:6042130/93‑1 (MQ=255)

                                                                         gAGGCGGTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCACTCCGCGCCATCga                >  2:3948596/1‑92 (MQ=255)

                                                                              ggTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCACTCCGGGCCATCGCAcccc           <  1:6064918/93‑1 (MQ=255)

                                                                               gTAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCACTCCGGGCCATCGCACCCCg          >  1:4145515/1‑93 (MQ=255)

                                                                                tAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCACTCCGGGCCATCGCACCCCGc         >  2:2771712/1‑93 (MQ=255)

                                                                                tAGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCACTCCGGGCCATCGCACCCCGc         >  1:1467669/1‑93 (MQ=255)

                                                                                 aGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCACTCCGGGCCATCGCACCCCGCt        >  1:3741984/1‑93 (MQ=255)

                                                                                 aGGCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCACTCCGGGCCATCGCACCCCGCt        >  1:2363962/1‑93 (MQ=255)

                                                                                   gCAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCACTCCGGGCCATCGCACCCCGCttt      <  1:5739279/93‑1 (MQ=255)

                                                                                    cAGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCACTCCGGGCCATCGCACCCCGCTTTg     <  2:4423277/93‑1 (MQ=255)

                                                                                     aGCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCACTCCGGGCCATCGCACCCCGCTTTgg    >  1:4294575/1‑93 (MQ=255)

                                                                                      gCTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCACTCCGGGCCATCGCACCCCGCTTTGGc   >  2:2265396/1‑93 (MQ=255)

                                                                                       cTCAGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCACTCCGGGCCATCGCACCCCGCTTTGGca  <  1:2824678/93‑1 (MQ=255)

TGACACCGATGGTAAATTGCCCGATGCCGCGATTTTTTTGTTTGAACCAGAAAAGTTGCTGCAAATTGTCCGCGAGGCGGTAGGCAGCTCGGCACTTTTCGCCGCCCGTTTTCGCGAATGCGCCGCGCGGGCATTATTAATGCCGGGGCGCACTCCGGGCCATCGCACCCCGCTTTGGCA  >  NC_000913/1729339‑1729518

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