BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	464,354	C→A	16.4%	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	16.4%	409.4 / 61.9	162	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 (64/70); new base A (8/19); total (72/90)
Fisher's exact test for biased strand distribution p-value = 9.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.26e-01

Only 100 of 171 total aligned reads displayed.

CGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACG  >  NC_000913/464256‑464453

cGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTATCCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAcc                                                                                                    <  2:2147344/100‑1 (MQ=255)

cGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAcc                                                                                                    <  1:1563844/100‑1 (MQ=255)

   ggATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGt                                                                                                 >  1:2764373/1‑100 (MQ=255)

     aTGGGCAATTCGCTGGTGGAACGTAATATGGTAGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTa                                                                                               <  1:1773588/100‑1 (MQ=255)

      tGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTAt                                                                                              >  2:2443467/1‑100 (MQ=255)

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

         gCAATTCGTTGGTGTAACTTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATg                                                                                             <  1:809111/98‑1 (MQ=255)

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

          cAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGAcc                                                                                          >  2:948923/1‑100 (MQ=255)

           aaTTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACct                                                                                         <  2:2950332/100‑1 (MQ=255)

            aTTCTCTGGTGGAACTTAATCTGGCGTTTTTAACTCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATg                                                                                             <  1:2274001/95‑1 (MQ=255)

            aTTCGCTGGTTGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTTCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctc                                                                                        <  1:1981441/100‑1 (MQ=255)

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

            aTTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTAAAACCAGTTATGACctc                                                                                        >  1:3866488/1‑100 (MQ=255)

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

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

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

               ctcTGGTGTAACGTAATCTGGCGGCTTTAACCCTGTAATTTATTTGCATAGTTTCAACCTTTTGCCATCCTGAAGATTTCATAACTGCTATGACctcct                                                                                      <  2:557298/97‑3 (MQ=255)

               cGCTGGTGGAACGTAATCTGTCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGt                                                                                     <  1:1872431/100‑1 (MQ=255)

               cGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTTCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctcttc                                                                                     <  1:348082/100‑3 (MQ=38)

               cGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGTCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctcttc                                                                                     <  1:1010504/100‑3 (MQ=38)

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

               cGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTAGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGt                                                                                     >  1:1077827/1‑100 (MQ=255)

                gCTGGTGGAACGTAATCTGGCGGTATTAATCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctcttca                                                                                    <  1:2076247/100‑4 (MQ=38)

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

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

                  tggtggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACt                                                                                  <  1:3801898/100‑1 (MQ=255)

                  tggtggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACt                                                                                  <  2:3332317/100‑1 (MQ=255)

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

                    gtggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTta                                                                                <  2:186057/100‑1 (MQ=255)

                     tggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTtat                                                                               <  1:1145926/100‑1 (MQ=255)

                     tggAACGTAATCTGGCGGTATTAACCATGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTtat                                                                               >  1:2374747/1‑100 (MQ=255)

                      ggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTtata                                                                              >  1:1349307/1‑100 (MQ=255)

                       gAACGTAATCTGGCTGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaa                                                                             <  1:595365/100‑1 (MQ=255)

                        aaCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaan                                                                            <  2:1436317/100‑2 (MQ=255)

                          cGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATGGTGGCAATTTTTTGCCAGCATGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaacaa                                                                          >  2:2742994/1‑100 (MQ=255)

                          cGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGATTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaacaa                                                                          <  1:984583/100‑1 (MQ=255)

                           gTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaacaac                                                                         <  2:2977357/100‑1 (MQ=255)

                             aaTCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaacaacaa                                                                       >  2:3462410/1‑100 (MQ=255)

                               tCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACg                                                                     <  1:2500583/100‑1 (MQ=255)

                               tCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACg                                                                     >  2:3882287/1‑100 (MQ=255)

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

                                  ggCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAAGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTaa                                                                  >  1:2192319/1‑100 (MQ=255)

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

                                   gCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAg                                                                 >  1:723134/1‑100 (MQ=255)

                                     ggTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGt                                                               <  2:2337363/100‑1 (MQ=255)

                                      gTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGtt                                                              <  1:1738967/100‑1 (MQ=255)

                                           aaCCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTgc                                                         >  2:1371128/1‑100 (MQ=255)

                                           aaCCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTgc                                                         >  1:2231800/1‑100 (MQ=255)

                                            aCCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTgcg                                                        >  1:407089/1‑100 (MQ=255)

                                                 gTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATg                                                   >  1:856197/1‑100 (MQ=255)

                                                 gTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATg                                                   >  2:1000599/1‑100 (MQ=255)

                                                  taattaatTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGc                                                  <  1:2702493/100‑1 (MQ=255)

                                                  taattaatTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACAGTTATGAACTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGc                                                  <  2:1080666/100‑1 (MQ=255)

                                                     ttaattTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCaaa                                               <  1:2024031/100‑1 (MQ=255)

                                                        attTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAacac                                            <  1:2952192/100‑1 (MQ=255)

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

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

                                                                   ggcaaaaaattGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGtt                                 >  2:848961/10‑100 (MQ=255)

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

                                                                    gCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGTAAACACAAATCAAAGTTc                                >  2:3189759/1‑100 (MQ=255)

                                                                      aaTTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGt                              >  2:1659012/1‑100 (MQ=255)

                                                                       aTTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTg                             <  2:615806/100‑1 (MQ=255)

                                                                         tttttGCCAGACTGAAGAGGTCATACAAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGATATACAAACACAAATCAAAGTTCGTgg                            >  2:1445775/1‑99 (MQ=255)

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

                                                                           tttGCCAGACTGAAGAGGTCATACCAGTTATGTCCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGata                         <  1:2722256/100‑1 (MQ=255)

                                                                           tttGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGata                         <  2:1196622/100‑1 (MQ=255)

                                                                            ttGCCAGACTGAAGAGGTTATACTAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGatat                        <  1:3541024/100‑1 (MQ=255)

                                                                             tgacAGACTGAAGAGGTCATACCAGTGATGAGCTCTGTACTTATAACAACAACGTAAAGTTATTGCCCTATGCAAACACGAATCAAAGTTCGTGGatat                        <  2:3562097/96‑1 (MQ=255)

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

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

                                                                                cAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATatcatc                    >  2:3123104/1‑100 (MQ=255)

                                                                                cAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATatcatc                    >  1:1699921/1‑100 (MQ=255)

                                                                                cAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATatcat                     >  1:3608433/1‑99 (MQ=255)

                                                                                 aGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCAtct                   <  1:2474539/100‑1 (MQ=255)

                                                                                    cTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGa                >  2:1784993/1‑100 (MQ=255)

                                                                                     tacagagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCACTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCg                 >  2:258267/4‑98 (MQ=38)

                                                                                     tGAAGAGGTCATATCAGTTATGACCTCTGTACTTATAATAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGAc               <  1:3330270/100‑1 (MQ=255)

                                                                                     tGAAGAGGTCATACCAGTTACGACCTCTGTACTTATACCAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCTAATTTCGTGGATATCATCTCGAc               <  2:567760/100‑1 (MQ=255)

                                                                                      acagagGTCATAACTGGTATGACCTCTGTACTAATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGAc               >  1:3728774/3‑99 (MQ=38)

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

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

                                                                                       aaGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGt             <  2:2284737/100‑1 (MQ=255)

                                                                                         gagGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGttt           <  2:2433831/100‑1 (MQ=255)

                                                                                         gagGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGctt           >  2:1846570/1‑100 (MQ=255)

                                                                                         gagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGttt           >  2:2359763/1‑100 (MQ=255)

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

                                                                                            gTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATACAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTAcc        >  1:2879746/1‑100 (MQ=255)

                                                                                            gTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGATATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTAcc        <  1:1306797/100‑1 (MQ=255)

                                                                                             tCATACCAGTTATGATCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCa       <  2:395701/100‑1 (MQ=255)

                                                                                             tCATACCAGTTATGACCTTTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTAcc        <  1:2780872/99‑1 (MQ=255)

                                                                                             tCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCa       >  1:1453367/1‑100 (MQ=255)

                                                                                              cATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGAAGTTTACCa       >  2:3757581/1‑99 (MQ=255)

                                                                                               aTACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGc     <  2:1217238/100‑1 (MQ=255)

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

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

                                                                                                 aCCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCAc   >  2:1199934/1‑100 (MQ=255)

                                                                                                  ccAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCTTGGATATCATCTCGACGTTTACCAGCACg  <  2:596888/100‑1 (MQ=255)

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

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

CGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACG  >  NC_000913/464256‑464453

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