BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,805,574	A→G	T51A (ACT→GCT)	b1724 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,805,574	0	A	G	95.3%	362.0 / 6.2	128	T51A (ACT→GCT)	b1724	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (3/3); new base G (57/65); total (60/68)
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.60e-01

Only 100 of 129 total aligned reads displayed.

CTGAAATTAACGAGTTGAATAATTACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAACTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCGGGGCTTACCCAAAGTGCAACCGGCAAAGA  >  NC_000913/1805485‑1805662

cTGAAATTAACGAGTTGAATAATTACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTa                                                                                       >  1:482236/1‑93 (MQ=255)

 tGAAATTAACGAGTTGAATAATTACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTaa                                                                                      <  1:1003920/93‑1 (MQ=255)

   aaaTTAACGAGTTGAATAATTACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAaca                                                                                    >  1:3007171/1‑93 (MQ=255)

    aaTTAACGAGTTGAATAATTACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAacac                                                                                   >  2:3951767/1‑93 (MQ=255)

    aaTTAACGAGTTGAATAATTACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAacac                                                                                   <  1:4152118/93‑1 (MQ=255)

     aTTAACGAGTTGAATAATTACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACt                                                                                  <  1:2442996/93‑1 (MQ=255)

      ttAACGAGTTGAATAATTACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACtt                                                                                 >  1:1367487/1‑93 (MQ=255)

       tAACGAGTTGAATAATTACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACttt                                                                                >  1:1280575/1‑93 (MQ=255)

        aaCGAGTTGAATAATTACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACtttt                                                                               >  1:1982544/1‑93 (MQ=255)

        aaCGAGTTGAATAATTACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACtttt                                                                               >  1:584462/1‑93 (MQ=255)

        aaCGAGTTGAATAATTACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACtttt                                                                               >  2:2578025/1‑93 (MQ=255)

           gAGTTGAATAATTACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGa                                                                            <  2:1210598/93‑1 (MQ=255)

           gAGTTGAATAATTACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGa                                                                            <  1:4440873/93‑1 (MQ=255)

           gAGTTGAATAATTACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGa                                                                            <  2:2987653/93‑1 (MQ=255)

           gAGTTGAATAATTACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGa                                                                            <  2:2684183/93‑1 (MQ=255)

           gAGTTGAATAATTACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAACTAACACTTTTGGa                                                                            >  2:2180556/1‑93 (MQ=255)

            aGTTGAATAATTACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGAc                                                                           <  1:816184/93‑1 (MQ=255)

             gTTGAATAATTACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACt                                                                          >  1:3198595/1‑93 (MQ=255)

               tGAATAATTACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTa                                                                        <  2:4062010/93‑1 (MQ=255)

                gAATAATTACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTAt                                                                       <  1:4417900/93‑1 (MQ=255)

                gAATAATTACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTAt                                                                       >  2:3413158/1‑93 (MQ=255)

                 aataatTACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATc                                                                      <  2:2463582/93‑1 (MQ=255)

                  ataatTACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCa                                                                     <  2:3665579/93‑1 (MQ=255)

                    aatTACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGt                                                                   >  1:4172865/1‑93 (MQ=255)

                     atTACCTGACTCGCGACGGTTTTTTGCTGACTATGACGGATGATGAGGGCAATATCCATGAGCTGGGAGCTAGCACTTTTGGACTTATCAGTa                                                                  >  2:2751724/1‑93 (MQ=255)

                     atTACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTa                                                                  >  1:4008613/1‑93 (MQ=255)

                      ttACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTAc                                                                 <  1:4130789/93‑1 (MQ=255)

                      ttACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTAc                                                                 <  2:2332443/93‑1 (MQ=255)

                      ttACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTAc                                                                 >  1:4683904/1‑93 (MQ=255)

                      ttACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTAGACTTATCAGTAc                                                                 >  1:2676508/1‑93 (MQ=255)

                        aCCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTAccc                                                               >  2:639018/1‑93 (MQ=255)

                         ccTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCa                                                              >  2:625400/1‑93 (MQ=255)

                          cTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCaa                                                             >  2:2372248/1‑93 (MQ=255)

                           tGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCaaa                                                            >  2:2628437/1‑93 (MQ=255)

                           tGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCaaa                                                            <  2:3319643/93‑1 (MQ=255)

                           tGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTGTTGGACTTATCAGTACCCaaa                                                            <  2:2399176/93‑1 (MQ=255)

                            gACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAg                                                           <  1:168933/93‑1 (MQ=255)

                             aCTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGt                                                          >  1:273021/1‑93 (MQ=255)

                             aCTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGt                                                          <  2:825123/93‑1 (MQ=255)

                              cTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTg                                                         >  2:569253/1‑93 (MQ=255)

                               tCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTga                                                        >  2:1213112/1‑93 (MQ=255)

                               tCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTga                                                        <  2:3191712/93‑1 (MQ=255)

                                  cgACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTgaaga                                                     <  2:3366610/93‑1 (MQ=255)

                                   gACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTgaagaa                                                    >  2:1500985/1‑93 (MQ=255)

                                    aCGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTgaagaag                                                   <  2:4595806/93‑1 (MQ=255)

                                     cGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTgaagaaga                                                  >  2:200838/1‑93 (MQ=255)

                                      ggctttttGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGAgaagaagaa                                                 >  2:4179666/4‑93 (MQ=255)

                                      ggTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTgaagaagaa                                                 >  1:1809437/1‑93 (MQ=255)

                                       gTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGaaa                                                <  2:3484788/93‑1 (MQ=255)

                                        ttttttGCTCATTATGACGGATGATGAGGGAAATATCAATGAGCTGGGAGCTAACACTTCTGGACTTATCAGTTCCCAAAGTGAAGAAGCAAt                                               <  2:1659154/93‑1 (MQ=255)

                                         tttttGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAAtt                                              >  2:300430/1‑93 (MQ=255)

                                          ttttGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTa                                             <  2:4215776/93‑1 (MQ=255)

                                          ttttGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTa                                             <  1:2930091/93‑1 (MQ=255)

                                          ttttGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTa                                             <  2:1587480/93‑1 (MQ=255)

                                           tttGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTag                                            <  2:2939731/93‑1 (MQ=255)

                                            ttGCTTAATAAGACGGATGAATAAGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTaga                                           <  2:3384244/93‑1 (MQ=255)

                                            ttGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTaga                                           <  2:4631795/93‑1 (MQ=255)

                                            ttGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTaga                                           <  1:3448157/93‑1 (MQ=255)

                                              gCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAgaga                                         <  1:3046734/93‑1 (MQ=255)

                                              gCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAgaga                                         <  2:4418481/93‑1 (MQ=255)

                                              gCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAgaga                                         <  1:3504981/93‑1 (MQ=255)

                                               cTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGaa                                        <  1:2069087/93‑1 (MQ=255)

                                                tCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAAc                                       <  1:4322740/93‑1 (MQ=255)

                                                tCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAAc                                       <  1:3594308/93‑1 (MQ=255)

                                                 cACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACt                                      <  1:4049157/93‑1 (MQ=255)

                                                 cACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACt                                      <  1:2681329/93‑1 (MQ=255)

                                                 cACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACt                                      <  1:4332519/93‑1 (MQ=255)

                                                  aCTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTg                                     >  1:1962325/1‑93 (MQ=255)

                                                    tATGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGt                                   >  2:2752345/1‑93 (MQ=255)

                                                     aTGACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGtt                                  <  1:2903879/93‑1 (MQ=255)

                                                       gACGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTc                                >  2:554907/1‑93 (MQ=255)

                                                         cGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCgg                              <  1:863069/93‑1 (MQ=255)

                                                         cGGATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCgg                              <  2:3397896/93‑1 (MQ=255)

                                                          ggATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCggg                             <  2:3627628/93‑1 (MQ=255)

                                                          ggATGATGAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCggg                             >  2:3291971/1‑93 (MQ=255)

                                                            atgatgAGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCGGGGc                           <  2:4081443/93‑1 (MQ=255)

                                                             tgatgaGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCGGGGCt                          <  1:106966/93‑1 (MQ=255)

                                                              gatgaGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCGGGGCtt                         >  2:3615643/1‑93 (MQ=255)

                                                                tgaGGGAAATATCCATGAGCTGGGGGCTAACACTTTTGGACTTATCAGTACCCAAAGTGCAGAAGAAATTAGAGAACTGGGTTCGGGGCTTAc                       >  1:2762079/1‑93 (MQ=255)

                                                                tgaGGGAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCGGGGCTTAc                       <  1:2209413/93‑1 (MQ=255)

                                                                  aGGGAAATATCCATGAGCTGGGAACTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCGGGGCTTAccc                     >  2:2700832/1‑93 (MQ=255)

                                                                   gggAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCGGGGCTTACCCa                    >  1:2201932/1‑93 (MQ=255)

                                                                   gggAAATATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCGGGGCTTACCCa                    <  2:4457207/93‑1 (MQ=255)

                                                                      aaaTATCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCGGGGCTTACCCAAAg                 <  1:4663914/93‑1 (MQ=255)

                                                                        atatCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCGGGGCTTACCCAAAGTg               >  2:1587972/1‑93 (MQ=255)

                                                                        atatCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCGGGGCTTACCCAAAGTg               <  2:4247617/93‑1 (MQ=255)

                                                                           tCCATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCGGGGCTTACCCAAAGTGCaa            >  2:2980369/1‑93 (MQ=255)

                                                                             cATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCGGGGCTTACCCAAAGTGCAAcc          <  1:1980085/93‑1 (MQ=255)

                                                                             cATGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCGGGGCTTACCCAAAGTGCAAcc          >  1:2981996/1‑93 (MQ=255)

                                                                              aTGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCGGGGCTTACCCAAAGTGCAACCg         <  1:2609204/93‑1 (MQ=255)

                                                                              aTGAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCGGGGCTTACCCAAAGTGCAACCg         >  2:2777103/1‑93 (MQ=255)

                                                                                gAGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCGGGGCTTACCCAAAGTGCAACCGGc       <  1:2551505/93‑1 (MQ=255)

                                                                                gAGCTGGGAACTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCGGGGCTTACCCAAAGTGCAACCGGc       <  1:2730038/93‑1 (MQ=255)

                                                                                 aGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCGGGGCTTACCCAAAGTGCAACCGGCa      <  1:3019765/93‑1 (MQ=255)

                                                                                 aGCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCGGGGCTTACCCAAAGTGCAACCGGCa      <  1:438648/93‑1 (MQ=255)

                                                                                  gCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCGGGGCTTACCCAAAGTGCAACCGGCaa     <  1:2988521/93‑1 (MQ=255)

                                                                                  gCTGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCGGGGCTTACCCAAAGTGCAACCGGCaa     >  2:644064/1‑93 (MQ=255)

                                                                                   cTGGTAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCGGGGCTTACCCAAAGTGCAACCGGCaaa    <  1:1728794/93‑1 (MQ=255)

                                                                                    tGGGAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCGGGGCTTACCCAAAGTGCAACCGGCaacg   >  1:3947285/1‑91 (MQ=255)

                                                                                     gggAGCTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCGGGGCTTACCCAAAGTGCAACCGGCAAAGa  >  2:1315721/1‑93 (MQ=255)

CTGAAATTAACGAGTTGAATAATTACCTGACTCGCGACGGTTTTTTGCTCACTATGACGGATGATGAGGGAAATATCCATGAGCTGGGAACTAACACTTTTGGACTTATCAGTACCCAAAGTGAAGAAGAAATTAGAGAACTGGTTTCGGGGCTTACCCAAAGTGCAACCGGCAAAGA  >  NC_000913/1805485‑1805662

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