BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,617,966	T→C	N125N (AAT→AAC)	marA →	MarA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,617,966	0	T	C	99.4%	996.9 / ‑3.8	306	N125N (AAT→AAC)	marA	MarA
Reads supporting (aligned to +/- strand): ref base T (0/2); new base C (143/161); total (143/163)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

Only 100 of 311 total aligned reads displayed.

CCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCT  >  NC_000913/1617876‑1618055

ccGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCa                                                                                         >  1:5030490/1‑93 (MQ=255)

ccGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAAACa                                                                                         >  2:1857772/1‑93 (MQ=255)

 cGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCAt                                                                                        <  2:2656842/93‑1 (MQ=255)

  gAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCAtt                                                                                       >  2:1641422/1‑93 (MQ=255)

   aaCCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTa                                                                                      <  2:3244759/93‑1 (MQ=255)

   aaCCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTa                                                                                      >  1:2525827/1‑93 (MQ=255)

    aCCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTac                                                                                     <  2:4027359/93‑1 (MQ=255)

     ccTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTaca                                                                                    >  2:7102485/1‑93 (MQ=255)

      cTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTacaa                                                                                   <  1:7013091/93‑1 (MQ=255)

      cTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTacaa                                                                                   >  1:9385071/1‑93 (MQ=255)

      cTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTacaa                                                                                   <  1:4320776/93‑1 (MQ=255)

       ttCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTacaac                                                                                  <  1:3506246/93‑1 (MQ=255)

       ttCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTacaac                                                                                  >  1:8932949/1‑93 (MQ=255)

        tCAAAAATTCCTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTacaaca                                                                                 <  2:3463998/93‑1 (MQ=255)

        tCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTacaaca                                                                                 >  1:6313066/1‑93 (MQ=255)

         cAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAg                                                                                <  2:5535451/93‑1 (MQ=255)

          aaaaaTTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGc                                                                               <  1:5298842/93‑1 (MQ=255)

          aaaaaTTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGc                                                                               <  2:758801/93‑1 (MQ=255)

               ttACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGtt                                                                          <  2:9524498/93‑1 (MQ=255)

               ttACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGtt                                                                          <  2:8911260/93‑1 (MQ=255)

                   tttGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGaaa                                                                      >  1:9812555/1‑93 (MQ=255)

                    ttGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGaaaa                                                                     <  2:8041344/93‑1 (MQ=255)

                     tGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAAc                                                                    >  2:7206778/1‑93 (MQ=255)

                      gATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACg                                                                   <  2:5504338/93‑1 (MQ=255)

                      gATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACg                                                                   >  1:4079567/1‑93 (MQ=255)

                       aTGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGt                                                                  <  2:6504380/93‑1 (MQ=255)

                       aTGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGt                                                                  >  1:8641311/1‑93 (MQ=255)

                         gTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGa                                                                <  1:10008338/93‑1 (MQ=255)

                          ttCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGAc                                                               >  1:2746883/1‑93 (MQ=255)

                          ttCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGAc                                                               >  1:6208511/1‑93 (MQ=255)

                          ttCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGAc                                                               >  2:9080111/1‑93 (MQ=255)

                          tgcggccgcATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGAc                                                               <  2:4673607/89‑1 (MQ=255)

                            ccgccgCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACaa                                                             <  1:2105458/93‑1 (MQ=255)

                            ccgccgCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACaa                                                             <  1:8122879/93‑1 (MQ=255)

                             cgccgcATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAAc                                                            <  1:7640208/93‑1 (MQ=255)

                              gccgcATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACA‑CCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGt                                                          >  2:7245345/1‑93 (MQ=255)

                               ccgcATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGt                                                          <  1:8505685/93‑1 (MQ=255)

                                 gcATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCa                                                        <  1:8693690/93‑1 (MQ=255)

                                 gcATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCa                                                        <  1:7529764/93‑1 (MQ=255)

                                 gcATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCa                                                        >  2:2456335/1‑93 (MQ=255)

                                     aaaTACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGa                                                    <  1:4886247/93‑1 (MQ=255)

                                       aTACCGGATGACCAATATGCAGGGCGAATCGAGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAgg                                                  >  1:184517/1‑93 (MQ=255)

                                         aCCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCa                                                >  2:4615339/1‑93 (MQ=255)

                                          ccGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCaa                                               <  2:2393090/93‑1 (MQ=255)

                                          ccGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCaa                                               <  2:9659049/93‑1 (MQ=255)

                                           cGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAAt                                              <  1:3618988/93‑1 (MQ=255)

                                            ggATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCAc                                                       <  2:390168/83‑1 (MQ=255)

                                            ggATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATc                                             <  2:4408739/93‑1 (MQ=255)

                                             gATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCa                                            <  1:3556956/93‑1 (MQ=255)

                                              aTGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCAt                                           <  2:7891108/93‑1 (MQ=255)

                                                gACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGa                                         <  2:2849885/93‑1 (MQ=255)

                                                 aCCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGaa                                        <  1:6505385/93‑1 (MQ=255)

                                                  ccAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGaaa                                       <  1:8328892/93‑1 (MQ=255)

                                                   cAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAAc                                      >  2:1333492/1‑93 (MQ=255)

                                                    aaTATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAAcc                                     <  2:2292526/93‑1 (MQ=255)

                                                    aaTATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAAcc                                     <  1:4380018/93‑1 (MQ=255)

                                                    aaTATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAAcc                                     >  1:1611014/1‑93 (MQ=255)

                                                    aaTATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAAcc                                     <  2:3064277/93‑1 (MQ=255)

                                                    aaTATGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAAcc                                     <  1:684672/93‑1 (MQ=255)

                                                     atatGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTATAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCa                                    <  2:5751905/93‑1 (MQ=255)

                                                     atatGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCa                                    <  1:6029696/93‑1 (MQ=255)

                                                      tatGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAGCCAc                                   >  1:5882533/1‑93 (MQ=255)

                                                      tatGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCAc                                   >  2:4242198/1‑93 (MQ=255)

                                                       atGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACt                                  <  1:6634635/93‑1 (MQ=255)

                                                        tGCAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACtt                                 <  2:1461410/93‑1 (MQ=255)

                                                          cAGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAAAGTGACAACGTCACTGAGGCAATCATGAAACCACTTTc                               >  2:4964835/1‑93 (MQ=255)

                                                           aGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCa                              <  2:9170062/93‑1 (MQ=255)

                                                           aGGGCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCa                              >  2:6055710/1‑93 (MQ=255)

                                                           aGGGCGAATCGCAATTTTTACATCCATTAAGCCAGTACAGCAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCa                              <  2:220027/93‑1 (MQ=255)

                                                            gggCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCAt                             >  2:6727034/1‑93 (MQ=255)

                                                             ggCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATc                            <  1:3931480/93‑1 (MQ=255)

                                                             ggCGAATCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATc                            >  1:2009826/1‑93 (MQ=255)

                                                                 aaTCGCGCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCa                        <  1:5113349/93‑1 (MQ=255)

                                                                    cgcgCTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATa                     <  2:5679522/93‑1 (MQ=255)

                                                                     gctcTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATag                    <  1:3377823/90‑1 (MQ=255)

                                                                     gcgcTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATag                    >  2:449085/1‑93 (MQ=255)

                                                                      cgcTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATagc                   <  1:2991034/93‑1 (MQ=255)

                                                                       gcTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATagca                  >  2:8678256/1‑93 (MQ=255)

                                                                        cTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATagcag                 >  2:6175106/1‑93 (MQ=255)

                                                                        cTTTTTACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATagcag                 <  1:5403564/93‑1 (MQ=255)

                                                                         tttttACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATagcagc                >  2:6887724/1‑93 (MQ=255)

                                                                         tttttACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATagcagc                >  2:1050744/1‑93 (MQ=255)

                                                                         tttttACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATagcagc                >  1:3831600/1‑93 (MQ=255)

                                                                         tttttACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCAGCAATagcagc                >  1:3441378/1‑93 (MQ=255)

                                                                          ttttACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCt               >  2:8229940/1‑93 (MQ=255)

                                                                          ttttACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCt               >  2:5577896/1‑93 (MQ=255)

                                                                          ttttACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCt               >  2:4998433/1‑93 (MQ=255)

                                                                          ttttACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCt               >  1:4601346/1‑93 (MQ=255)

                                                                          ttttACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCt               >  1:6703898/1‑93 (MQ=255)

                                                                          ttttACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCt               <  2:4151258/93‑1 (MQ=255)

                                                                             tACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTgcg            >  2:7714414/1‑93 (MQ=255)

                                                                             tACATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTgcg            <  2:3750466/93‑1 (MQ=255)

                                                                              aCATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTgcgc           >  1:1497585/1‑93 (MQ=255)

                                                                               cATCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCt          <  1:167162/93‑1 (MQ=255)

                                                                                aTCCATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCtt         >  2:5274601/1‑93 (MQ=255)

                                                                                  ccATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTAt       <  1:5480309/93‑1 (MQ=255)

                                                                                  ccATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTAt       <  1:9126694/93‑1 (MQ=255)

                                                                                   cATTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTAtt      <  2:2987748/93‑1 (MQ=255)

                                                                                    aTTAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATtc     >  1:7296769/1‑93 (MQ=255)

                                                                                      tAAACCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATtctc   <  1:5096346/93‑1 (MQ=255)

CCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCT  >  NC_000913/1617876‑1618055

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