BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,142,917	T→C	D225G (GAC→GGC)	rne ←	RNase E

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,142,917	0	T	C	87.0%	542.9 / 70.8	223	D225G (GAC→GGC)	rne	RNase E
Reads supporting (aligned to +/- strand): ref base T (17/12); new base C (97/97); total (114/109)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.02e-01

Only 100 of 227 total aligned reads displayed.

TTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTG  >  NC_000913/1142827‑1143009

 tgctgcTGAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCt                                                                                           >  1:1874512/1‑93 (MQ=255)

   ctgctgAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGa                                                                                         >  2:613512/1‑93 (MQ=255)

   ctgctgAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGa                                                                                         <  2:1503829/93‑1 (MQ=255)

    tgctgAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGAc                                                                                        <  1:45221/93‑1 (MQ=255)

    tgctgAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGAc                                                                                        <  1:1308251/93‑1 (MQ=255)

     gctgAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACg                                                                                       >  1:7481054/1‑93 (MQ=255)

        gAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACgtaa                                                                                    <  2:5054246/93‑1 (MQ=255)

        gAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACgtaa                                                                                    <  2:1521108/93‑1 (MQ=255)

         aaaTCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACgtaag                                                                                   <  1:1969868/93‑1 (MQ=255)

         aaaTCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACgtaag                                                                                   <  1:2464817/93‑1 (MQ=255)

          aaTCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACgtaagt                                                                                  >  1:888067/1‑93 (MQ=255)

               gggCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATcg                                                                             <  2:6612040/93‑1 (MQ=255)

               gggCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATcg                                                                             <  2:4153666/93‑1 (MQ=255)

                ggCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATcgc                                                                            <  2:740915/93‑1 (MQ=255)

                   gACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGa                                                                         >  1:7200239/1‑93 (MQ=255)

                    aCCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGaa                                                                        <  2:5659694/93‑1 (MQ=255)

                     ccTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAAt                                                                       >  1:6794245/1‑93 (MQ=255)

                     ccTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAAt                                                                       >  2:6779252/1‑93 (MQ=255)

                     ccTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAAt                                                                       >  1:472177/1‑93 (MQ=255)

                      cTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATg                                                                      >  1:4285778/1‑93 (MQ=255)

                       tAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATgc                                                                     >  2:4936904/1‑93 (MQ=255)

                          tGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGcgcg                                                                  >  2:1231359/1‑93 (MQ=255)

                          tGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGcgcg                                                                  >  2:2861129/1‑93 (MQ=255)

                           gCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGa                                                                 <  1:6850161/93‑1 (MQ=255)

                           gCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGa                                                                 <  2:2203376/93‑1 (MQ=255)

                            cAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGaa                                                                <  1:648955/93‑1 (MQ=255)

                            cAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCAGGTTATCGATAAGGATTTCGCCTATGCCCTGACGTAAGTAATCGCGGAATGCGCGaa                                                                <  2:3352119/93‑1 (MQ=255)

                             aGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAAc                                                               <  1:4910348/93‑1 (MQ=255)

                             aGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAAc                                                               <  2:2336636/93‑1 (MQ=255)

                             aGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAAc                                                               <  1:1122479/93‑1 (MQ=255)

                                gATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGAt                                                            <  2:5916123/93‑1 (MQ=255)

                                 atatGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGAtt                                                           >  2:2622144/1‑93 (MQ=255)

                                 atatGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGAtt                                                           <  1:1237174/93‑1 (MQ=255)

                                  tatGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTa                                                          <  2:7190841/93‑1 (MQ=255)

                                  tatGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTa                                                          >  1:7056198/1‑93 (MQ=255)

                                  tatGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTa                                                          >  2:3843686/1‑93 (MQ=255)

                                   atGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTAc                                                         >  1:4793363/1‑93 (MQ=255)

                                    tGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACg                                                        >  2:4491548/1‑93 (MQ=255)

                                     gCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGt                                                       <  2:54946/93‑1 (MQ=255)

                                      cTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGtt                                                      <  1:988385/93‑1 (MQ=255)

                                       tGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTg                                                     >  1:5659268/1‑93 (MQ=255)

                                        gACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGc                                                    >  1:6678547/1‑93 (MQ=255)

                                        gACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGc                                                    <  2:3637525/93‑1 (MQ=255)

                                        gACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGc                                                    <  2:911880/93‑1 (MQ=255)

                                          cGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGctc                                                  <  2:7130338/93‑1 (MQ=255)

                                          cGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGctc                                                  >  2:347825/1‑93 (MQ=255)

                                          cGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGctc                                                  <  1:306923/93‑1 (MQ=255)

                                           gTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGctct                                                 >  1:2325254/1‑93 (MQ=255)

                                            tGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCtctc                                                >  2:5027204/1‑93 (MQ=255)

                                            tGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCACTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCtctc                                                >  1:3399961/1‑93 (MQ=255)

                                               cAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCtg                                             >  2:4359381/1‑93 (MQ=255)

                                               cAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCtg                                             >  1:5657216/1‑93 (MQ=255)

                                                aGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCtga                                            <  2:7504023/93‑1 (MQ=255)

                                                aGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCtga                                            <  2:1973158/93‑1 (MQ=255)

                                                 gTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCtgat                                           >  1:361769/1‑93 (MQ=255)

                                                   tCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCtgatga                                         <  1:397916/93‑1 (MQ=255)

                                                    cGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGaa                                        <  1:2720288/93‑1 (MQ=255)

                                                     gAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAAt                                       >  1:6924818/1‑93 (MQ=255)

                                                     gAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAAt                                       >  1:2610207/1‑93 (MQ=255)

                                                     gAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAAGCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAAt                                       >  1:5145357/1‑93 (MQ=255)

                                                      aGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATc                                      >  1:3192701/1‑93 (MQ=255)

                                                      aGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATc                                      >  2:458226/1‑93 (MQ=255)

                                                       gCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCa                                     <  2:7038547/93‑1 (MQ=255)

                                                       gCACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCa                                     >  1:6657166/1‑93 (MQ=255)

                                                        cACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAg                                    >  2:62810/1‑93 (MQ=255)

                                                        cACTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAg                                    >  1:216809/1‑93 (MQ=255)

                                                          cTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGa                                  <  2:2543992/93‑1 (MQ=255)

                                                          cTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGa                                  <  1:4940791/93‑1 (MQ=255)

                                                          cTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGa                                  <  1:4530441/93‑1 (MQ=255)

                                                          cTTTCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGa                                  <  1:152340/93‑1 (MQ=255)

                                                           tttCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGaa                                 <  2:6711103/93‑1 (MQ=255)

                                                           tttCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGaa                                 >  2:1803925/1‑93 (MQ=255)

                                                           tttCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGaa                                 >  2:7043640/1‑93 (MQ=255)

                                                            ttCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAAc                                >  2:6508343/1‑93 (MQ=255)

                                                             tCGGGTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACg                               <  2:3872999/93‑1 (MQ=255)

                                                               gggTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACggg                             <  1:3701955/93‑1 (MQ=255)

                                                               gggTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACggg                             >  2:1436181/1‑93 (MQ=255)

                                                               gggTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACggg                             >  2:3236872/1‑93 (MQ=255)

                                                               gggTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACggg                             <  1:2855538/93‑1 (MQ=255)

                                                                 gTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGc                           >  1:622807/1‑93 (MQ=255)

                                                                 gTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGc                           <  1:552374/93‑1 (MQ=255)

                                                                 gTTATCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGc                           >  2:1036854/1‑93 (MQ=255)

                                                                    aTCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCgg                        <  2:5293462/93‑1 (MQ=255)

                                                                    aTCGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCgg                        >  2:2977074/1‑93 (MQ=255)

                                                                      cGATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGggc                      <  2:2786709/93‑1 (MQ=255)

                                                                       gATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGggcg                     >  1:5785626/1‑93 (MQ=255)

                                                                       gATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGggcg                     <  2:5324173/93‑1 (MQ=255)

                                                                       gATAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGggcg                     >  1:3130627/1‑93 (MQ=255)

                                                                         tAAGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGggcggc                   >  2:3034591/1‑93 (MQ=255)

                                                                          aaGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCt                  >  2:4088481/1‑93 (MQ=255)

                                                                           aGGATTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCtt                 <  2:5248619/93‑1 (MQ=255)

                                                                              aTTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCa              <  2:3109880/93‑1 (MQ=255)

                                                                              aTTTCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCa              >  1:4480468/1‑93 (MQ=255)

                                                                               tttCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAg             >  1:73242/1‑93 (MQ=255)

                                                                               tttCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAg             >  2:2466323/1‑93 (MQ=255)

                                                                                 tCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCg           <  2:1878445/93‑1 (MQ=255)

                                                                                 tCGCCGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCg           >  1:2814801/1‑93 (MQ=255)

                                                                                     cGATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCtt       <  1:6393711/93‑1 (MQ=255)

                                                                                      gATGCCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCttt      >  1:1041793/1‑93 (MQ=255)

                                                                                       aTGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCtttt     >  2:3279950/1‑93 (MQ=255)

TTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCTCTCCTGATGAATCAGGAACGGGGCCGGGCGGCTTTCAGCGGCTTTTTTG  >  NC_000913/1142827‑1143009

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