BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	342,821	C→T	F238F (TTC→TTT)	yahK →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	342,821	0	C	T	99.6%	918.2 / ‑5.8	246	F238F (TTC→TTT)	yahK	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (132/113); total (132/114)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 255 total aligned reads displayed.

GCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCTGGTTGG  >  NC_000913/342732‑342910

gCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGAt                                                                                        <  2:1601485/93‑1 (MQ=255)

 cAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGAtt                                                                                       <  2:2471956/93‑1 (MQ=255)

   aaaCGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTc                                                                                     <  2:5895038/93‑1 (MQ=255)

     aCGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGCTTTGATTTCAt                                                                                   >  2:5621426/1‑93 (MQ=255)

     aCGCGAAGAGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGCGATGGCGGCTCATCTGCAGAGTTTTGATTTCAt                                                                                   >  1:2170785/1‑93 (MQ=255)

      cgcgAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCAtt                                                                                  <  2:5982611/93‑1 (MQ=255)

       gcgAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCAttt                                                                                 <  2:3771559/93‑1 (MQ=255)

        cgAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCAtttt                                                                                <  1:3166753/93‑1 (MQ=255)

        cgAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCAtttt                                                                                <  1:1783017/93‑1 (MQ=255)

         gAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTg                                                                               >  2:1488261/1‑93 (MQ=255)

            gCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAAt                                                                            >  1:5424365/1‑93 (MQ=255)

             cGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATa                                                                           >  2:1558198/1‑93 (MQ=255)

               gCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACa                                                                         >  2:4494155/1‑93 (MQ=255)

               gCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACa                                                                         <  2:5228976/93‑1 (MQ=255)

                cAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAg                                                                        >  2:2095848/1‑93 (MQ=255)

                 aaaaGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGt                                                                       <  1:5730730/93‑1 (MQ=255)

                  aaaGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTa                                                                      >  2:4250347/1‑93 (MQ=255)

                    aGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGc                                                                    >  2:5758270/1‑93 (MQ=255)

                    aGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCCGAAGAGTTTTGATTTCATTTTGAATACAGTAGc                                                                    >  1:6595/1‑93 (MQ=255)

                     gCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCt                                                                   <  2:3645094/93‑1 (MQ=255)

                     gCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCt                                                                   >  2:5413226/1‑93 (MQ=255)

                       ccTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTgc                                                                 <  1:5271199/93‑1 (MQ=255)

                       ccTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGGTTTGATTTCATTTTGAATACAGTAGATgc                                                                 >  1:3351192/1‑93 (MQ=255)

                        cTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTGGCTgcg                                                                >  1:421510/1‑93 (MQ=255)

                        cTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTgcg                                                                >  2:3193878/1‑93 (MQ=255)

                          gggggCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGcc                                                              <  2:260299/93‑1 (MQ=255)

                          gggggCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGcc                                                              >  1:3753538/1‑93 (MQ=255)

                            gggCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCcac                                                            >  2:6734590/1‑93 (MQ=255)

                             ggCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCcaca                                                           >  1:202704/1‑93 (MQ=255)

                              gCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACAt                                                          >  2:5306726/1‑93 (MQ=255)

                                 gATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATaaa                                                       >  2:6436004/1‑92 (MQ=255)

                                 gATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAAt                                                       >  2:653263/1‑93 (MQ=255)

                                  aTGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAAtc                                                      <  1:1155024/93‑1 (MQ=255)

                                  aTGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAAtc                                                      <  1:2859745/93‑1 (MQ=255)

                                  aTGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAAtc                                                      <  1:1841875/93‑1 (MQ=255)

                                   tGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAAtct                                                     >  1:5347146/1‑93 (MQ=255)

                                   tGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAAtct                                                     >  1:1213916/1‑93 (MQ=255)

                                    gAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAAtttc                                                    >  1:4118549/1‑93 (MQ=255)

                                    gAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAAtctc                                                    <  1:5846111/93‑1 (MQ=255)

                                    gAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAAtctc                                                    >  1:6047296/1‑93 (MQ=255)

                                    gAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAAtctc                                                    >  2:4991790/1‑93 (MQ=255)

                                       gttgttATCTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTcgac                                                 <  2:4002391/93‑1 (MQ=255)

                                       gttgttAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTcgac                                                 >  2:6431969/1‑93 (MQ=255)

                                        ttgttAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTcgacg                                                >  1:5018918/1‑93 (MQ=255)

                                        ttgttAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTcgacg                                                <  1:2379904/93‑1 (MQ=255)

                                          gttAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGAt                                              >  1:6720575/1‑93 (MQ=255)

                                           ttAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGAtt                                             <  1:6383943/93‑1 (MQ=255)

                                           ttAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGAtt                                             >  1:3828341/1‑93 (MQ=255)

                                             aaCTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGAtttt                                           <  2:4383276/93‑1 (MQ=255)

                                              aCTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTa                                          <  2:677984/93‑1 (MQ=255)

                                              aCTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTa                                          >  2:2693827/1‑93 (MQ=255)

                                              aCTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTa                                          <  2:6157094/93‑1 (MQ=255)

                                              aCTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTa                                          >  1:5482075/1‑93 (MQ=255)

                                               cTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTac                                         <  1:2704889/93‑1 (MQ=255)

                                               cTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTac                                         <  2:6263787/93‑1 (MQ=255)

                                                tCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTacc                                        <  2:1308423/93‑1 (MQ=255)

                                                tCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTacc                                        >  1:4873246/1‑93 (MQ=255)

                                                tCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTacc                                        >  1:6227199/1‑93 (MQ=255)

                                                 cACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTacca                                       <  1:3944629/93‑1 (MQ=255)

                                                 cACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTacca                                       <  1:3078295/93‑1 (MQ=255)

                                                  aCGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTaccac                                      >  2:5806016/1‑93 (MQ=255)

                                                   cGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTaccacc                                     >  1:6677046/1‑93 (MQ=255)

                                                    gCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCt                                    <  2:2299704/93‑1 (MQ=255)

                                                    gCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCt                                    <  1:6269228/93‑1 (MQ=255)

                                                      aaTGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTg                                  >  2:486756/1‑93 (MQ=255)

                                                       aTGCCGATGAGATGGCGGCTCATCTGAAGATTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGc                                 <  2:3695374/93‑1 (MQ=255)

                                                       aTGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGc                                 <  2:3881465/93‑1 (MQ=255)

                                                       aTGCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGc                                 <  1:75580/93‑1 (MQ=255)

                                                         gCCGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCGTGCTg                               >  1:894922/1‑93 (MQ=255)

                                                          ccGATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGa                              <  2:2840048/93‑1 (MQ=255)

                                                            gATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAg                            >  2:6652529/1‑93 (MQ=255)

                                                            gATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAg                            <  1:2251038/93‑1 (MQ=255)

                                                            gATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAg                            >  2:2004599/1‑93 (MQ=255)

                                                            gATGAGATGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAg                            >  1:1530754/1‑93 (MQ=255)

                                                                agaTGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTg                        >  2:3459855/1‑93 (MQ=255)

                                                                 gaTGGCGGCTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGa                       >  1:1054047/1‑93 (MQ=255)

                                                                     gcggcTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGc                   <  1:3794618/93‑1 (MQ=255)

                                                                      cggcTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCa                  <  1:6548218/93‑1 (MQ=255)

                                                                        gcTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCAcc                >  2:3475038/1‑93 (MQ=255)

                                                                         cTCATCTGAATAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCa               <  1:5095464/93‑1 (MQ=255)

                                                                         cTCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCa               >  2:800078/1‑93 (MQ=255)

                                                                          tCATCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCAt              <  1:2752650/93‑1 (MQ=255)

                                                                            aTCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGa            <  1:2300946/93‑1 (MQ=255)

                                                                            aTCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGa            >  2:2966227/1‑93 (MQ=255)

                                                                            aTCTGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGa            <  1:4294854/93‑1 (MQ=255)

                                                                               tGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGc         >  1:2524151/1‑93 (MQ=255)

                                                                               tGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGc         >  2:2175763/1‑93 (MQ=255)

                                                                               tGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGc         >  2:1648150/1‑93 (MQ=255)

                                                                               tGAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGc         <  2:623606/93‑1 (MQ=255)

                                                                                gAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGTACCATGACGCt        >  2:3849043/1‑93 (MQ=255)

                                                                                gAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCt        >  1:4421475/1‑93 (MQ=255)

                                                                                gAAGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCt        >  1:1012792/1‑93 (MQ=255)

                                                                                 aaGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCtg       >  2:4650936/1‑93 (MQ=255)

                                                                                 aaGAGTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCtg       <  2:4531637/93‑1 (MQ=255)

                                                                                  agagTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCtgg      >  1:4091064/1‑93 (MQ=255)

                                                                                   gagTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCtggt     >  1:4525223/1‑93 (MQ=255)

                                                                                   gagTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCtggt     >  1:5592668/1‑93 (MQ=255)

                                                                                   gagTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCtggt     <  1:139963/93‑1 (MQ=255)

                                                                                   gagTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCtggt     >  2:4295066/1‑93 (MQ=255)

                                                                                     gTTTTGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCtggttg   >  2:1889305/1‑93 (MQ=255)

GCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCTGGTTGG  >  NC_000913/342732‑342910

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