BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,767,507	(C)_7→8	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	93.1%	425.3 / 24.1	131	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 . (6/3); new base C (58/63); total (65/66)
Fisher's exact test for biased strand distribution p-value = 3.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 141 total aligned reads displayed.

ATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑T‑‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCT  >  NC_000913/2767401‑2767594

aTGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAt‑t‑t                                                                                                      <  1:1127774/100‑1 (MQ=255)

   atatTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑T‑‑c‑c‑c                                                                                               <  1:1599332/100‑1 (MQ=255)

       tAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTTACAACTTCCTGCATAGGTAACACGTTTTCCCCTATGCACCCACGACTAAAT‑T‑C‑‑c‑c‑c‑ccc                                                                                           >  2:921105/1‑99 (MQ=255)

        aaCTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGACAACCCGTTTTCCCCGGGGCACCCCCGACTAAAT‑T‑T‑‑c‑c‑c‑cccc                                                                                          >  1:1950787/1‑99 (MQ=255)

        aaCTTGTCCAATATGATCAAATAGCATTAACCACCCCTCACAACATCCTGCAGAGGGAGCACGTTTTCCGCTGTGTACCCCCGACTAGAt‑t‑tt                                                                                                     >  2:256290/1‑92 (MQ=255)

         aCTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑T‑‑c‑c‑c‑ccccct                                                                                        <  2:114307/100‑3 (MQ=255)

          cTTGTCCAATATGATCAAATAGCATTAACCCCCACTTACAACACCATGATGCGGTAACAAGCTCTCTCCTGTGCAACCACGACTACAT‑T‑c‑‑c‑c                                                                                                 >  2:1365719/1‑92 (MQ=255)

              tCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTaa‑a                                                                                   <  2:2284448/100‑1 (MQ=255)

               ccAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑Ac                                                                                  <  1:670521/100‑1 (MQ=255)

                cAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACt                                                                                 <  1:653193/100‑1 (MQ=255)

                 aaTATGATCAAATAGCATTAACCCCCCCTCACAACGTCATGCATAGGGAACACGTTTTCCCCTGTGAACACACTACTAAAT‑T‑T‑‑C‑C‑C‑CCCCCTTta‑a                                                                                   >  2:2520158/1‑96 (MQ=255)

                  atatGATCAAATAGTATTAACCCCCCATCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTtc                                                                               <  2:2686171/100‑1 (MQ=255)

                   tatGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTtct                                                                              <  1:2815406/100‑1 (MQ=255)

                       aTCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAAGGAACACGCTTTCCCCTGTGCACCCACAACTAAAT‑T‑T‑‑C‑C‑C‑CCCCATCAACACTTCTCaa                                                                           >  2:308017/1‑99 (MQ=255)

                       aTCAAATAGCATTAACCCCCCCTAAAGATGTCGTGAATAGGGAACACGATATCCCCTGAGCACAGACGAGTAAAT‑T‑T‑‑c‑c‑c‑cacattctc                                                                                     >  2:724612/1‑81 (MQ=255)

                       aTCAAATAGCATCAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCaat                                                                          <  2:1519118/100‑1 (MQ=255)

                         cAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCaataa                                                                        <  1:472811/100‑1 (MQ=255)

                          aaaTAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCaataat                                                                       <  1:2107186/100‑1 (MQ=255)

                          aaaTAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑T‑‑CAC‑C‑CCCCTTTAA‑ACTTCTCaataat                                                                       >  2:2382849/1‑100 (MQ=255)

                            aTAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑AGTTCTCAATAATGt                                                                     <  1:2095374/100‑1 (MQ=255)

                            aTAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGt                                                                     <  1:2229585/100‑1 (MQ=255)

                             tAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTc                                                                    <  1:2389511/100‑1 (MQ=255)

                               gCATTACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCAc                                                                  <  2:1951256/100‑1 (MQ=255)

                               gCATTAACCCCCCCTCACAACGTCCTTCATAGGGAACACATTTTCCCCCCTGCACCCACACCTAAAT‑T‑TC‑C‑C‑C‑CCCCTTCAA‑ACTTCTCAATACTGTc                                                                    >  2:2347188/1‑98 (MQ=255)

                               gCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCAc                                                                  <  2:1540068/100‑1 (MQ=255)

                               gCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCAc                                                                  <  1:863366/100‑1 (MQ=255)

                               gCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCAc                                                                  <  2:959229/100‑1 (MQ=255)

                               gCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGCTTTCCCCTGTGCACCCACGACTCAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTAAATAATGTCAc                                                                  >  2:633328/1‑100 (MQ=255)

                               gCATTAAACCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCAAGTCCACCCACGACTCAAT‑T‑A‑‑C‑CTC‑CGGCTTTAA‑CCTTCCCaataa                                                                        >  2:609503/1‑94 (MQ=255)

                                cATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCAcc                                                                 <  2:578713/100‑1 (MQ=255)

                                cATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCAcc                                                                 <  1:243565/100‑1 (MQ=255)

                                cATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCAcc                                                                 <  1:2775699/100‑1 (MQ=255)

                                 aTTAACCCCCCCTCACAACGTCCTGCCTAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTAAATAATGTAACCa                                                                >  2:661341/1‑100 (MQ=255)

                                    atcccccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACTACTAAAT‑TGC‑‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCa                                                                <  2:2964939/95‑1 (MQ=255)

                                      cccccccTCACAACGTCCTGCATAGGGCACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAAtt                                                             <  1:33925/98‑1 (MQ=255)

                                      cccccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAAtttt                                                           >  2:608173/1‑100 (MQ=255)

                                       ccccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTc                                                          <  1:300884/100‑1 (MQ=255)

                                        cccccTCACAACTTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCt                                                         <  1:2612372/100‑1 (MQ=255)

                                        cccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCt                                                         <  1:745065/100‑1 (MQ=255)

                                           ccTCACAACTGCCCGCATAGGGAACACGTTTTCCCCTGTTCACCCACGACTAAATGT‑C‑‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGGCACCAATTGTCTg                                                        <  2:1158533/98‑1 (MQ=255)

                                             tCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGaaa                                                    <  2:1608878/100‑1 (MQ=255)

                                              cacaACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAAt                                                   >  2:112004/1‑100 (MQ=255)

                                               acaACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATc                                                  >  2:822777/1‑100 (MQ=255)

                                               acaACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATc                                                  >  1:2265053/1‑100 (MQ=255)

                                                caACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATcc                                                 <  1:2519123/100‑1 (MQ=255)

                                                caACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATcc                                                 >  1:2640848/1‑100 (MQ=255)

                                                    gTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatg                                             >  2:1994724/1‑100 (MQ=255)

                                                    gTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatg                                             >  1:1315266/1‑100 (MQ=255)

                                                     tCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatga                                            >  1:2744513/1‑100 (MQ=255)

                                                       cTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgatg                                          >  2:2804667/1‑100 (MQ=255)

                                                       cTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgatg                                          <  1:2318726/100‑1 (MQ=255)

                                                       cTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgatg                                          <  1:2320254/100‑1 (MQ=255)

                                                        tGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGc                                         <  1:980598/100‑1 (MQ=255)

                                                        tGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGc                                         <  1:2914270/100‑1 (MQ=255)

                                                        tGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTAACCAATTTTCTGAAAATCCTatgct                                           >  2:2565289/1‑96 (MQ=255)

                                                         gCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGcc                                        <  2:2561169/100‑1 (MQ=255)

                                                         gCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGcc                                        >  1:135703/1‑100 (MQ=255)

                                                         gCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGcc                                        <  1:1927829/100‑1 (MQ=255)

                                                         gCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGcc                                        <  2:1663742/100‑1 (MQ=255)

                                                          cATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCt                                       >  1:1473067/1‑100 (MQ=255)

                                                          cATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCt                                       >  1:1754987/1‑100 (MQ=255)

                                                            tAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCttt                                     >  1:734748/1‑100 (MQ=255)

                                                            tAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCttt                                     >  1:2882505/1‑100 (MQ=255)

                                                             aGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc                                    >  1:1211183/1‑100 (MQ=255)

                                                              gggAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑T‑‑C‑C‑CTCCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTcc                                   >  2:1684721/1‑100 (MQ=255)

                                                               ggAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTccc                                  <  1:1016761/100‑1 (MQ=255)

                                                                gAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCg                                 >  1:1630992/1‑100 (MQ=255)

                                                                 aaCACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TCCC‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCg                                 >  1:2058495/1‑100 (MQ=255)

                                                                 aaCACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCACTAATGTCACCAATTTTATGGAAATCCAATGATACCTTTCCAgg                                >  1:254825/1‑100 (MQ=255)

                                                                 aaCACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCgg                                <  2:677841/100‑1 (MQ=255)

                                                                  acacGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCGCCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGt                               <  1:2327840/100‑1 (MQ=255)

                                                                  acacGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGt                               >  1:419939/1‑100 (MQ=255)

                                                                  acacGTTTTCCCCGGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑GCTTCTAAAGAATATCACCAATTTTCTGGAAATCCTATGAGGCCTTTCCCGGt                               >  1:3075650/1‑100 (MQ=255)

                                                                   cacGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TCCC‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTccc                                  >  2:346405/1‑97 (MQ=255)

                                                                   cacGTTTTCCCCTGCGCACCCATGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTa                              >  2:2383620/1‑100 (MQ=255)

                                                                      ttccgatcttgtgCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACTAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAAt                           <  2:1895987/91‑1 (MQ=255)

                                                                      gTTTTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAAt                           >  1:720563/1‑100 (MQ=255)

                                                                      gTTCTCCCCTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAAt                           <  1:2617902/100‑1 (MQ=255)

                                                                           ccccTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGc                      >  1:974313/1‑100 (MQ=255)

                                                                           ccccTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGc                      >  1:1021148/1‑100 (MQ=255)

                                                                           ccccTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGc                      <  1:2646248/100‑1 (MQ=255)

                                                                            cccTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGc                      >  2:54941/1‑99 (MQ=255)

                                                                              cTGTGCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgt                   >  2:354016/1‑100 (MQ=255)

                                                                               tgtgCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCATTAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgtg                  <  2:2169253/100‑1 (MQ=255)

                                                                               tgtgCACCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAAtagatcggaa                  >  1:1896008/1‑91 (MQ=255)

                                                                                    aCCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTc             >  1:704199/1‑100 (MQ=255)

                                                                                    aCCCACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTc             <  2:135703/100‑1 (MQ=255)

                                                                                     cccACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTc             <  2:3061454/99‑1 (MQ=255)

                                                                                     cccACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCa            <  2:1211173/100‑1 (MQ=255)

                                                                                     cccACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCa            >  2:3068108/1‑100 (MQ=255)

                                                                                       cACGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGa          >  1:2582791/1‑100 (MQ=255)

                                                                                         cGACTAAGT‑T‑TG‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGaaa        >  2:1621599/1‑100 (MQ=255)

                                                                                         cGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGaaa        <  2:2965015/100‑1 (MQ=255)

                                                                                         cGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGaaa        <  1:633330/100‑1 (MQ=255)

                                                                                         cGACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGaaa        >  1:19403/1‑100 (MQ=255)

                                                                                          gACTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAg       >  2:5330/1‑100 (MQ=255)

                                                                                           aCTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTCCTCAATAATGCCACCAATTTTCTGGAAATCCTATGAAGCCTTTCCCCATAAATGGTGCTGTAAGTTCAGaa         >  2:2902306/1‑97 (MQ=255)

                                                                                            cTAAAT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGAt     <  2:121543/100‑1 (MQ=255)

                                                                                              aaaT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTc   >  1:1023776/1‑100 (MQ=255)

                                                                                               aaT‑T‑TC‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCt  <  2:1630974/100‑1 (MQ=255)

ATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑T‑T‑‑C‑C‑C‑CCCCTTTAA‑ACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCT  >  NC_000913/2767401‑2767594

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: `‑`