BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,873,986	T→C	T199A (ACG→GCG)	yidA ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,873,986	0	T	C	99.8%	1411.0 / ‑5.3	408	T199A (ACG→GCG)	yidA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (193/214); total (194/214)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.80e-01

Only 100 of 419 total aligned reads displayed.

CAATCATTGCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATAC  >  NC_000913/3873896‑3874075

cAATCATTGCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCAc                                                                                         >  2:572485/1‑93 (MQ=255)

cAATCATTGCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCAc                                                                                         <  2:5555434/93‑1 (MQ=255)

 aaTCATTGCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCAcc                                                                                        >  2:1927038/1‑93 (MQ=255)

  aTCATTGCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCt                                                                                       <  1:7710267/93‑1 (MQ=255)

  aTCATTGCGATATCGTTTTCCTGATCGCCAATCGCCATCATTTCTTCCGGTTTAATACCTAACACGGCGGCCACTGATTTCACCCCCGCACCt                                                                                       >  1:5524528/1‑93 (MQ=255)

   tCATTGCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCtt                                                                                      >  1:3396640/1‑93 (MQ=255)

    cATTGCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCttt                                                                                     <  1:6687768/93‑1 (MQ=255)

    cATTGCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCttt                                                                                     <  1:286410/93‑1 (MQ=255)

     aTTGCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTg                                                                                    <  2:7128708/93‑1 (MQ=255)

       tGCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGtt                                                                                  <  1:5639775/93‑1 (MQ=255)

        gCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTa                                                                                 >  2:2471534/1‑93 (MQ=255)

         cGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTaa                                                                                <  2:894151/93‑1 (MQ=255)

          gATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAAc                                                                               <  1:4597130/93‑1 (MQ=255)

          gATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAAc                                                                               >  1:5172122/1‑93 (MQ=255)

          gATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAAc                                                                               <  2:4026341/93‑1 (MQ=255)

          gATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAAc                                                                               <  2:418867/93‑1 (MQ=255)

           atatCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAAcg                                                                              >  1:165785/1‑93 (MQ=255)

           atatCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAAcg                                                                              >  1:7510819/1‑93 (MQ=255)

            tatCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAAcgc                                                                             <  2:3090983/93‑1 (MQ=255)

             atCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAAcgcg                                                                            >  2:2821387/1‑93 (MQ=255)

              tCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGt                                                                           >  1:2006335/1‑93 (MQ=255)

               cGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGtt                                                                          <  2:1572004/93‑1 (MQ=255)

                gTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGttt                                                                         >  2:1415857/1‑93 (MQ=255)

                gTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGttt                                                                         >  2:1771173/1‑93 (MQ=255)

                gTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGttt                                                                         <  1:4010528/93‑1 (MQ=255)

                 ttttCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTa                                                                        >  2:5257247/1‑93 (MQ=255)

                  tttCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTAt                                                                       <  1:147213/93‑1 (MQ=255)

                   ttCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACGCGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATc                                                                      <  2:2580479/93‑1 (MQ=255)

                     ccTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGa                                                                    <  2:7032669/93‑1 (MQ=255)

                      cTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAg                                                                   <  2:5090846/93‑1 (MQ=255)

                        gATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGa                                                                 <  2:6552045/93‑1 (MQ=255)

                          tcgccaatcgccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCGCCTTTGTTAACGCGTTTATCGAGGAtt                                                               <  1:3688182/93‑1 (MQ=255)

                          tcgccaatcgccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGAtt                                                               <  1:3084183/93‑1 (MQ=255)

                           cgccaatcgccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGAttt                                                              <  1:2139815/93‑1 (MQ=255)

                           cgccaatcgccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGAttt                                                              <  2:6828894/93‑1 (MQ=255)

                            gccaatcgccaTGATTTCTTCCGGTTTAATGCCTAACACGTCGGCCAGTGATTTCACCCCCGCCCCTTTGTTAACGCGTTTATCGAGGATTTc                                                             <  1:7295077/93‑1 (MQ=255)

                            gccaatcgccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTc                                                             <  1:6232395/93‑1 (MQ=255)

                             ccaatcgccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCGCCTTTGTTAACGCGTTTATCGAGGATTTCg                                                            <  2:4578570/93‑1 (MQ=255)

                              caatcgccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGa                                                           >  1:7780204/1‑93 (MQ=255)

                              caatcgccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGa                                                           <  1:3603173/93‑1 (MQ=255)

                               aatcgccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAg                                                          >  2:7666659/1‑93 (MQ=255)

                               aatcgccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAg                                                          >  2:5418655/1‑93 (MQ=255)

                               aatcgccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGGTTATCGAGGATTTCGAg                                                          >  2:3270319/1‑93 (MQ=255)

                                atcgccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCGCCTTTGTTAACGCGTTTATCGAGGATTTCGAgg                                                         >  1:1323196/1‑93 (MQ=255)

                                atcgccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAga                                                         >  1:832687/1‑92 (MQ=255)

                                  cgccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGaa                                                       <  1:8776794/93‑1 (MQ=255)

                                   gccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAg                                                      >  1:2275035/1‑93 (MQ=255)

                                   gccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAg                                                      <  2:4876003/93‑1 (MQ=255)

                                    ccaTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGt                                                     >  1:4781855/1‑93 (MQ=255)

                                     caTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTa                                                    >  2:1662408/1‑93 (MQ=255)

                                     caTGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTa                                                    >  1:946358/1‑93 (MQ=255)

                                       tGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACg                                                  <  1:1857743/93‑1 (MQ=255)

                                       tGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACg                                                  <  2:4319500/93‑1 (MQ=255)

                                       tGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACg                                                  <  2:4714145/93‑1 (MQ=255)

                                        gATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACgg                                                 >  2:6223702/1‑93 (MQ=255)

                                          tttCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGgcg                                               >  1:5019772/1‑93 (MQ=255)

                                           ttcttcCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGgcgc                                              >  1:107956/1‑93 (MQ=255)

                                            tcttcCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCa                                             <  2:3631667/93‑1 (MQ=255)

                                             cttcCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCAc                                            <  2:3511493/93‑1 (MQ=255)

                                              ttcCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTCCGGCGCACt                                           <  2:5819355/93‑1 (MQ=255)

                                              ttcCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACt                                           >  1:6139534/1‑93 (MQ=255)

                                               tcCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCGCCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACtt                                          <  1:6458372/93‑1 (MQ=255)

                                                   gTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCa                                      <  1:4390499/93‑1 (MQ=255)

                                                   gTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCa                                      <  1:4204317/93‑1 (MQ=255)

                                                   gTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCa                                      <  2:8017261/93‑1 (MQ=255)

                                                    tttAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAg                                     >  1:4709717/1‑93 (MQ=255)

                                                      tAATACCTAACACGTCGGCCAGTGATTTCCCCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCa                                   <  1:2599965/93‑1 (MQ=255)

                                                      tAATACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCa                                   <  1:4072036/93‑1 (MQ=255)

                                                         tACCTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACgg                                <  2:2054568/93‑1 (MQ=255)

                                                            cTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGtat                             <  1:3382558/93‑1 (MQ=255)

                                                            cTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGtat                             >  1:7029858/1‑93 (MQ=255)

                                                            cTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGtat                             <  1:1694785/93‑1 (MQ=255)

                                                            cTAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGtat                             >  2:1080628/1‑93 (MQ=255)

                                                             tAATACGTCGACCAGTGATTTCACCCCAGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGtata                            <  2:5926237/93‑1 (MQ=255)

                                                             tAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGtata                            >  2:241901/1‑93 (MQ=255)

                                                             tAACACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGtata                            <  2:4464357/93‑1 (MQ=255)

                                                              aaCACGTCGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTatat                           >  2:3772328/1‑93 (MQ=255)

                                                              aaCACGTCGGCCAGTGATTTAACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGAACGGCGCAATTTTAAGCACGGGatat                           >  2:7419448/1‑93 (MQ=255)

                                                              aaCACGTAGGCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTatat                           >  2:8643738/1‑93 (MQ=255)

                                                                    tCGGCCAGTGATTTCCCCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTctct                     <  2:1189878/93‑1 (MQ=255)

                                                                       gCCAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTc                  <  1:4752353/93‑1 (MQ=255)

                                                                        ccAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCa                 >  1:8511211/1‑93 (MQ=255)

                                                                        ccAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCa                 <  2:8390361/93‑1 (MQ=255)

                                                                         cAGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCAc                <  1:3796143/93‑1 (MQ=255)

                                                                          aGTGCTTTCTCCCCCGCACCTTTGTTAACGCGTTTATCGAGGTTTTCGAGGAAGTGCGGCACCCTTTTCCGCACGGTATATTTCTCTTTCACt               <  2:6675569/93‑1 (MQ=255)

                                                                          aGTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACt               >  2:8210660/1‑93 (MQ=255)

                                                                           gTGATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACtt              <  2:2713827/93‑1 (MQ=255)

                                                                            ggATTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTc             <  2:4887065/92‑1 (MQ=255)

                                                                              aTTTCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCt           <  2:7283088/93‑1 (MQ=255)

                                                                               tttCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTg          <  1:2456802/93‑1 (MQ=255)

                                                                                ttCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGc         <  2:5622034/93‑1 (MQ=255)

                                                                                 tCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCg        <  2:3913911/93‑1 (MQ=255)

                                                                                 tCACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCg        <  2:2327114/93‑1 (MQ=255)

                                                                                  cACCCCCGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCgg       <  1:5279292/93‑1 (MQ=255)

                                                                                    cccccGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGaa     <  2:6078450/93‑1 (MQ=255)

                                                                                      cccGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATa   <  2:3196920/93‑1 (MQ=255)

                                                                                       ccGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCCCGGCATATTTCTCTTTCACTTCCTGCGGAATAc  >  2:7579662/1‑93 (MQ=255)

                                                                                       ccGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATAc  <  2:2281690/93‑1 (MQ=255)

                                                                                       ccGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATAc  >  1:5964041/1‑93 (MQ=255)

                                                                                       ccGCACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATAc  >  2:1008136/1‑93 (MQ=255)

CAATCATTGCGATATCGTTTTCCTGATCGCCAATCGCCATGATTTCTTCCGGTTTAATACCTAACACGTCGGCCAGTGATTTCACCCCCGTACCTTTGTTAACGCGTTTATCGAGGATTTCGAGGAAGTACGGCGCACTTTTCAGCACGGTATATTTCTCTTTCACTTCCTGCGGAATAC  >  NC_000913/3873896‑3874075

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