BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,767,507:1	(C)_7→8	100%	intergenic (+152/‑203)	rnlB → / → yfjP	CP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; 50S ribosome‑binding GTPase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,767,500	1	.	C	85.6%	315.4 / 52.6	111	intergenic (+145/‑210)	rnlB/yfjP	CP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; 50S ribosome‑binding GTPase family protein
Reads supporting (aligned to +/- strand): ref base . (8/8); new base C (43/49); total (54/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 119 total aligned reads displayed.

GATATTAACTTGTCCAATATGATCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTT‑‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTC  >  NC_000913/2767403‑2767593

gATATTAACTTGTCCAATATGATCAAATAGCATTAA‑CCCCCCCTTCCAACATCATCCATAAGAAACACATTTTCCCCTGCGCACCCGCCACAA‑AATTT‑‑c‑c                                                                                                 >  1:2869011/1‑100 (MQ=255)

 atatTAACTTGTCCAATATGATCAAATAGCATTAA‑TAAACCCTACCAACTTCCTGCAAAGGGAATAAGTCTTCCCCTCTGCACCCACGACAA‑AAtct‑‑c                                                                                                   >  1:891988/1‑98 (MQ=255)

 acatTAACTTGTCCAATATGATCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTT‑‑c‑c‑c                                                                                               <  1:1853851/98‑1 (MQ=255)

   atTAACTTGTCCAATATGATCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTT‑‑c‑c‑cc‑c                                                                                            <  1:3026293/100‑1 (MQ=255)

     tAACTTGTCCAATATGATCAAATAGCATTAA‑CTCACTCTTTCAACGTCCTGCATAAGGAATACTTTTTCCCCTGTGCACCCACGACTA‑AATTT‑‑c‑a‑ccatc                                                                                           >  1:302248/1‑97 (MQ=255)

      aaCTTGTCCAATATGATCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTT‑‑c‑c‑cc‑c‑ccc                                                                                        <  1:2294571/100‑2 (MQ=255)

       aCTTGTCCAATATGATCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTT‑‑c‑c‑cc‑c‑ccct                                                                                       <  2:2938817/100‑3 (MQ=255)

       aCTTGTCCAATATGATCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGGTTGCCCCTGTGCACCCCCGACTA‑AATGT‑‑c‑c‑cc‑c‑ccc                                                                                        <  1:1884654/99‑2 (MQ=255)

          tGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTT‑‑C‑C‑CC‑C‑CC‑cttt                                                                                    <  2:622584/100‑2 (MQ=255)

          tGTCCAATATGATCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTT‑‑C‑C‑CC‑C‑CC‑cttta                                                                                   <  2:1837435/100‑3 (MQ=255)

           gTCCAATATGATCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTaa                                                                                   <  2:870422/100‑1 (MQ=255)

           gTCCAATATGATCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTaa                                                                                   <  2:349739/100‑1 (MQ=255)

           gTCCAATATGATCAAATAGCATTAA‑CCCCCCCTCACAACATCATGCATAGGGAACACGTTTTCCACAGTGCACTCACGACTA‑AATTTA‑C‑C‑CC‑C‑CC‑TTTaa                                                                                   >  1:228912/1‑100 (MQ=255)

             ccAATATGATCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAAc                                                                                 <  2:853018/100‑1 (MQ=255)

                  atGATCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTtctc                                                                            <  1:125167/100‑1 (MQ=255)

                   tGATCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTCCC‑C‑CC‑C‑CC‑TTTAAACTtct                                                                             <  2:2405831/99‑1 (MQ=255)

                    gATCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCaa                                                                          <  1:143130/100‑1 (MQ=255)

                     aTCAAATAGCATTAC‑CCCCCCCCCACAACGCCCTCCATAGGCAACACGTTTGCCCCCCTGCACCCACGAATA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAAcccccc                                                                            >  1:2173331/1‑92 (MQ=255)

                        aaaTAGCATTAA‑CCCCCCCTCACTACGTCCTTCAAAAGGAACACGTGTTCCCCTGTACACCCACGACTC‑AATTT‑‑CAC‑CC‑C‑CC‑CTTCACCGTCTCaata                                                                        >  1:2192794/1‑98 (MQ=255)

                        aaaTAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAATCTTCTCaataat                                                                      <  2:2007464/100‑1 (MQ=255)

                        aaaTAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCaataat                                                                      <  2:942232/100‑1 (MQ=255)

                        aaaTAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAAGGAACACATTTTCCCCTGTGCACCCACAAATA‑AATTT‑‑C‑C‑CC‑C‑CC‑CTTAAACTCCCCAATAATg                                                                     >  1:1580254/1‑100 (MQ=255)

                         aaTAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATGTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCaataa                                                                       <  1:918774/99‑1 (MQ=255)

                            aGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCa                                                                  <  2:2858/100‑1 (MQ=255)

                            aGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCa                                                                  <  2:37532/100‑1 (MQ=255)

                             gCATTAA‑CCCCCCCTCACAACTTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCCCGACTA‑AATTT‑‑C‑C‑CC‑CACC‑TTTAAAATTCTCcataatct                                                                    >  1:324965/1‑95 (MQ=255)

                             gCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTG‑C‑C‑CC‑C‑CC‑TTTAAACTTCCCAATAATGt                                                                    >  1:1233154/1‑97 (MQ=255)

                              cATTAA‑CCCCCCCTCACAACGTCCTCCATAGGGAACACGTTTTCACCTGTGCACCCACGACTA‑AATTT‑‑C‑CGCC‑C‑CC‑TTTAACCTTCTACATAATGTCAcc                                                                >  1:2995736/1‑100 (MQ=255)

                               aTTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTCCTCAATCATGACACCa                                                               >  2:1587025/1‑100 (MQ=255)

                                  aa‑CCCCCCCTCACACCGTCCTGCATCGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCCCAATAAAGTTAcc                                                                >  1:1300517/1‑96 (MQ=255)

                                       cccccTCACAACGTCCTGCATAGGAAACACGTTTTCCCCTGTGCACCCACTACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCt                                                        <  1:1511759/100‑1 (MQ=255)

                                         cccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTgg                                                      >  1:1431844/1‑100 (MQ=255)

                                          ccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGa                                                     <  1:2336063/100‑1 (MQ=255)

                                           cacccaACGCCCTGCATAGGGAACACGTTTTCCCCTGTGTCCATACGTCTACTACTC‑‑C‑C‑CC‑C‑CC‑TTTCAACTTCTCAATAATGTCTGCAATTGTCTg                                                       <  1:251892/93‑1 (MQ=255)

                                           cTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGaa                                                    >  2:94027/1‑100 (MQ=255)

                                           cTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGaa                                                    >  2:1806473/1‑100 (MQ=255)

                                            tCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGaaa                                                   >  2:185367/1‑100 (MQ=255)

                                             cacaACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACCACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAAt                                                  <  1:2565059/100‑1 (MQ=255)

                                             cacaACGTCCTGCATACGGAACACCTTTTCCGCTGTGCACCCACGACTA‑AAATT‑‑C‑CGCC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGaaa                                                   <  1:1097180/99‑1 (MQ=255)

                                              acaACGTCCTGCATCGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTACACTTCTCAAGAATGTCACCAATATTCTGAAAATc                                                 >  1:527294/1‑100 (MQ=255)

                                              acaACGTCCTGCATCGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCGCAATAATGTCACCAATTTTCTGGaa                                                    >  1:822212/1‑97 (MQ=255)

                                              acaACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATc                                                 <  2:891993/100‑1 (MQ=255)

                                               caACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATcc                                                >  2:1695419/1‑100 (MQ=255)

                                                  cGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTa                                              >  2:611934/1‑99 (MQ=255)

                                                   gTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatg                                            <  2:1805161/100‑1 (MQ=255)

                                                   gTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatg                                            <  1:1717517/100‑1 (MQ=255)

                                                    tCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatga                                           <  1:564964/100‑1 (MQ=255)

                                                     ccTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgat                                          >  1:2104220/1‑100 (MQ=255)

                                                     ccTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgat                                          <  2:1424418/100‑1 (MQ=255)

                                                       tGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGc                                        >  2:330402/1‑100 (MQ=255)

                                                       tGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGc                                        >  1:141618/1‑100 (MQ=255)

                                                        gCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AGTTTC‑C‑C‑CC‑C‑CC‑TTTACACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatga                                           >  1:2875620/1‑96 (MQ=255)

                                                        gCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGcc                                       >  2:2353340/1‑100 (MQ=255)

                                                        gCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGcc                                       >  2:2834045/1‑100 (MQ=255)

                                                        gCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGcc                                       >  1:348336/1‑100 (MQ=255)

                                                         cATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCCATTTTCTGGAAATCCTATGATGCCt                                      >  2:718264/1‑100 (MQ=255)

                                                          aTAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCtt                                     >  1:1484026/1‑100 (MQ=255)

                                                           tAGGTAACCCGTTGTCCCCTGTGCACCCACTACTA‑AATTC‑‑C‑CTCC‑C‑CC‑GTTAAACTGCGCAATAATGTCACCAATTTTCTGGAAATCCTatgatg                                         <  1:3012936/95‑1 (MQ=255)

                                                           tAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCttt                                    >  2:1818324/1‑100 (MQ=255)

                                                           tAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTAACCAATTTTCTGGAAATCCTATGATGCCttt                                    >  1:2894318/1‑100 (MQ=255)

                                                              ggAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTccc                                 <  2:328253/100‑1 (MQ=255)

                                                              ggAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTccc                                 <  2:151299/100‑1 (MQ=255)

                                                              ggAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTGTCTGGAAATCCTATGATGCCTTTccc                                 <  1:735369/100‑1 (MQ=255)

                                                              ggAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTT‑‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCg                                >  2:2842439/1‑100 (MQ=255)

                                                               gAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCg                                >  1:1402232/1‑100 (MQ=255)

                                                                aaCACGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCgg                               >  2:537335/1‑100 (MQ=255)

                                                                  cacGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTa                             <  2:422628/100‑1 (MQ=255)

                                                                  cacGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTa                             >  2:2479243/1‑100 (MQ=255)

                                                                  cacGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTa                             >  1:1273709/1‑100 (MQ=255)

                                                                   acGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTaa                            >  2:1567813/1‑100 (MQ=255)

                                                                   acGTTTTCCCCTGTGCACCAACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCTATTTTCTTTATATCCTATGATGCCTTTCCCGGTa                             <  1:1018769/99‑1 (MQ=255)

                                                                    cGTTTTCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTaaa                           <  2:76280/100‑1 (MQ=255)

                                                                    cGTTTTCCCCTGTGCACCCACGACTA‑AATTT‑‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAAt                          >  1:642648/1‑100 (MQ=255)

                                                                      ttttCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATg                         <  1:2428214/100‑1 (MQ=255)

                                                                       nttCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATgg                        >  2:2036218/2‑100 (MQ=255)

                                                                         tCCCCTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTg                      <  2:2238463/100‑1 (MQ=255)

                                                                          ccccTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTACACTTCTCAATAATCTCACAACTTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGc                     >  2:1987900/1‑100 (MQ=255)

                                                                          ccccTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGc                     <  1:412280/100‑1 (MQ=255)

                                                                          ccccTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGc                     <  1:733067/100‑1 (MQ=255)

                                                                          ccccTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGc                     <  1:2289141/100‑1 (MQ=255)

                                                                           cccTGTGCCCCCTCGACTA‑AATTT‑‑C‑C‑CCGC‑CC‑TTTAAACTTCTCATTAATGTCACCAATTTTCTGTAAATCCTATGATGCCTTTCCCGGTAAATGGTGc                     <  1:1836334/99‑1 (MQ=255)

                                                                           cccTGTGCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGc                     <  1:2758549/99‑1 (MQ=255)

                                                                            ccTGTGCACCCACGACTA‑AATTT‑‑C‑C‑CC‑C‑CCATTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtg                   <  1:1466067/100‑1 (MQ=255)

                                                                            ccTGTGCACCCACGACTA‑AATTT‑‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTaagatcggaaga                  >  2:1694793/1‑90 (MQ=255)

                                                                              tgtgCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgtg                 <  2:1300523/100‑1 (MQ=255)

                                                                              tgtgCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgcg                 <  2:2076089/100‑3 (MQ=255)

                                                                               gtgCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgtg                 <  2:1702191/99‑1 (MQ=255)

                                                                                 gCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGt              >  1:213628/1‑100 (MQ=255)

                                                                                 gCACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGt              >  2:1343757/1‑100 (MQ=255)

                                                                                  cACCCCCGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTAAAACTTCTCAACAATTTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATgctgcttttattt             >  1:2266337/1‑93 (MQ=255)

                                                                                  cACCCACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGtt             <  2:2680337/100‑1 (MQ=255)

                                                                                    cccACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAGTAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCa           >  2:2149248/1‑100 (MQ=255)

                                                                                     ccACGACTA‑ACTTTC‑C‑C‑CC‑C‑CC‑TGTGAACTTCTCAATAATGTCACCAGTTTTCTGGAAATCCTATGATGCCTTTCCCGGTACATGGTGCTGTGAGTTCa           <  1:826566/99‑1 (MQ=255)

                                                                                     ccACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAg          <  2:2661300/100‑1 (MQ=255)

                                                                                     ccACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAg          <  1:2914804/100‑1 (MQ=255)

                                                                                      cACGACTA‑AATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGa         >  1:1376505/1‑100 (MQ=255)

                                                                                       actaCTA‑AATTT‑‑C‑CTCC‑C‑CC‑TTTAATCTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAg          <  1:536726/95‑1 (MQ=255)

                                                                                       aCGACTA‑TATTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAGTTTTCTGGAAAGCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTc            <  1:2341768/96‑1 (MQ=255)

                                                                                             a‑aaTTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGAtt   >  1:718825/1‑99 (MQ=255)

                                                                                             a‑aaTTTC‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAATTCAGAAAGATTc  >  2:2234183/1‑100 (MQ=255)

GATATTAACTTGTCCAATATGATCAAATAGCATTAA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTA‑AATTT‑‑C‑C‑CC‑C‑CC‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTC  >  NC_000913/2767403‑2767593

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