BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	464,358	T→G	15.8%	intergenic (+50/‑44)	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,358	0	T	G	15.8%	276.5 / 45.4	108	intergenic (+50/‑44)	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 T (38/53); new base G (5/12); total (43/65)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

Only 100 of 109 total aligned reads displayed.

ATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAA  >  NC_000913/464261‑464457

aTGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTa                                                                                                   <  1:2532897/100‑1 (MQ=255)

aTGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTa                                                                                                   <  1:2888454/100‑1 (MQ=255)

aTGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTa                                                                                                   <  2:1089466/100‑1 (MQ=255)

aTGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTa                                                                                                   <  2:1640208/100‑1 (MQ=255)

  gggCAATTCGCTGGTGGAACGTAATCTTGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGTTCATACCAGTTATg                                                                                                 <  2:737152/100‑1 (MQ=255)

  gggCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATg                                                                                                 <  1:67147/100‑1 (MQ=255)

  gggCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATg                                                                                                 <  2:2193415/100‑1 (MQ=255)

  gggCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATg                                                                                                 <  2:3040910/100‑1 (MQ=255)

  gggCAATTCGCTGGTGGAACGTAATCGGGCGGTATTAACCCCGTAATTAATTTGCATAGTGGCACTTTATTGCCAGACTGAAGAGGTCATACCAGTTAt                                                                                                  >  2:803703/1‑99 (MQ=255)

    gCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGATGTCATACGAGTTATg                                                                                                 >  2:805146/1‑98 (MQ=255)

    gCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGAc                                                                                               <  2:2215234/100‑1 (MQ=255)

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

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

       aTTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCCGACTGAAGAGGTCATACCAGTTATGACctc                                                                                            <  2:2047748/100‑1 (MQ=255)

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

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

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

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

         tCGCTGGTGGAATGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTATTGCCAGACTGAAGAGGTCATAACTGGTATGACctctt                                                                                          <  2:858029/100‑2 (MQ=38)

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

          cGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctcttc                                                                                         <  2:740619/100‑3 (MQ=39)

            cTGGTTGAACGTAATTTGGCGTTATTTTCTGTGTAATTAATTTTTATAGTGGCTATTTTTTGCCAGACGGAAGAGGTCATAACTGGTATGACctcttc                                                                                         <  2:1851571/98‑3 (MQ=16)

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

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

                  gAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctcttcagtctggc                                                                                 <  1:512346/100‑11 (MQ=25)

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

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

                   aaCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaac                                                                                >  2:1710059/1‑100 (MQ=255)

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

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

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

                          tCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACGACAACg                                                                         <  2:1193339/100‑1 (MQ=255)

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

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

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

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

                               cGGTATTAACTCTGTAATTAATTTGCATAGTGGCAATTTTTTCCCAGACTGAAGAGGTCACACCAGTTATGCCCTCTGTACTTATAACAACAACGTAAgg                                                                    <  2:1228361/100‑1 (MQ=255)

                               cGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACTAGTTATGACCTCTGTACTTATAACAACAACGTAAgg                                                                    <  1:2196768/100‑1 (MQ=255)

                               cGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAgg                                                                    <  1:1167487/100‑1 (MQ=255)

                                 gTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGt                                                                   <  2:168562/99‑1 (MQ=255)

                                  tATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTa                                                                 <  1:1973247/100‑1 (MQ=255)

                                     tAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTg                                                              <  1:15974/100‑1 (MQ=255)

                                     tAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTg                                                              <  1:1592395/100‑1 (MQ=255)

                                     tAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTg                                                              <  2:2607275/100‑1 (MQ=255)

                                      aaCCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTgc                                                             <  1:2865366/100‑1 (MQ=255)

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

                                       aCCCTGTAATTAATTTGCATAGTGGCAATTTTTTGACATATTGAAAAGGACATATCAGTTATGACTTCTGTACTTGTACCAACAACGTAAGtttctat                                                              >  2:254573/1‑94 (MQ=255)

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

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

                                        cccTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGCTATTgcgc                                                           >  1:2081797/1‑100 (MQ=255)

                                          cTGTATTTAAGTTGCATAGTGGCAATTTTTTGCCATACTTAAGAGTTCATACCAGTTATTACTTCTTTACTTATAACAACAGCGTAAGGTTATTGCGCt                                                          <  2:143199/99‑1 (MQ=255)

                                          cTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTa                                                         >  2:2046492/1‑100 (MQ=255)

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

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

                                                  aattTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAaca                                                 >  1:116771/1‑100 (MQ=255)

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

                                                    tttGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAAcaca                                               <  2:1908061/100‑1 (MQ=255)

                                                      ctcATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACaaa                                             <  1:1132402/98‑1 (MQ=255)

                                                        cATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACaaaag                                           >  2:1132397/1‑98 (MQ=255)

                                                            gTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAg                                       >  2:1317250/1‑100 (MQ=255)

                                                              ggCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGtt                                     <  1:1189529/100‑1 (MQ=255)

                                                                  aTTTTTTGCCTGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTg                                 <  2:2035639/100‑1 (MQ=255)

                                                                  aTTTTTTGCCAGACTGAAGAGGTCATACCAGTGATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTg                                 >  1:2946537/1‑100 (MQ=255)

                                                                  aTTTTTTGCCAGACTCAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTg                                 <  1:3061308/100‑1 (MQ=255)

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

                                                                    tttttGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGa                               <  1:2541663/100‑1 (MQ=255)

                                                                     ttttGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGat                              >  1:2847539/1‑100 (MQ=255)

                                                                         gCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATatca                          <  1:827069/100‑1 (MQ=255)

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

                                                                             gACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCAtctc                      <  2:1766151/100‑1 (MQ=255)

                                                                             gACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCAtctc                      <  1:3082426/100‑1 (MQ=255)

                                                                              aCTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCAtctc                      >  1:916900/1‑99 (MQ=255)

                                                                               gtacagagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGATATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGa                    >  2:2992727/5‑100 (MQ=38)

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

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

                                                                                 gAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATAAAAGTTCGTGGATATCATCTCGACg                  >  1:1099199/1‑100 (MQ=255)

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

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

                                                                                  aaGAGGTCATACGAGTTATGACCTCTTTACTTCTAACAAAAACCTAAGCTTATTGTGCTATGCAAACCCAATTCAAAGTTCGTGGATATCATCTAGACg                  <  2:1919445/99‑1 (MQ=255)

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

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

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

                                                                                    gagGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGttt               >  1:1778745/1‑100 (MQ=255)

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

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

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

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

                                                                                        tCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCATACACAAATCAAAGTTCGTTGATATCATCTCGACGTTTACCa           >  2:2477686/1‑100 (MQ=255)

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

                                                                                         cATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAg          >  2:2508524/1‑100 (MQ=255)

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

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

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

                                                                                            aCCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCAc       <  1:877100/100‑1 (MQ=255)

                                                                                            aCCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCAc       <  2:1777262/100‑1 (MQ=255)

                                                                                               aGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTc    >  1:2057502/1‑100 (MQ=255)

                                                                                               aGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTc    >  2:526538/1‑100 (MQ=255)

                                                                                               aGTTATGACCTCTGTACTTATAACAAAAACGTAAGGTTATTGCGCTATGCAAACACAAAGCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTc    <  2:2762159/100‑1 (MQ=255)

                                                                                                gTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTca   <  1:2084547/100‑1 (MQ=255)

                                                                                                 ttATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTcaa  <  2:1556212/100‑1 (MQ=255)

ATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAA  >  NC_000913/464261‑464457

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