BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,016,822	G→A	L127L (CTG→CTA)	ubiE →	UbiE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,016,822	0	G	A	99.5%	1292.4 / ‑4.7	371	L127L (CTG→CTA)	ubiE	UbiE
Reads supporting (aligned to +/- strand): ref base G (1/1); new base A (165/203); total (166/205)
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.40e-01

Only 100 of 383 total aligned reads displayed.

TGAATCCATGCCCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGCGTTCAATGTATC  >  NC_000913/4016732‑4016913

  aaTCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGt                                                                                         >  2:4321308/1‑93 (MQ=255)

   aTCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGtt                                                                                        >  2:4057617/1‑93 (MQ=255)

    tCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTc                                                                                       >  1:7461036/1‑93 (MQ=255)

    tCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTc                                                                                       >  1:6661635/1‑93 (MQ=255)

      cATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTccc                                                                                     <  1:7058598/93‑1 (MQ=255)

       aTGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCg                                                                                    >  1:608138/1‑93 (MQ=255)

       aTGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCg                                                                                    >  1:4670984/1‑93 (MQ=255)

       aTGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGAGCTACCGTTCCCg                                                                                    >  2:4190539/1‑93 (MQ=255)

        tGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCgg                                                                                   >  1:7188068/1‑93 (MQ=255)

        tGCTCAAAATGGGCCGCGAGAAGCTGCGGAATATCGGTGTGATTGGCAACGTTGAGTGGGTTCAGGCGTACGCGGAGGCGCTACCGGTCCCgg                                                                                   >  1:7020490/1‑93 (MQ=255)

         gctcAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGa                                                                                  <  1:1301058/90‑1 (MQ=255)

           tcAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATa                                                                                <  1:2167193/92‑1 (MQ=255)

           tcAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATa                                                                                >  2:4151261/2‑93 (MQ=255)

            cAAAATGGGAGGCGAGAAGCTGCGTAATATCGGTGGCATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTTAGGCGCTACCGTTCCCTGATaa                                                                               <  1:7383321/93‑1 (MQ=255)

             aaaaTGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAac                                                                              >  2:6494115/1‑93 (MQ=255)

             aaaaTGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAac                                                                              <  1:7916145/93‑1 (MQ=255)

             aaaaTGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAac                                                                              >  1:3345648/1‑93 (MQ=255)

               aaTGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAacac                                                                            >  1:2726659/1‑93 (MQ=255)

               aaTGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAacac                                                                            >  1:6436624/1‑93 (MQ=255)

                aggggCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACAcc                                                                           <  2:109622/91‑1 (MQ=255)

                aTGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACAcc                                                                           <  1:5838666/93‑1 (MQ=255)

                aTGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACAcc                                                                           >  2:1213858/1‑93 (MQ=255)

                 tGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCt                                                                          >  2:6315312/1‑93 (MQ=255)

                 tGGGACGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCt                                                                          <  1:3134606/93‑1 (MQ=255)

                   ggCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCttt                                                                        <  2:1700925/93‑1 (MQ=255)

                   ggCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCttt                                                                        <  1:1014656/93‑1 (MQ=255)

                     ccGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGa                                                                      <  1:1701131/93‑1 (MQ=255)

                      cgcgAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGAt                                                                     <  2:6931998/93‑1 (MQ=255)

                       gcgAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGAtt                                                                    <  2:7164843/93‑1 (MQ=255)

                       gcgAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGAtt                                                                    >  2:2637000/1‑93 (MQ=255)

                          agaAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCa                                                                 <  2:3606703/93‑1 (MQ=255)

                           gaAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCAt                                                                >  1:5258563/1‑93 (MQ=255)

                             aGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCa                                                              <  1:4898007/93‑1 (MQ=255)

                              gCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCAc                                                             <  2:781804/93‑1 (MQ=255)

                                tGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCa                                                           <  2:7286587/93‑1 (MQ=255)

                                 gagTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCAt                                                          <  1:5139883/91‑1 (MQ=255)

                                 gCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCAt                                                          >  1:342942/1‑93 (MQ=255)

                                  cGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCAtt                                                         <  2:5741185/93‑1 (MQ=255)

                                  cGTAATACCGGTGTGATTGGCAACGTTGAGTGTGTTCAGCCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCAtt                                                         <  1:5023261/93‑1 (MQ=255)

                                   gTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCAttt                                                        >  2:4663632/1‑93 (MQ=255)

                                   gTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCAttt                                                        <  2:1308418/93‑1 (MQ=255)

                                    tAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTc                                                       <  2:7412581/93‑1 (MQ=255)

                                    tAATATCGGTGTGATTGGCAACGTTGAGTGGGTTCAGGCGAACGCTGAGGCGCTAACGGTCCCTGATAACACCTTTGATTGCATCACCATTTc                                                       >  2:2104528/1‑93 (MQ=255)

                                     aaTATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCg                                                      >  2:7710761/1‑93 (MQ=255)

                                      atatCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGt                                                     <  2:5255715/93‑1 (MQ=255)

                                       tatCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGtt                                                    >  2:3430234/1‑93 (MQ=255)

                                       tatCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGtt                                                    >  2:561578/1‑93 (MQ=255)

                                       tatCGGTGTGATTGGCAACGTTGAATGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGtt                                                    >  1:5446940/1‑93 (MQ=255)

                                        atCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGttt                                                   <  1:1218839/93‑1 (MQ=255)

                                        atCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGttt                                                   <  1:3794666/93‑1 (MQ=255)

                                        atCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGttt                                                   <  2:2328636/93‑1 (MQ=255)

                                          cGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTgg                                                 <  2:6984801/93‑1 (MQ=255)

                                          cGGTGTGATTGGCAACGTTGAGTGTGTTAAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTgg                                                 >  2:2240274/1‑93 (MQ=255)

                                            gtgtGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTc                                               >  1:6743861/1‑93 (MQ=255)

                                            gtgtGATTGGCAACGTTGAGTGTGTTCAGGAGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTc                                               >  1:2379111/1‑93 (MQ=255)

                                             tgtgATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCt                                              <  2:684393/93‑1 (MQ=255)

                                               tgATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGc                                            <  1:5743920/93‑1 (MQ=255)

                                               tgATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGc                                            <  2:4433397/93‑1 (MQ=255)

                                               tgATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGc                                            >  2:5464622/1‑93 (MQ=255)

                                                gATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCg                                           >  1:2873110/1‑93 (MQ=255)

                                                 aTTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGATGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGt                                          >  1:441308/1‑93 (MQ=255)

                                                 aTTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGt                                          <  2:3366073/93‑1 (MQ=255)

                                                    ggCAACTTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAAc                                       <  1:4547934/93‑1 (MQ=255)

                                                    ggCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAAc                                       <  1:750283/93‑1 (MQ=255)

                                                       aaCGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTc                                    <  2:82410/93‑1 (MQ=255)

                                                       aaCGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTc                                    >  1:6518360/1‑93 (MQ=255)

                                                         cGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCAc                                  <  2:3170611/93‑1 (MQ=255)

                                                         cGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCAc                                  <  2:6414896/93‑1 (MQ=255)

                                                          gTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCAcc                                 <  2:3681288/93‑1 (MQ=255)

                                                          gTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCAcc                                 >  2:4780580/1‑93 (MQ=255)

                                                           ttGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCg                                >  2:5896729/1‑93 (MQ=255)

                                                            tGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACATCACCGa                               <  2:2666870/93‑1 (MQ=255)

                                                                 gtGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAg                          <  1:6393063/92‑1 (MQ=255)

                                                                 gtGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAg                          >  2:4557147/2‑93 (MQ=255)

                                                                  tGTTCATGCGAACGCTGAGGCGCTACCGTTCACGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGAGAAAGa                         <  2:7932688/93‑1 (MQ=255)

                                                                  tGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGa                         >  1:1057422/1‑93 (MQ=255)

                                                                   gTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGAt                        >  2:4915460/1‑93 (MQ=255)

                                                                   gTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGAt                        >  2:3070905/1‑93 (MQ=255)

                                                                   gTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGAt                        >  2:6739010/1‑93 (MQ=255)

                                                                   gTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGAt                        >  2:617063/1‑93 (MQ=255)

                                                                    ttCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATa                       <  2:1613050/93‑1 (MQ=255)

                                                                     tCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATaa                      >  1:1139379/1‑93 (MQ=255)

                                                                      cAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATaaa                     <  2:3846988/93‑1 (MQ=255)

                                                                      cAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATaaa                     <  2:277708/93‑1 (MQ=255)

                                                                       aGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATTACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAg                    <  1:1903615/93‑1 (MQ=255)

                                                                        ggCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGc                   >  2:2434260/1‑93 (MQ=255)

                                                                        ggCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGc                   <  1:4122034/93‑1 (MQ=255)

                                                                            aaCGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTg               <  1:143324/93‑1 (MQ=255)

                                                                            aaCGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTg               >  1:7660131/1‑93 (MQ=255)

                                                                             aCGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGc              <  2:2848706/93‑1 (MQ=255)

                                                                             aCGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGc              <  1:181728/93‑1 (MQ=255)

                                                                             aCGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGc              <  2:71189/93‑1 (MQ=255)

                                                                               gCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGCGt            >  2:7440817/1‑93 (MQ=255)

                                                                                 tGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGCGTTc          <  2:1781448/93‑1 (MQ=255)

                                                                                  gAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGCGTTCa         >  2:5522091/1‑93 (MQ=255)

                                                                                   aGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGCGTTCaa        >  1:3631321/1‑93 (MQ=255)

                                                                                   aGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGCGTTCaa        <  1:6883966/93‑1 (MQ=255)

                                                                                     gcgcTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGCGTTCAATg      <  1:97702/93‑1 (MQ=255)

                                                                                      cgcTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGCGTTCAATGt     <  1:561522/93‑1 (MQ=255)

                                                                                       gcTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGCGTTCAATGTa    <  1:3236312/93‑1 (MQ=255)

TGAATCCATGCCCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGCGTTCAATGTATC  >  NC_000913/4016732‑4016913

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