BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,162,448	A→G	T411A (ACA→GCA)	btuB →	BtuB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,162,448	0	A	G	99.5%	729.6 / ‑4.5	211	T411A (ACA→GCA)	btuB	BtuB
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (109/101); total (109/102)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.94e-01

Only 100 of 216 total aligned reads displayed.

GATAACTCACAGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGACATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTGGGAAGGC  >  NC_000913/4162358‑4162537

  tAACTCACAGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATc                                                                                       >  2:450496/1‑93 (MQ=255)

   aaCTCACAGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCt                                                                                      <  2:856940/93‑1 (MQ=255)

   aaCTCACAGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCt                                                                                      <  1:176786/93‑1 (MQ=255)

   aaCTCACAGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCt                                                                                      <  2:166327/93‑1 (MQ=255)

   aaCTCACAGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCt                                                                                      >  1:6159330/1‑93 (MQ=255)

    aCTCACAGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCtt                                                                                     <  2:2968373/93‑1 (MQ=255)

      tCACAGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTtat                                                                                   <  1:5616025/93‑1 (MQ=255)

       cacaGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTtata                                                                                  <  1:6198072/93‑1 (MQ=255)

        acaGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATaa                                                                                 <  1:4152662/93‑1 (MQ=255)

          aGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAgg                                                                               >  2:5658471/1‑93 (MQ=255)

           gTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGc                                                                              >  2:5748003/1‑93 (MQ=255)

            tttGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCa                                                                             >  1:1022631/1‑93 (MQ=255)

               ggTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTCCCTACGGGGCATCTTATAAGGCACCa                                                                          >  2:1230570/1‑93 (MQ=255)

                gtcgtcATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCaa                                                                         >  2:4301833/1‑93 (MQ=255)

                gtcgtcATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCaa                                                                         >  2:3316593/1‑93 (MQ=255)

                 tcgtcATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCaaa                                                                        >  1:5170719/1‑93 (MQ=255)

                   gtcATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATc                                                                      <  1:4980215/93‑1 (MQ=255)

                     cATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTg                                                                    >  1:5622709/1‑93 (MQ=255)

                      aTGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTgg                                                                   <  1:2851746/93‑1 (MQ=255)

                      aTGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTgg                                                                   <  2:6459141/93‑1 (MQ=255)

                       tGGAACCTGGCAAACCAGCGCCGGGTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTggg                                                                  >  2:1686512/1‑93 (MQ=255)

                        ggAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTgggg                                                                 <  1:3888926/93‑1 (MQ=255)

                        ggAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTgggg                                                                 <  2:256377/93‑1 (MQ=255)

                        ggAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTgggg                                                                 >  1:6277861/1‑93 (MQ=255)

                         gacccTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGc                                                                <  1:3086900/90‑1 (MQ=255)

                         gAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGc                                                                <  2:2935267/93‑1 (MQ=255)

                          aaCCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCa                                                               <  1:3539352/93‑1 (MQ=255)

                           aCCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCaa                                                              >  1:5037306/1‑93 (MQ=255)

                            ccTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAAc                                                             >  2:4898656/1‑93 (MQ=255)

                             cTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACt                                                            <  2:6424187/93‑1 (MQ=255)

                              tGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTg                                                           <  1:1217652/93‑1 (MQ=255)

                                gCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTa                                                         <  1:6028635/93‑1 (MQ=255)

                                 cAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTAt                                                        <  1:2300405/93‑1 (MQ=255)

                                 cAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTAt                                                        >  1:700928/1‑93 (MQ=255)

                                  aaaCCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATg                                                       <  1:4265631/93‑1 (MQ=255)

                                  aaaCCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATg                                                       <  1:5976980/93‑1 (MQ=255)

                                  aaaCCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATg                                                       >  2:1086007/1‑93 (MQ=255)

                                    aCCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGc                                                     >  2:1780437/1‑93 (MQ=255)

                                     ccAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCt                                                    <  2:2089275/93‑1 (MQ=255)

                                      cAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCtt                                                   >  2:1088093/1‑93 (MQ=255)

                                       aGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTc                                                  <  1:4868207/93‑1 (MQ=255)

                                       aGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTc                                                  >  2:4428245/1‑93 (MQ=255)

                                       aGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTc                                                  >  2:900532/1‑93 (MQ=255)

                                        gcgcCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCt                                                 >  2:3212695/1‑93 (MQ=255)

                                         cgcCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTa                                                <  2:1340127/93‑1 (MQ=255)

                                         cgcCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGACATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTa                                                <  2:3045471/93‑1 (MQ=255)

                                          gcCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTAc                                               <  2:5219511/93‑1 (MQ=255)

                                          gcCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTAc                                               <  2:6403315/93‑1 (MQ=255)

                                           ccGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACg                                              >  1:2238457/1‑93 (MQ=255)

                                            cGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACgg                                             >  2:5582472/1‑93 (MQ=255)

                                             ggTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGa                                            <  2:497683/93‑1 (MQ=255)

                                             ggTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGa                                            <  2:2477082/93‑1 (MQ=255)

                                              gTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGaa                                           >  2:1109141/1‑93 (MQ=255)

                                                tGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAAt                                         <  1:5305287/93‑1 (MQ=255)

                                                 gggAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATc                                        >  1:3607003/1‑93 (MQ=255)

                                                  ggAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATcc                                       <  1:4559670/93‑1 (MQ=255)

                                                   gAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCg                                      >  2:4493202/1‑93 (MQ=255)

                                                   gAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCg                                      <  2:4598801/93‑1 (MQ=255)

                                                    aaTTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGa                                     <  1:5530975/93‑1 (MQ=255)

                                                     aTTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGaa                                    >  1:1990791/1‑93 (MQ=255)

                                                     aTTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGaa                                    <  1:2693521/93‑1 (MQ=255)

                                                     aTTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTCCTACGGAAATCCGaa                                    <  2:3109326/93‑1 (MQ=255)

                                                        cATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCt                                 <  1:2490095/93‑1 (MQ=255)

                                                        cATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCt                                 >  2:2716667/1‑93 (MQ=255)

                                                        cATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCt                                 >  2:4582536/1‑93 (MQ=255)

                                                        cATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCt                                 <  1:3836765/93‑1 (MQ=255)

                                                         aTCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTg                                <  2:3572711/93‑1 (MQ=255)

                                                         aTCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTg                                >  1:1518755/1‑93 (MQ=255)

                                                         aTCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAAGCACCAAATCTGGGGCAGCTGTATGGCTTCTACGGAAATCCAGATCTg                                >  1:3443389/1‑93 (MQ=255)

                                                            gAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGAc                             >  1:5803118/1‑93 (MQ=255)

                                                             aaGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGAcc                            <  2:2552824/93‑1 (MQ=255)

                                                             aaGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGAcc                            >  2:2348016/1‑93 (MQ=255)

                                                              aGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGAccc                           >  1:4708360/1‑93 (MQ=255)

                                                               ggTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCg                          <  1:548399/93‑1 (MQ=255)

                                                                 ttATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCggg                        >  2:479898/1‑92 (MQ=255)

                                                                  tATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGgag                       >  1:6430723/1‑93 (MQ=255)

                                                                   aTCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGgaga                      >  2:3641183/1‑93 (MQ=255)

                                                                   aTCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGgaga                      <  2:1175757/93‑1 (MQ=255)

                                                                     cGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGaaa                    >  2:5518256/1‑93 (MQ=255)

                                                                       cTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAg                  >  2:46489/1‑93 (MQ=255)

                                                                        ttCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGc                 >  2:185787/1‑93 (MQ=255)

                                                                         tCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCa                >  2:3795211/1‑93 (MQ=255)

                                                                           aTTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCaaa              >  2:124235/1‑93 (MQ=255)

                                                                           aTTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCaaa              >  1:2048642/1‑93 (MQ=255)

                                                                            ttGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAAc             >  1:349040/1‑93 (MQ=255)

                                                                            ttGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAAc             <  1:4005814/93‑1 (MQ=255)

                                                                            ttGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAAc             >  1:5419143/1‑93 (MQ=255)

                                                                            ttGCTTCCTACGGGGCATCTTATAAGACTCCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAAc             <  2:5265418/93‑1 (MQ=255)

                                                                              gCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAg           >  2:5058091/1‑93 (MQ=255)

                                                                              gCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAg           >  2:1466632/1‑93 (MQ=255)

                                                                               cTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGt          >  2:5814158/1‑93 (MQ=255)

                                                                                ttCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTg         >  2:2510534/1‑93 (MQ=255)

                                                                                ttCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTg         >  1:5773172/1‑93 (MQ=255)

                                                                                  ccTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTggg       <  2:3896862/93‑1 (MQ=255)

                                                                                  ccTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTggg       <  1:1448347/93‑1 (MQ=255)

                                                                                    tACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTGGGaa     <  2:922670/93‑1 (MQ=255)

                                                                                    tACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTGGGaa     <  2:5378464/93‑1 (MQ=255)

                                                                                     aCGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTGGGAAg    <  2:3667364/93‑1 (MQ=255)

                                                                                     aCGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTGGGAAg    >  2:4089054/1‑93 (MQ=255)

                                                                                       ggggCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTGGGAAGgc  <  1:823721/93‑1 (MQ=255)

GATAACTCACAGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGACATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTGGGAAGGC  >  NC_000913/4162358‑4162537

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