BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,133,750	C→T	intergenic (+244/‑85)	metF → / → katG	5,10‑methylenetetrahydrofolate reductase/hydroperoxidase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,133,750	0	C	T	100.0%	1272.4 / NA	304	intergenic (+244/‑85)	metF/katG	5,10‑methylenetetrahydrofolate reductase/hydroperoxidase I
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (152/152); total (152/152)

Only 100 of 304 total aligned reads displayed.

TTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACGTCAG  >  NC_000913/4133652‑4133847

tttATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAATCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTa                                                                                                 <  2:5449433/101‑1 (MQ=255)

  tataAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATc                                                                                               >  1:5449925/1‑101 (MQ=255)

  tataAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATc                                                                                               <  2:5449488/101‑1 (MQ=255)

  tataAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATc                                                                                               >  2:5449928/1‑101 (MQ=255)

    taAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGc                                                                                             >  1:5449932/1‑101 (MQ=255)

    taAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGc                                                                                             <  1:5449479/101‑1 (MQ=255)

      aaGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCAt                                                                                           >  1:5449936/1‑101 (MQ=255)

       aGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATc                                                                                          >  2:5449937/1‑101 (MQ=255)

         cATTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCg                                                                                        <  2:5449680/101‑1 (MQ=255)

          aTTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGt                                                                                       >  1:5449941/1‑101 (MQ=255)

           tttGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTg                                                                                      >  2:5449943/1‑101 (MQ=255)

             tGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGa                                                                                    <  1:5449534/101‑1 (MQ=255)

                ccGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTa                                                                                 <  2:5449249/101‑1 (MQ=255)

                 cGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTaa                                                                                <  2:5449692/101‑1 (MQ=255)

                  gAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAAt                                                                               >  1:5449948/1‑101 (MQ=255)

                   aaTTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAAtt                                                                              >  2:5449950/1‑101 (MQ=255)

                   aaTTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAAtt                                                                              >  1:5449952/1‑101 (MQ=255)

                    aTTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTc                                                                             <  2:5449551/101‑1 (MQ=255)

                    aTTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTc                                                                             <  2:5449572/101‑1 (MQ=255)

                     ttCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCa                                                                            <  1:5449587/101‑1 (MQ=255)

                       cGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAAt                                                                          >  1:5449956/1‑101 (MQ=255)

                         gACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATta                                                                        >  1:5449959/1‑101 (MQ=255)

                         gACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATta                                                                        >  1:5449962/1‑101 (MQ=255)

                          aCATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCa                                                                            >  1:5449965/1‑96 (MQ=255)

                          aCATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATtat                                                                       <  2:5449467/101‑1 (MQ=255)

                          aCATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATtat                                                                       <  1:5449409/101‑1 (MQ=255)

                          aCATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATtat                                                                       <  1:5449604/101‑1 (MQ=255)

                           cATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATtata                                                                      <  1:5449497/101‑1 (MQ=255)

                           cATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATtata                                                                      <  1:5449620/101‑1 (MQ=255)

                           cATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATtata                                                                      <  1:5449520/101‑1 (MQ=255)

                            aTAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATaa                                                                     <  1:5449688/101‑1 (MQ=255)

                            aTAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATaa                                                                     >  2:5449968/1‑101 (MQ=255)

                             tAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAAc                                                                    <  2:5449566/101‑1 (MQ=255)

                              aaTCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACt                                                                   <  1:5449344/101‑1 (MQ=255)

                              aaTCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACt                                                                   >  1:5449972/1‑101 (MQ=255)

                               aTCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACtt                                                                  <  1:5449810/101‑1 (MQ=255)

                               aTCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACtt                                                                  >  2:5449974/1‑101 (MQ=255)

                                tCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTtc                                                                 <  1:5449775/101‑1 (MQ=255)

                                 cAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTtct                                                                >  1:5449977/1‑101 (MQ=255)

                                   aaaaaGCTTAATTAAGATCAATTTGCTCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTtctct                                                              >  2:5449980/1‑101 (MQ=255)

                                    aaaaGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTtctctc                                                             >  1:5449983/1‑101 (MQ=255)

                                    aaaaGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTtctctc                                                             <  1:5449681/101‑1 (MQ=255)

                                     aaaGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTctctct                                                            >  2:5449985/1‑101 (MQ=255)

                                      aaGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTa                                                           <  2:5449521/101‑1 (MQ=255)

                                      aaGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTa                                                           <  1:5449643/101‑1 (MQ=255)

                                         cTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACg                                                        >  1:5449989/1‑101 (MQ=255)

                                         cTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACg                                                        <  2:5449511/101‑1 (MQ=255)

                                            aattaaGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCtg                                                     >  1:5449992/1‑101 (MQ=255)

                                             attaaGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCtgt                                                    <  1:5449445/101‑1 (MQ=255)

                                             attaaGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCtgt                                                    <  2:5449546/101‑1 (MQ=255)

                                               taaGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTgtgt                                                  <  2:5449741/101‑1 (MQ=255)

                                               taaGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTgtgt                                                  >  1:5449995/1‑101 (MQ=255)

                                                aaGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTa                                                 >  2:5449999/1‑101 (MQ=255)

                                                   aTCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCg                                              <  2:5449586/101‑1 (MQ=255)

                                                   aTCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCg                                              >  1:5450000/1‑101 (MQ=255)

                                                     cAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTa                                            >  1:5450005/1‑101 (MQ=255)

                                                      aaTTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTaa                                           >  1:5450008/1‑101 (MQ=255)

                                                        tttGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACg                                         >  2:5450011/1‑101 (MQ=255)

                                                          tGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGt                                       <  2:5449818/101‑1 (MQ=255)

                                                           gATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTa                                      <  2:5449646/101‑1 (MQ=255)

                                                           gATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTa                                      >  1:5450014/1‑101 (MQ=255)

                                                           gATCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTa                                      <  2:5449384/101‑1 (MQ=255)

                                                            aTCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTaa                                     >  1:5450017/1‑101 (MQ=255)

                                                             tCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGt                                       >  1:5450021/1‑98 (MQ=255)

                                                             tCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAac                                    <  1:5449552/101‑1 (MQ=255)

                                                             tCTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAac                                    >  2:5450019/1‑101 (MQ=255)

                                                              cTACATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAaca                                   >  1:5450025/1‑101 (MQ=255)

                                                                aCATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACt                                 >  1:5450028/1‑101 (MQ=255)

                                                                aCATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACt                                 <  1:5449629/101‑1 (MQ=255)

                                                                 cATCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTg                                >  2:5450031/1‑101 (MQ=255)

                                                                  aTCTCTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGt                               <  2:5449605/101‑1 (MQ=255)

                                                                   tctcTTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTa                              <  2:5449782/101‑1 (MQ=255)

                                                                    ctctTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTag                             <  1:5449658/101‑1 (MQ=255)

                                                                     tctTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTaga                            >  1:5450034/1‑101 (MQ=255)

                                                                      ctTTAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTagag                           <  1:5449594/101‑1 (MQ=255)

                                                                       tttAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAgg                          <  1:5449737/101‑1 (MQ=255)

                                                                       tttAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAgg                          <  2:5449704/101‑1 (MQ=255)

                                                                       tttAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAgg                          >  1:5450037/1‑101 (MQ=255)

                                                                        ttAACCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAggg                         <  1:5449901/101‑1 (MQ=255)

                                                                          aaCCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGa                       >  2:5450043/1‑101 (MQ=255)

                                                                          aaCCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGa                       >  1:5450042/1‑101 (MQ=255)

                                                                          aaCCAACAATATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGa                       <  2:5449720/101‑1 (MQ=255)

                                                                             caacaaTATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGca                    >  1:5450048/1‑101 (MQ=255)

                                                                              aacaaTATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGcac                   <  2:5449631/101‑1 (MQ=255)

                                                                                caaTATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACAt                 >  2:5450052/1‑101 (MQ=255)

                                                                                caaTATGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACAt                 >  1:5450051/1‑101 (MQ=255)

                                                                                  atatGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATtg               >  2:5450058/1‑101 (MQ=255)

                                                                                  atatGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATtg               >  1:5450057/1‑101 (MQ=255)

                                                                                  atatGTAAGATCTCAATTAGCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCCCATtg               <  2:5449348/101‑1 (MQ=255)

                                                                                    atGTAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGa                       <  1:5449718/91‑1 (MQ=255)

                                                                                     tttAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATtgatg            >  2:5450062/3‑101 (MQ=255)

                                                                                       tAAGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCGCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAg          <  1:5449576/101‑1 (MQ=255)

                                                                                         aGATCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCa        <  1:5449699/101‑1 (MQ=255)

                                                                                           aTCTCAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACg      >  1:5450066/1‑101 (MQ=255)

                                                                                             ctcAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAgggg                        >  1:5450070/1‑81 (MQ=255)

                                                                                             ctcAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACGTc    <  2:5449618/101‑1 (MQ=255)

                                                                                             ctcAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACGTc    <  2:5449640/101‑1 (MQ=255)

                                                                                              tcAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACGTCa   >  1:5450073/1‑101 (MQ=255)

                                                                                              tcAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACGTCa   <  2:5449667/101‑1 (MQ=255)

                                                                                               cAATTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACGTCAg  <  2:5449580/101‑1 (MQ=255)

TTTATAAAGCATTTGTCCGAATTCGGACATAATCAAAAAAGCTTAATTAAGATCAATTTGATCTACATCTCTTTAACCAACAATATGTAAGATCTCAACTATCGCATCCGTGGATTAATTCAATTATAACTTCTCTCTAACGCTGTGTATCGTAACGGTAACACTGTAGAGGGGAGCACATTGATGAGCACGTCAG  >  NC_000913/4133652‑4133847

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