BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,500,309	Δ1 bp	coding (333/363 nt)	b3372 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,500,308	0	C	.	99.6%	1073.9 / ‑3.4	230	coding (332/363 nt)	b3372	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (1/0); new base . (133/96); total (134/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 251 total aligned reads displayed.

CGATGAACATATGCGTCGCGAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCGCCGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAACTGTGTGAGTACGCGGAAAA  >  NC_000913/3500216‑3500396

ggATGAACATATGCGTCGCGAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTccgc                                                                                          <  2:6105745/92‑1 (MQ=255)

 gATGAACATATGCGTCGCGAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTccgcg                                                                                         >  1:1297023/1‑92 (MQ=255)

 gATGAACATATGCGTCGCGAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTccgcg                                                                                         <  1:5086596/93‑2 (MQ=255)

 gATGAACATATGCGTCGCGAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTccgcg                                                                                         >  2:5316666/1‑92 (MQ=255)

  aTGCACATATGCGTCGCGAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTccgcgg                                                                                        >  1:5330552/1‑91 (MQ=255)

  aTGAACATATGCGTCGCGAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTccgcgg                                                                                        >  1:6256135/1‑91 (MQ=255)

    gAACATATGCGTCGCGAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTccgcggcc                                                                                      >  2:4071722/1‑89 (MQ=255)

     aaCATATGCGTCGCGAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGccc                                                                                    >  1:1373491/1‑93 (MQ=255)

     aaCATATGCGTCGCGAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGccc                                                                                    <  2:3368959/93‑1 (MQ=255)

      aCATATGCGTCGCGAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCGCGGGcca                                                                                    >  1:1647973/1‑92 (MQ=255)

        atatGCGTCGCGAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCAcg                                                                                 >  1:6118986/1‑93 (MQ=255)

        atatGCGTCGCGAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCAcg                                                                                 >  2:1422680/1‑93 (MQ=255)

         tatGCGTCGCGAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCAcgc                                                                                <  2:2117427/93‑1 (MQ=255)

         tatGCGTCGCGAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCAcgc                                                                                >  1:388514/1‑93 (MQ=255)

           tGCGTCGCGAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCgg                                                                              >  2:680228/1‑93 (MQ=255)

            gCGTCGCGAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCggg                                                                             >  2:3607976/1‑93 (MQ=255)

              gTCGCGAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCt                                                                           <  1:3017544/93‑1 (MQ=255)

               tCGCGAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTa                                                                          <  2:2327382/93‑1 (MQ=255)

                cgcgAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTAt                                                                         <  2:1638066/93‑1 (MQ=255)

                 gcgAAAGCCTCGACATGATCAAGCTGGCGATGTATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTAtc                                                                        <  1:2083715/93‑1 (MQ=255)

                 gcgAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTAtc                                                                        <  1:908709/93‑1 (MQ=255)

                  cgAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTAtct                                                                       <  1:1392316/93‑1 (MQ=255)

                   gAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTAtctc                                                                      <  2:3561168/93‑1 (MQ=255)

                   gAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTAtctc                                                                      <  2:4907874/93‑1 (MQ=255)

                    aaaGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCa                                                                     >  1:3459636/1‑93 (MQ=255)

                     aaGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCAc                                                                    <  2:220866/93‑1 (MQ=255)

                      aGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACg                                                                   >  2:3108936/1‑93 (MQ=255)

                       gCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGc                                                                  <  1:2771077/93‑1 (MQ=255)

                         cTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCa                                                                >  1:6181582/1‑93 (MQ=255)

                          tCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCAc                                                               >  2:2644863/1‑93 (MQ=255)

                           cGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTTACGCCAcc                                                              <  1:6298506/93‑1 (MQ=255)

                           cGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCAcc                                                              <  1:4990091/93‑1 (MQ=255)

                            gACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCt                                                             >  2:1731623/1‑93 (MQ=255)

                            gACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCt                                                             <  1:983920/93‑1 (MQ=255)

                            gACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCt                                                             >  1:2331934/1‑93 (MQ=255)

                              cATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTaa                                                           <  2:5585459/93‑1 (MQ=255)

                                 gATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGt                                                        <  1:4631312/93‑1 (MQ=255)

                                 gATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGt                                                        <  2:5468946/93‑1 (MQ=255)

                                  aTCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGtt                                                       >  2:1588216/1‑93 (MQ=255)

                                     aaGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATc                                                    >  2:2404929/1‑93 (MQ=255)

                                     aaGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATc                                                    <  2:4775510/93‑1 (MQ=255)

                                     aaGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATc                                                    >  2:5993440/1‑93 (MQ=255)

                                      aGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACTCGGGCTATCTCACGCCACCTAATGTTATCt                                                   >  1:3512376/1‑93 (MQ=255)

                                      aGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCt                                                   >  1:5597825/1‑93 (MQ=255)

                                       gCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTg                                                  >  1:668236/1‑93 (MQ=255)

                                       gCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTg                                                  >  2:5085414/1‑93 (MQ=255)

                                        cTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTgg                                                 <  1:6041475/93‑1 (MQ=255)

                                         tGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTggg                                                <  2:3889397/93‑1 (MQ=255)

                                             gATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACg                                            >  2:1203962/1‑93 (MQ=255)

                                              aTGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACgg                                           <  2:2553152/93‑1 (MQ=255)

                                              aTGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACgg                                           <  2:4182984/93‑1 (MQ=255)

                                               tGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGc                                          <  1:4385773/93‑1 (MQ=255)

                                               tGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGc                                          >  2:4873317/1‑93 (MQ=255)

                                                ggATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCt                                         >  1:5527210/1‑93 (MQ=255)

                                                ggATATGGAAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCt                                         >  1:3299375/1‑93 (MQ=255)

                                                 gATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTg                                        >  2:5683095/1‑93 (MQ=255)

                                                 gATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTg                                        <  2:3305759/93‑1 (MQ=255)

                                                 gATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTg                                        >  1:2605970/1‑93 (MQ=255)

                                                   tatGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGc                                      >  1:6141998/1‑93 (MQ=255)

                                                     tGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCag                                    >  1:1952823/1‑93 (MQ=255)

                                                      ggCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCaga                                   <  1:4022584/93‑1 (MQ=255)

                                                      ggCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCaga                                   >  1:3959969/1‑93 (MQ=255)

                                                         aaaaGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGGTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGaa                                >  1:84716/1‑93 (MQ=255)

                                                         aaaaGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGaa                                >  1:5646551/1‑93 (MQ=255)

                                                          aaaGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAAc                               <  2:4166092/93‑1 (MQ=255)

                                                          aaaGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAAc                               <  2:6155267/93‑1 (MQ=255)

                                                           aaGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCGCCTAATGTTATCTGGGGACGGCTGGCAGAGAAcc                              <  2:1828638/93‑1 (MQ=255)

                                                           aaGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAAct                              >  2:890146/1‑92 (MQ=255)

                                                           aaGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAAcc                              >  1:4776896/1‑93 (MQ=255)

                                                           aaGAGATGAACGCGGGTTATACGCTGATTTCAG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGCGGAcc                              >  1:290737/1‑93 (MQ=255)

                                                            agagATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCt                             >  2:5268359/1‑93 (MQ=255)

                                                            agagATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCt                             >  1:308592/1‑93 (MQ=255)

                                                              agaTGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTaa                           <  2:2936125/93‑1 (MQ=255)

                                                              agaTGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTaa                           <  1:4822524/93‑1 (MQ=255)

                                                              agaTGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTaa                           <  1:3656197/93‑1 (MQ=255)

                                                               gaTGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAg                          <  2:3681425/93‑1 (MQ=255)

                                                                aTGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGc                         >  2:3380641/1‑93 (MQ=255)

                                                                   aaCGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGaa                      >  2:447941/1‑93 (MQ=255)

                                                                    aCGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAAc                     <  1:4730552/93‑1 (MQ=255)

                                                                    aCGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAAc                     >  1:3738518/1‑93 (MQ=255)

                                                                    aCGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCAGGCTATCTCACCCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGCAc                     >  1:758424/1‑93 (MQ=255)

                                                                     cgcgGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAACt                    <  2:640609/93‑1 (MQ=255)

                                                                     cgcgGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAACt                    <  2:1472592/93‑1 (MQ=255)

                                                                        gggTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAACtgtg                 <  2:4288124/93‑1 (MQ=255)

                                                                         ggTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAACtgtgt                <  1:4358770/93‑1 (MQ=255)

                                                                         ggTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAACtgtgt                <  2:1691796/93‑1 (MQ=255)

                                                                          gTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAACtgtgtg               <  2:4112883/93‑1 (MQ=255)

                                                                          gTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAACtgtgtg               <  2:1947617/93‑1 (MQ=255)

                                                                           ttATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAACTGTGTGa              <  1:2169642/93‑1 (MQ=255)

                                                                            tataCGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAACTGTGTGAg             >  1:415791/1‑93 (MQ=255)

                                                                             ataCGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAACTGTGTGAGt            <  1:3229759/93‑1 (MQ=255)

                                                                             ataCGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAACTGTGTGAGt            <  1:1031684/93‑1 (MQ=255)

                                                                              taCGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAACTGTGTGAGTa           <  1:5710913/93‑1 (MQ=255)

                                                                              taCGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAACTGTGTGAGTa           <  2:3295403/93‑1 (MQ=255)

                                                                                cGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAACTGTGTGAGTAcg         >  2:5915180/1‑93 (MQ=255)

                                                                                cGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAACTGTGTGAGTAcg         <  1:6378090/93‑1 (MQ=255)

                                                                                   tGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAACTGTGTGAGTACGCgg      >  1:5184091/1‑93 (MQ=255)

                                                                                     aTTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAACTGTGTGAGTACGCGGaa    >  2:5810896/1‑93 (MQ=255)

                                                                                      tttCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGCTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAACTGTGCGAGTACGCGGaaa   >  1:6020185/1‑93 (MQ=255)

                                                                                       ttCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAACTGTGTGAGTACGCGGaaaa  <  2:1285141/93‑1 (MQ=255)

CGATGAACATATGCGTCGCGAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCGCCGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAACTGTGTGAGTACGCGGAAAA  >  NC_000913/3500216‑3500396

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