BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,455,292	T→C	W427R (TGG→CGG)	b3325 →	YheF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,455,292	0	T	C	99.6%	779.0 / ‑3.7	225	coding (1279/1953 nt)	b3325	YheF
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (119/105); minor base . (0/1); total (119/106)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.28e-01

Only 100 of 233 total aligned reads displayed.

AAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTCGCG  >  NC_000913/3455202‑3455381

aaaaaGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCgg                                                                                         <  1:4266072/93‑1 (MQ=255)

  aaaGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCgggg                                                                                       <  1:4048868/93‑1 (MQ=255)

  aaaGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCgggg                                                                                       <  2:4101647/93‑1 (MQ=255)

   aaGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGa                                                                                      <  2:427687/93‑1 (MQ=255)

     gAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGGGt                                                                                    >  1:6034416/1‑93 (MQ=255)

     gAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAga                                                                                    >  2:2196427/1‑92 (MQ=255)

      aaTGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTa                                                                                   <  1:324345/93‑1 (MQ=255)

       aTGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTAc                                                                                  >  1:3691584/1‑93 (MQ=255)

       aTGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTAc                                                                                  <  2:1135/93‑1 (MQ=255)

       aTGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTAc                                                                                  >  2:1236635/1‑93 (MQ=255)

       aTGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTAc                                                                                  <  1:637052/93‑1 (MQ=255)

       aTGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTAc                                                                                  >  1:1245176/1‑93 (MQ=255)

       aTGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTAc                                                                                  <  1:2668454/93‑1 (MQ=255)

       aTGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGCAc                                                                                  >  1:864200/1‑93 (MQ=255)

        tggtggGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACt                                                                                 >  2:5139430/1‑93 (MQ=255)

           tggGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCt                                                                              <  1:4902132/93‑1 (MQ=255)

            gggATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCATACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCtt                                                                             <  1:4806706/93‑1 (MQ=255)

              gATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTAc                                                                           >  1:5516311/1‑93 (MQ=255)

              gATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTAc                                                                           <  1:992332/93‑1 (MQ=255)

               aTCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTatc                                                                          <  2:1369879/93‑3 (MQ=255)

               aTCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTAcc                                                                          >  1:408509/1‑93 (MQ=255)

               aTCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTAcc                                                                          >  2:826953/1‑93 (MQ=255)

               aTCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTAcc                                                                          >  2:2555106/1‑93 (MQ=255)

                 cACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACcgc                                                                        <  1:2149586/93‑1 (MQ=255)

                   ccAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCgcgc                                                                      >  1:3910542/1‑93 (MQ=255)

                   ccAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCgcgc                                                                      <  1:4111901/93‑1 (MQ=255)

                     aGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTg                                                                    <  1:1788009/93‑1 (MQ=255)

                      gcgcGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTgg                                                                   >  2:5353674/1‑93 (MQ=255)

                        gcgAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGcc                                                                 >  1:6027302/1‑93 (MQ=255)

                        gcgAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGcc                                                                 <  1:5234762/93‑1 (MQ=255)

                        gcgAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGcc                                                                 <  2:3342895/93‑1 (MQ=255)

                        gcgAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGcc                                                                 <  2:2919478/93‑1 (MQ=255)

                         cgAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCa                                                                >  1:466345/1‑93 (MQ=255)

                         cgAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCa                                                                <  2:2007607/93‑1 (MQ=255)

                         cgAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCa                                                                <  2:5650668/93‑1 (MQ=255)

                           aaTCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGt                                                              >  1:16664/1‑93 (MQ=255)

                           aaTCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGt                                                              >  2:619491/1‑93 (MQ=255)

                            aTCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTa                                                             <  1:4189472/93‑1 (MQ=255)

                            aTCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTa                                                             >  2:2171172/1‑93 (MQ=255)

                              cctgcctgGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAAc                                                           <  2:47181/93‑1 (MQ=255)

                              cctgcctgGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAAc                                                           >  1:3081873/1‑93 (MQ=255)

                                 gcctgGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAAt                                                        <  2:3054033/93‑1 (MQ=255)

                                 gcctgGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAAt                                                        <  2:211268/93‑1 (MQ=255)

                                 gcctgGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAAt                                                        <  1:914828/93‑1 (MQ=255)

                                 gcctgGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAAt                                                        <  1:1288496/93‑1 (MQ=255)

                                  cgngngnTATGTTTAGCCCCCNNAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATa                                                       <  2:2489267/90‑1 (MQ=255)

                                  cctgGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATa                                                       >  1:3028617/1‑93 (MQ=255)

                                  cctgGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATa                                                       >  1:2532932/1‑93 (MQ=255)

                                  cctgGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATa                                                       >  1:5338459/1‑93 (MQ=255)

                                   ctgGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATaa                                                      >  2:5030400/1‑93 (MQ=255)

                                     gggATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCACTGGCGACCGTGGAGTACTGCTTAACGCGCTGGCCAGTAACAATaaaa                                                    >  2:3918814/1‑93 (MQ=255)

                                     gggATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATaaaa                                                    <  2:4380907/93‑1 (MQ=255)

                                     gggATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATaaaa                                                    >  1:3948008/1‑93 (MQ=255)

                                     gggATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATaaaa                                                    >  2:6322849/1‑93 (MQ=255)

                                      ggATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATaaaaa                                                   >  1:329356/1‑93 (MQ=255)

                                       gATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAAt                                                  <  2:5940697/93‑1 (MQ=255)

                                       gATATGTGTAGCGCCTACAATGGCATGGCCGCTGGCTTCTTCAATGGCGACCGGGGAGTCCTGCTTACCGCGCTGGCCAGTAACAATAAAAAt                                                  <  2:6409099/93‑1 (MQ=255)

                                        atatGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATg                                                 >  1:5272091/1‑93 (MQ=255)

                                        atatGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATg                                                 <  2:1413115/93‑1 (MQ=255)

                                         tatGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGa                                                >  1:5046281/1‑93 (MQ=255)

                                           tGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACa                                              >  2:6115513/1‑93 (MQ=255)

                                            gTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACAt                                             <  1:6255219/93‑1 (MQ=255)

                                            gTTGAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACAt                                             <  2:2523483/93‑1 (MQ=255)

                                               tAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTCCTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCt                                          <  1:3161637/93‑1 (MQ=255)

                                                aGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTc                                         >  2:675282/1‑93 (MQ=255)

                                                   gcCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGcc                                      >  1:2546538/1‑93 (MQ=255)

                                                    ccTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCa                                     >  2:3472617/1‑93 (MQ=255)

                                                     cTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCAc                                    >  1:5892871/1‑93 (MQ=255)

                                                      tACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCAcc                                   >  2:5473875/1‑93 (MQ=255)

                                                       aCAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCAccc                                  <  2:1801600/93‑1 (MQ=255)

                                                       aCAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCAccc                                  >  2:1120654/1‑93 (MQ=255)

                                                       aCAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCAccc                                  <  2:3766024/93‑1 (MQ=255)

                                                        cAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCAcccc                                 <  1:1477586/93‑1 (MQ=255)

                                                         aaTGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCa                                <  2:2805285/93‑1 (MQ=255)

                                                         aaTGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCa                                <  2:3995827/93‑1 (MQ=255)

                                                             gcatggcCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCa                            <  1:542200/93‑1 (MQ=255)

                                                               atggcCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATc                          <  1:1615810/93‑1 (MQ=255)

                                                                tggcCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCg                         <  2:5710269/93‑1 (MQ=255)

                                                                 ggcCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGt                        >  1:1022548/1‑93 (MQ=255)

                                                                    cGCAGGCTTCTTCAACGGCGACGGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAAc                     <  1:1169463/93‑1 (MQ=255)

                                                                     gCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACg                    >  2:1452001/1‑93 (MQ=255)

                                                                      cAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGc                   <  1:5651056/93‑1 (MQ=255)

                                                                      cAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGc                   <  1:806092/93‑1 (MQ=255)

                                                                       aGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCt                  >  1:2298120/1‑93 (MQ=255)

                                                                           ttcttcAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGat              <  2:4851183/93‑1 (MQ=255)

                                                                           ttcttcAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGat              <  2:706139/93‑1 (MQ=255)

                                                                           ttcttcAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGat              <  2:1832032/93‑1 (MQ=255)

                                                                           ttcttcAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGat              >  2:5972757/1‑93 (MQ=255)

                                                                           ttcttcAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGat              >  1:1954004/1‑93 (MQ=255)

                                                                            tcttcAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGata             <  2:2636921/93‑1 (MQ=255)

                                                                             cttcAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGataa            >  2:584847/1‑93 (MQ=255)

                                                                              ttcAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGataat           >  1:5775911/1‑93 (MQ=255)

                                                                                 aaTGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATaaa        >  1:2275961/1‑93 (MQ=255)

                                                                                  aTGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAAc       <  1:1375915/93‑1 (MQ=255)

                                                                                  aTGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAAc       >  2:3808491/1‑93 (MQ=255)

                                                                                    ggCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTc     <  1:3349365/93‑1 (MQ=255)

                                                                                      cGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTcgc   <  1:4223306/93‑1 (MQ=255)

                                                                                       gACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTcgcg  <  2:6049036/93‑1 (MQ=255)

                                                                                       gACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTcgcg  >  1:4429975/1‑93 (MQ=255)

                                                                                       gACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTcgcg  >  2:806405/1‑93 (MQ=255)

AAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTCGCG  >  NC_000913/3455202‑3455381

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