BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	464,354	C→A	15.0%	intergenic (+46/‑48)	ybaV → / → fadM	putative competence‑suppressing periplasmic helix‑hairpin‑helix DNA‑binding protein/long‑chain acyl‑CoA thioesterase III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	464,354	0	C	A	15.0%	453.4 / 60.9	165	intergenic (+46/‑48)	ybaV/fadM	putative competence‑suppressing periplasmic helix‑hairpin‑helix DNA‑binding protein/long‑chain acyl‑CoA thioesterase III
Reads supporting (aligned to +/- strand): ref base C (63/77); new base A (10/15); total (73/92)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

Only 100 of 174 total aligned reads displayed.

CCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACG  >  NC_000913/464255‑464453

ccGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAc                                                                                                     <  2:3113907/100‑1 (MQ=255)

 cGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAcc                                                                                                    >  1:2633523/1‑100 (MQ=255)

   gggATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAg                                                                                                  <  2:3139873/100‑1 (MQ=255)

         ggCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGa                                                                                            >  2:193676/1‑100 (MQ=255)

           cAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGAcc                                                                                          <  2:796077/100‑1 (MQ=255)

             aTTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACctc                                                                                        <  1:476834/100‑1 (MQ=255)

             aTTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctc                                                                                        <  1:3600064/100‑1 (MQ=255)

              ttCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctct                                                                                       <  2:3647614/100‑1 (MQ=255)

              ttCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctct                                                                                       <  1:1211387/100‑1 (MQ=255)

              ttCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctct                                                                                       <  2:2791418/100‑1 (MQ=255)

              ttCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctct                                                                                       <  1:3133818/100‑1 (MQ=255)

               tCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTg                                                                                      <  2:395798/100‑1 (MQ=255)

               tCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctctt                                                                                      <  1:1234577/100‑2 (MQ=39)

                cGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGCCTGAAGAGGTCATACCAGTTATGACCTCTGt                                                                                     >  1:3527824/1‑100 (MQ=255)

                cGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGt                                                                                     >  1:2067266/1‑100 (MQ=255)

                cGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGt                                                                                     >  1:2522906/1‑100 (MQ=255)

                cGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGt                                                                                     >  2:1863964/1‑100 (MQ=255)

                cGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGt                                                                                     <  2:3507502/100‑1 (MQ=255)

                cGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctcttc                                                                                     <  1:3193550/100‑3 (MQ=38)

                 gCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTa                                                                                    <  1:1133156/100‑1 (MQ=255)

                 gCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctcttca                                                                                    <  2:3070652/100‑4 (MQ=38)

                  cTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGTATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCTGGTATGACctcttcag                                                                                   <  1:874557/100‑5 (MQ=39)

                  cTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTTCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTAc                                                                                   >  2:3463789/1‑100 (MQ=255)

                  cTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTAc                                                                                   <  1:662890/100‑1 (MQ=255)

                  cTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctcttcag                                                                                   <  2:1020894/100‑5 (MQ=38)

                  cTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctcttcag                                                                                   <  1:3215526/100‑5 (MQ=38)

                    ggtggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACtt                                                                                 >  1:649878/1‑100 (MQ=255)

                       ggAACGTAATCTGGCGGTATTGACCCTGTAATTAATTTGTATAGTGGCTATTTTTTGCCAGACTGAAGAGGTAATCCCAGTTATGACCTGTGTACTtata                                                                              <  2:3064425/100‑1 (MQ=255)

                       ggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTtata                                                                              <  1:194478/100‑1 (MQ=255)

                       ggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTtata                                                                              <  1:2686723/100‑1 (MQ=255)

                         aaCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaac                                                                            <  1:6197/100‑1 (MQ=255)

                         aaCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaa                                                                             <  1:909861/99‑1 (MQ=255)

                          aCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaaca                                                                           <  1:1084223/100‑1 (MQ=255)

                          aCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaaca                                                                           >  1:3721636/1‑100 (MQ=255)

                          aCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaaca                                                                           <  1:1402059/100‑1 (MQ=255)

                          aCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaaca                                                                           >  1:2796837/1‑100 (MQ=255)

                           cGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaagat                                                                          >  1:605978/1‑97 (MQ=255)

                             tAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaacaaca                                                                        >  1:3448416/1‑100 (MQ=255)

                                 cTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGt                                                                    >  1:530321/1‑100 (MQ=255)

                                  tGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTa                                                                   <  2:1621169/100‑1 (MQ=255)

                                   ggCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTaa                                                                  >  2:1848857/1‑100 (MQ=255)

                                    gCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAg                                                                 <  1:1848867/100‑1 (MQ=255)

                                     cGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAgg                                                                >  1:2435321/1‑100 (MQ=255)

                                     cGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAgg                                                                >  1:1331207/1‑100 (MQ=255)

                                           tAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGATCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTg                                                          <  2:2538044/100‑1 (MQ=255)

                                             aCCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTgcg                                                        <  2:2053694/100‑1 (MQ=255)

                                              cccTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTgcgc                                                       <  2:2183850/100‑1 (MQ=255)

                                                 tGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTAt                                                    <  2:338137/100‑1 (MQ=255)

                                                 tGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTAt                                                    <  2:3136560/100‑1 (MQ=255)

                                                 tGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTAt                                                    <  1:3679383/100‑1 (MQ=255)

                                                    aattaattTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCa                                                 <  1:1863975/100‑1 (MQ=255)

                                                    aattaattTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCa                                                 >  1:3344727/1‑100 (MQ=255)

                                                      ttaattTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCaaa                                               >  1:3797978/1‑100 (MQ=255)

                                                       taattTGCATAGTGGCAATTTTTTGCCAGACTGAACAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAac                                              <  2:784249/100‑1 (MQ=255)

                                                          tttGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAAcaca                                           >  2:2662168/1‑100 (MQ=255)

                                                          tttGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAAcaca                                           >  2:675077/1‑100 (MQ=255)

                                                             tccgatctgCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAAt                                        <  2:1766476/92‑1 (MQ=255)

                                                              cATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATc                                       <  1:2658509/100‑1 (MQ=255)

                                                              cATAGTGGCAATCTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATc                                       >  1:2514868/1‑100 (MQ=255)

                                                               aTAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCa                                      >  1:259720/1‑100 (MQ=255)

                                                                 aGTGGCAATGTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCaaa                                    <  2:592218/100‑1 (MQ=255)

                                                                  gTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCaacg                                   <  2:1706277/100‑3 (MQ=255)

                                                                  gTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAg                                   <  1:2206308/100‑1 (MQ=255)

                                                                   tggcaaaaaattGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGt                                  >  2:2201346/11‑100 (MQ=255)

                                                                    ggCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGtt                                 >  2:1263847/1‑100 (MQ=255)

                                                                     gCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAAtagatcgga                                >  1:1766486/1‑92 (MQ=255)

                                                                     gCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTc                                >  2:2624402/1‑100 (MQ=255)

                                                                      cAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCg                               >  1:784155/1‑100 (MQ=255)

                                                                         ttttttGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTgg                            <  2:2489802/100‑1 (MQ=255)

                                                                         ttttttGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTgg                            <  1:975896/100‑1 (MQ=255)

                                                                         tcttctGCCAGACTGAAGAGGTCATACCAGTTATTACCTCTGTACTTATAACCACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTgg                            <  2:2309432/95‑1 (MQ=255)

                                                                          tttttGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGGTCGTGGa                           <  1:2001815/100‑1 (MQ=255)

                                                                           ttttGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGat                          <  2:307400/100‑1 (MQ=255)

                                                                            tttGCCAGACTGAAGAGGTTATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGata                         <  1:3085243/100‑1 (MQ=255)

                                                                               gCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTATTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATatca                      <  2:3205698/100‑1 (MQ=255)

                                                                               gCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATatca                      >  1:3461408/1‑100 (MQ=255)

                                                                                tataagtacagagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATatcat                     >  2:361085/10‑100 (MQ=37)

                                                                                ccAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATatcat                     >  2:1373339/1‑100 (MQ=255)

                                                                                   gACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCAtctc                  >  1:1196547/1‑100 (MQ=255)

                                                                                    aCTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCg                 >  2:1622196/1‑100 (MQ=255)

                                                                                       gAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACg              >  2:146358/1‑100 (MQ=255)

                                                                                       gAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACg              >  2:774391/1‑100 (MQ=255)

                                                                                       gAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACg              >  2:3601751/1‑100 (MQ=255)

                                                                                       gAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACg              <  1:1959207/100‑1 (MQ=255)

                                                                                       acagagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGGTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACg              >  1:1281328/3‑100 (MQ=38)

                                                                                        cagagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGt             >  1:2857429/2‑100 (MQ=39)

                                                                                        cagagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGt             >  2:3499159/2‑100 (MQ=39)

                                                                                        aaGAGGTCATACCAGTTATGACCTCTGTACTTATAACGACAATGTAAGGTTATGGCGCTATGCAAACACGAATCAAAGTTCGTGGATATCATCTCGACGt             <  2:1277359/100‑1 (MQ=255)

                                                                                         agagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGtt            >  2:3560450/1‑100 (MQ=255)

                                                                                            ggTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTAc         >  1:1077953/1‑100 (MQ=255)

                                                                                            ggTCATACCAGTTATGACCTCTGTACGTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTAc         >  2:754129/1‑100 (MQ=255)

                                                                                             gTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTAcc        <  1:2893614/100‑1 (MQ=255)

                                                                                              tCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCa       <  2:2716050/100‑1 (MQ=255)

                                                                                              tCATACCAGTTATGACCTCTGTACTTAGAACAACAACGTAAGGTTATTGCTCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCa       <  2:2824457/100‑1 (MQ=255)

                                                                                                 tACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCa    <  1:3423335/100‑1 (MQ=255)

                                                                                                 tACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCa    >  2:3739118/1‑100 (MQ=255)

                                                                                                 tACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCa    <  2:3438783/100‑1 (MQ=255)

                                                                                                 tACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCa    <  2:2956608/100‑1 (MQ=255)

                                                                                                   ccAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACg  >  1:3820075/1‑100 (MQ=255)

                                                                                                   ccAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACg  <  1:1482720/100‑1 (MQ=255)

CCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACG  >  NC_000913/464255‑464453

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