BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,117,911	A→G	V462A (GTA→GCA)	yghK ←	putative permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,117,911	0	A	G	86.0%	478.7 / 60.5	185	V462A (GTA→GCA)	yghK	putative permease
Reads supporting (aligned to +/- strand): ref base A (12/14); new base G (103/56); total (115/70)
Fisher's exact test for biased strand distribution p-value = 8.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

Only 100 of 210 total aligned reads displayed.

GCAGGGTGTCAGAGACGTTGATTTGCTGCGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAA  >  NC_000913/3117819‑3118000

   gggTGTCAGAGACGTTGATTTGCTGCGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATgcgc                                                                                        >  1:3508019/1‑93 (MQ=255)

   gggTGTCAGAGACGTTGATTTGCTGCGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATgcgc                                                                                        >  2:4637402/1‑93 (MQ=255)

   gggTGTCAGAGACGTTGATTTGCTGCGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATgcgc                                                                                        >  1:5335630/1‑93 (MQ=255)

   gggTGTCAGAGACGTTGATTTGCTGCGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGc                                                                                        >  1:48471/1‑93 (MQ=255)

    ggTGTCAGAGACGTTGATTTGCTGCGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGcc                                                                                       >  2:3019421/1‑93 (MQ=255)

     gtgtCAGAGACGTTGATTTGCTGCGCCGTGGTCGATTGCAGTGAACATAACAGGGCGTTAGAGGAGGTTTCCGAGCCGGTAAGGAATGCGccc                                                                                      <  2:4308726/93‑1 (MQ=255)

      tgtCAGAGACGTTGATTTGCTGCGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCcca                                                                                     <  2:4979540/93‑1 (MQ=255)

      tgtCAGAGACGTTGATTTGCTGCGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCcca                                                                                     <  1:1053742/93‑1 (MQ=255)

      tgtCAGAGACGTTGATTTGCTGCGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCcca                                                                                     <  1:3198754/93‑1 (MQ=255)

       gtCAGAGACGTTGATTTGCTGCGCCATGATCGATTGCAGTGAACCAAACGGGGCGTTGGAGGAGGTGTCGGAGCCGGTAAGGAATGCGCccag                                                                                    <  1:2162501/93‑1 (MQ=255)

        tCAGAGACGTTGATTTGCTGCGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCccagc                                                                                   <  2:6369354/93‑1 (MQ=255)

          agagACGTTGATTTGCTGCGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCccagcca                                                                                 <  2:6175286/93‑1 (MQ=255)

          agagACGTTGATTTGCTGCGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCccagcca                                                                                 <  2:5324027/93‑1 (MQ=255)

           gagaCGTTGATTTGCTGCGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCccagccag                                                                                >  1:4581023/1‑93 (MQ=255)

           gagaCGTTGATTTGCTGCGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCccagccag                                                                                >  1:1946769/1‑93 (MQ=255)

           gagaCGTTGATTTGCTGCGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCccagccag                                                                                >  1:6326688/1‑93 (MQ=255)

           gagaCGTTGATTTGCTGCGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCccagccac                                                                                >  1:5109502/1‑92 (MQ=255)

            agaCGTTGATTTGCTGCCCCGTGGTCGATTGCAGTGAGCCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCcagccagc                                                                               <  2:4119925/93‑1 (MQ=255)

              aCGTTGATTTGCTGCGCCGTGGTCGATTGCAGTGAACAAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCg                                                                             <  2:5129812/93‑1 (MQ=255)

               cGTTGATTTGCTGCGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCga                                                                            <  2:2199119/93‑1 (MQ=255)

               cGTTGATTTGCTGCGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCga                                                                            >  1:3271572/1‑93 (MQ=255)

                gTTGATTTGCTGCGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCgag                                                                           >  1:1362143/1‑93 (MQ=255)

                gTTGATTTGCTGCGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCgag                                                                           >  2:6530046/1‑93 (MQ=255)

                 ttGATTTGCTGCGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCgaga                                                                          >  2:806408/1‑93 (MQ=255)

                    aTTTGCGGCCCCGTGGTCGATTGCAGTAAACCAACCAGGGCTTTAGAGGAGGGGCCGGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAAc                                                                       <  2:163795/93‑1 (MQ=255)

                        gctgcGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTg                                                                   >  1:282140/1‑93 (MQ=255)

                        gctgcGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTg                                                                   <  2:4370877/93‑1 (MQ=255)

                          tgcGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTgag                                                                 >  1:4816926/1‑93 (MQ=255)

                          tgcGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTgag                                                                 >  2:6006134/1‑93 (MQ=255)

                           gcgcCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTgaga                                                                >  2:56375/1‑93 (MQ=255)

                            cgcTGTGGTCGATTGCAGTGACCCAAACGGGGCGTGAGAGGGGGTTTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGagaa                                                               <  2:4881890/93‑1 (MQ=25)

                            cgcCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGagaa                                                               <  1:4331234/93‑1 (MQ=255)

                            cgcCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAAAGGTGagaa                                                               <  1:6458721/93‑1 (MQ=255)

                             gcCGTGGTCGATTGCGGTGAACCAAACGGGCCGTTAGAGGAGGTGCCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGagaag                                                              <  2:1066363/93‑1 (MQ=255)

                             gcCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGagaag                                                              >  1:2467079/1‑93 (MQ=255)

                             gcCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGagaag                                                              >  2:1820786/1‑93 (MQ=255)

                              ccGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGagaaga                                                             >  1:5425625/1‑93 (MQ=255)

                              ccGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGagaaga                                                             >  1:4008096/1‑93 (MQ=255)

                               cGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAgaagaa                                                            >  1:4138334/1‑93 (MQ=255)

                               cGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAgaagaa                                                            >  1:6253446/1‑93 (MQ=255)

                                gTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAAc                                                           >  1:823964/1‑93 (MQ=255)

                                gTGGTCGATTGCAGGAACCCAAACGGGCCGTTAGAGGAGGTGCCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAAc                                                           <  2:158886/93‑1 (MQ=255)

                                 tggtggATTGCAGTGACCAAAACAGGGCGTGAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACg                                                          <  1:912833/93‑1 (MQ=255)

                                 tGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACg                                                          >  2:5793108/1‑93 (MQ=255)

                                  ggTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACgg                                                         >  1:1136135/1‑93 (MQ=255)

                                      gATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAAc                                                     >  1:968672/1‑93 (MQ=255)

                                       aTTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACa                                                    >  1:4198091/1‑93 (MQ=255)

                                       aTTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACa                                                    >  1:1672233/1‑93 (MQ=255)

                                       aTTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACa                                                    >  2:2969266/1‑93 (MQ=255)

                                        ttGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACAt                                                   >  1:5193632/1‑93 (MQ=255)

                                        ttGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACAt                                                   >  1:2380302/1‑93 (MQ=255)

                                         tGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATc                                                  >  1:3098814/1‑93 (MQ=255)

                                         tGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATc                                                  >  1:1900505/1‑93 (MQ=255)

                                          gCAGTGAACCAAACAGGGCTTTAGAGGGGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCa                                                 <  2:5524842/93‑1 (MQ=255)

                                          gCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTTCGAGCCGGTAAGGAATGCGCCCAGGCAGCCGAGAAACGGTGAGAAGAACGAGAACATCa                                                 >  2:3287418/1‑93 (MQ=255)

                                          gCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCa                                                 >  1:603504/1‑93 (MQ=255)

                                           cAGTGACCAAAACGGGCCGTTAGAGGAGGTGCCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCAc                                                <  1:6403471/93‑1 (MQ=25)

                                           cAGTGACCAAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCAc                                                <  1:5454606/93‑1 (MQ=21)

                                            aGTGAACAAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACg                                               <  2:2067492/93‑1 (MQ=255)

                                             gTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGc                                              >  2:2630515/1‑93 (MQ=255)

                                              taacccAAAGGGGGCGTGAGGGGGGGTCTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGcc                                             <  2:1451131/89‑1 (MQ=21)

                                              gaacccAAACGGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGcc                                             <  2:3814939/89‑1 (MQ=255)

                                                aaCCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTg                                           >  2:2898088/1‑93 (MQ=255)

                                                aaCCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTg                                           >  2:3443462/1‑93 (MQ=255)

                                                  ccAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTa                                         >  2:124541/1‑93 (MQ=255)

                                                  ccAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGCAAGGAATGCGCCCAGCCAGCCGAGACACGGTGAGAAGAACGAGCACATCACGCCTGTa                                         >  1:1810166/1‑93 (MQ=25)

                                                   cAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGGATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTAc                                        >  2:2540602/1‑93 (MQ=255)

                                                   cAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTAc                                        >  1:4662038/1‑93 (MQ=255)

                                                    aaaaAGGGTGTTAGGGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTAcc                                       <  1:1958748/93‑1 (MQ=25)

                                                     aaCAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCt                                      >  2:87799/1‑93 (MQ=255)

                                                     aaCAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCt                                      >  2:6010800/1‑93 (MQ=255)

                                                        aGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGcc                                   >  2:1872976/1‑93 (MQ=255)

                                                        aGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGcc                                   >  2:84538/1‑93 (MQ=255)

                                                          gggtttAGAGGAGGTGTCCGAGCCGTTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAg                                 <  2:2712080/89‑1 (MQ=37)

                                                             gTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAc                              >  2:3614032/1‑93 (MQ=255)

                                                             gTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAc                              >  2:5038465/1‑93 (MQ=255)

                                                                 gaggagGTGCCGGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAccagc                          <  2:5394554/93‑1 (MQ=255)

                                                                  aggaggTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATGACGCCTGTACCTGCCAGTAccagcg                         >  1:5513416/1‑93 (MQ=255)

                                                                  aggaggTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAccagcg                         >  1:4464902/1‑93 (MQ=255)

                                                                   ggaggTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAccagcgc                        >  1:2752082/1‑93 (MQ=255)

                                                                    gaggTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAccagcgcc                       >  2:6222839/1‑93 (MQ=255)

                                                                    gaggTACCCGACCCGTTAAGAAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAAAACATCACGCCTGTACCTGCCAGTAccagcgcc                       <  1:3803865/93‑1 (MQ=255)

                                                                     aggTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAccagcgcca                      >  2:3359855/1‑93 (MQ=255)

                                                                      ggTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAccagcgccag                     >  1:3610822/1‑93 (MQ=255)

                                                                         gcccGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCgt                  <  2:5235195/91‑1 (MQ=255)

                                                                          gcggAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCgtg                 <  2:4475702/90‑1 (MQ=255)

                                                                            cGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCgtggt               <  1:4498722/93‑1 (MQ=255)

                                                                            cGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCgtggt               <  2:657629/93‑1 (MQ=255)

                                                                            cGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCgtggt               <  1:6149973/93‑1 (MQ=255)

                                                                             gAGCGGGAAAGGAAGGCGCCCAGCCGGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCgtggtg              <  2:2399610/93‑1 (MQ=255)

                                                                             gAGCCGGTAAGGAATCCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTCCCTGCCAGTACCAGCGCCAGCgtggtg              <  2:4213903/93‑1 (MQ=255)

                                                                             gAGCCGGTAAGAAATCCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCgtggtg              <  1:3579091/93‑1 (MQ=255)

                                                                                cgggTAAGGAATGCGCACGGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGAc           <  1:743423/91‑1 (MQ=255)

                                                                                 cGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGACCGAGAACATCACGCTTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACa          <  2:5914295/93‑1 (MQ=255)

                                                                                 cGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCGGCGCCAGCGTGGTGGACa          >  1:5908432/1‑93 (MQ=255)

                                                                                 cGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACa          >  2:2813028/1‑93 (MQ=255)

                                                                                      aGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCa     >  2:3707345/1‑93 (MQ=255)

                                                                                      aGGAAGCCGCCCAGCCAGCCAAGAAACGGTGAGAAGAACGAGACCATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCa     <  2:1503600/93‑1 (MQ=255)

                                                                                       ggAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAg    >  2:4701764/1‑93 (MQ=255)

                                                                                        gAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGa   >  1:5409170/1‑93 (MQ=255)

GCAGGGTGTCAGAGACGTTGATTTGCTGCGCCGTGGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAA  >  NC_000913/3117819‑3118000

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