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.6%	846.4 / ‑5.0	246	T411A (ACA→GCA)	btuB	BtuB
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (127/118); total (128/118)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 257 total aligned reads displayed.

GATAACTCACAGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGACATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTGGGAAGGC  >  NC_000913/4162358‑4162537

gATAACTCACAGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCa                                                                                         >  1:2218636/1‑93 (MQ=255)

gATAACTCACAGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCa                                                                                         >  1:4660436/1‑93 (MQ=255)

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

     cTCACAGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTta                                                                                    >  1:1771414/1‑93 (MQ=255)

     cTCACAGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTta                                                                                    <  1:1039546/93‑1 (MQ=255)

     cTCACAGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTta                                                                                    >  2:375168/1‑93 (MQ=255)

      tCACAGTTTGGTCGTCATGGAACATGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGCGCATCTtat                                                                                   >  2:4134433/1‑93 (MQ=255)

       cacaGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCCTCTtata                                                                                  >  1:4204791/1‑93 (MQ=255)

       cacaGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTtata                                                                                  >  1:4139242/1‑93 (MQ=255)

        acaGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATaa                                                                                 >  1:2037191/1‑93 (MQ=255)

        acaGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATaa                                                                                 >  2:8552990/1‑93 (MQ=255)

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

          aGTTTGGTCCTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAgg                                                                               <  2:1873844/93‑1 (MQ=255)

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

             ttGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCAc                                                                            >  1:3337523/1‑93 (MQ=255)

             ttGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCAc                                                                            >  1:7571331/1‑93 (MQ=255)

               tgtcgtcATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCa                                                                          <  2:5697002/92‑1 (MQ=255)

                 tcgtcATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCagt                                                                        <  2:7694255/93‑3 (MQ=255)

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

                   gtcATTGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATc                                                                      <  1:6689487/93‑1 (MQ=255)

                   gtcATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATc                                                                      >  2:4929840/1‑93 (MQ=255)

                    tcATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCt                                                                     <  2:2122057/93‑1 (MQ=255)

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

                       tGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTggg                                                                  <  1:360308/93‑1 (MQ=255)

                        ggAACCTGGCAAACCAGCGCCGGTTTTGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTgggg                                                                 <  2:8806887/93‑1 (MQ=255)

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

                          aaCCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCa                                                               >  1:6409849/1‑93 (MQ=255)

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

                           aCCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCaa                                                              <  1:5122076/93‑1 (MQ=255)

                            ccTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAAc                                                             <  2:3457330/93‑1 (MQ=255)

                             cTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACt                                                            >  1:5930162/1‑93 (MQ=255)

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

                               ggCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGt                                                          <  2:7189315/93‑1 (MQ=255)

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

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

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

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

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

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

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

                                    aCCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGc                                                     <  2:2515366/93‑1 (MQ=255)

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

                                     ccAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCt                                                    >  2:8523958/1‑93 (MQ=255)

                                     ccAGCGCCGGTTGGGAATTCATCGAAGGTCATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCt                                                    <  2:1618020/93‑1 (MQ=255)

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

                                       aGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGAAACTGTATGGCttt                                                  >  1:6258059/1‑92 (MQ=255)

                                        gcgcCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCt                                                 <  1:7079271/93‑1 (MQ=255)

                                                 gggAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATc                                        <  2:488340/93‑1 (MQ=255)

                                                 gggAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATc                                        <  1:5591591/93‑1 (MQ=255)

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

                                                    aaTTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGa                                     >  1:3953440/1‑93 (MQ=255)

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

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

                                                      ttCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAAt                                   >  2:2378874/1‑93 (MQ=255)

                                                      ttCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAAt                                   <  2:1249609/93‑1 (MQ=255)

                                                      ttCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAAt                                   <  1:4874420/93‑1 (MQ=255)

                                                       tCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATc                                  <  2:3947348/93‑1 (MQ=255)

                                                       tCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATc                                  <  1:8283141/93‑1 (MQ=255)

                                                       tCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATc                                  <  2:3711011/93‑1 (MQ=255)

                                                       tCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATc                                  >  1:4750185/1‑93 (MQ=255)

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

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

                                                         aTCGAAGGTTATCGCTTCATTGCTTTCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTAAGGAAATCACAATCTg                                >  2:2043013/1‑93 (MQ=255)

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

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

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

                                                          tCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTgg                               >  2:1641874/1‑93 (MQ=255)

                                                          acGAAGGGTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTgg                               <  1:2913014/92‑1 (MQ=255)

                                                           cGAAGGTTATCGCTTCATTGCTTCCTACGGGGGATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGa                              >  1:8138494/1‑93 (MQ=255)

                                                           cGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGa                              <  1:8117899/93‑1 (MQ=255)

                                                            gAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAAGCCGAATCTGGAc                             >  2:7039700/1‑93 (MQ=255)

                                                              aGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGAccc                           <  1:4452751/93‑1 (MQ=255)

                                                                gTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCgg                         <  2:7018420/93‑1 (MQ=255)

                                                                gTTATCGCTTCATTGCTTCCTACGGGACATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCgg                         >  1:8120604/1‑93 (MQ=255)

                                                                 ttATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGga                        >  2:5472308/1‑93 (MQ=255)

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

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

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

                                                                     cGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGaaa                    <  1:5079355/93‑1 (MQ=255)

                                                                      gCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGaaaa                   >  1:8574676/1‑93 (MQ=255)

                                                                       cTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAg                  <  2:7983802/93‑1 (MQ=255)

                                                                       cTTCATTGCTTCCTACGGAGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAg                  <  1:1454040/93‑1 (MQ=255)

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

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

                                                                           aTTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCaaa              <  2:6111092/93‑1 (MQ=255)

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

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

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

                                                                            ttGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACACGGAGAACAGAAAAc             >  2:2507210/1‑93 (MQ=255)

                                                                              gCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAg           <  2:5764211/93‑1 (MQ=255)

                                                                               cTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGt          <  1:7800874/93‑1 (MQ=255)

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

                                                                                 tCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTgg        <  2:3162400/93‑1 (MQ=255)

                                                                                 tCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTgg        >  2:5988805/1‑93 (MQ=255)

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

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

                                                                                     aCGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTGGGaat    >  2:364065/1‑92 (MQ=255)

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

                                                                                      cGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTGGGAAgg   >  2:678241/1‑93 (MQ=255)

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

GATAACTCACAGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGACATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTGGGAAGGC  >  NC_000913/4162358‑4162537

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