BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,613,759	A→G	L48P (CTG→CCG)	perM ←	putative permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,613,759	0	A	G	87.7%	574.0 / 62.9	212	L48P (CTG→CCG)	perM	putative permease
Reads supporting (aligned to +/- strand): ref base A (7/19); new base G (87/98); total (94/118)
Fisher's exact test for biased strand distribution p-value = 6.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 223 total aligned reads displayed.

GACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCG  >  NC_000913/2613668‑2613849

gacgaAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATTGGCCAg                                                                                           <  1:5523126/93‑1 (MQ=255)

  cgaAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGca                                                                                         <  2:439527/93‑1 (MQ=255)

  cgaAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCAGca                                                                                         >  1:4246176/1‑93 (MQ=255)

   gaAACCCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGcac                                                                                        <  1:2598402/93‑1 (MQ=255)

    aaaaCCACCAATACAATCGACGTCGCCTAGCGGCGGGGGCGACCAATAGATCGCAGACGCACGGTTGCGCATTCCAGCACATAGGCCGGcaca                                                                                       >  2:4246465/1‑93 (MQ=255)

       accaccAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATa                                                                                    <  2:5023192/93‑1 (MQ=255)

        ccaccaATACAATCGGCGTCGCCCAGCGGCGGGAGCGGCCAATAGATTGCAGAGGCACGGCTGGCCATTCCACCAAACAAGCCACcacaccaa                                                                                   >  2:3596335/1‑89 (MQ=255)

         caccaATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAAGCCGGCACAATAGc                                                                                  >  2:268778/1‑93 (MQ=255)

          accaATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGcc                                                                                 <  1:564212/93‑1 (MQ=255)

           ccaATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGcca                                                                                >  1:6321372/1‑93 (MQ=255)

            caATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGccac                                                                               >  2:1138286/1‑93 (MQ=255)

                   aTCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAAc                                                                        >  2:2924887/1‑93 (MQ=255)

                    tCGACGCCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACg                                                                       <  1:4050553/93‑1 (MQ=255)

                     cGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACgg                                                                      <  2:1103654/93‑1 (MQ=255)

                     cGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACgg                                                                      <  2:2395076/93‑1 (MQ=255)

                      gACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGa                                                                     >  1:592653/1‑93 (MQ=255)

                       aCGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGag                                                                    >  1:5276087/1‑93 (MQ=255)

                         gTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGagca                                                                  >  2:5373746/1‑93 (MQ=255)

                          tCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGagcaa                                                                 >  2:368896/1‑93 (MQ=255)

                           cGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGagcaag                                                                <  2:4835105/93‑1 (MQ=255)

                            gCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGagcaagc                                                               <  2:410218/93‑1 (MQ=255)

                             cccAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGagcaagca                                                              <  1:5120964/93‑1 (MQ=255)

                              ccAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAg                                                             <  1:3137134/93‑1 (MQ=255)

                              ccAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCCGCACAATAGCCACCAGCAACGGAGCAAGCAg                                                             <  1:5306443/93‑1 (MQ=255)

                              ccAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAAGCCGGCACAATAGCCACCAGCAACGGAGCACGCAg                                                             >  2:1008596/1‑93 (MQ=255)

                               cAGCGGCGGGGGCAGCCAATAGATTTAAGACGCACGGTTGACCATTCCAGCAAATAGGCAGACACAATAGCCACCAGCAACGGAGCAAGCAgg                                                            <  1:3072991/93‑1 (MQ=255)

                               cAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAgg                                                            >  2:6527621/1‑93 (MQ=255)

                                 gcggcgGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAAGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGcc                                                          >  2:6131657/1‑93 (MQ=255)

                                 gcggcgGGAGCAGCCAATAGATTGCAGACGCACGGGTGGCCATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGcc                                                          >  2:1417297/1‑93 (MQ=255)

                                    gcggGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGGCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACt                                                       >  1:2695910/1‑93 (MQ=255)

                                     cggGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTa                                                      <  2:2471516/93‑1 (MQ=255)

                                      gggAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTaa                                                     <  1:2089771/93‑1 (MQ=255)

                                       ggAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTaaa                                                    >  2:6337168/1‑93 (MQ=255)

                                        gAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAg                                                   <  2:1878280/93‑1 (MQ=255)

                                         agcagcCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGa                                                  >  1:728910/1‑93 (MQ=255)

                                         agcagcCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGa                                                  <  2:4885571/93‑1 (MQ=255)

                                          gcagcCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGaa                                                 <  2:1823494/93‑1 (MQ=255)

                                             gcCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCCGGCCACTAAAGacaaa                                              >  2:3861680/1‑93 (MQ=255)

                                             gcCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGaaaaa                                              >  2:1384931/1‑93 (MQ=255)

                                             gcCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGaaaaa                                              <  2:3529344/93‑1 (MQ=255)

                                             gcCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGaaaaa                                              >  1:3013508/1‑93 (MQ=255)

                                               cAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGa                                            <  2:2289930/93‑1 (MQ=255)

                                                 aTAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATa                                          <  2:5956237/93‑1 (MQ=255)

                                                   aGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAAt                                        <  2:3933656/93‑1 (MQ=255)

                                                   aGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAAt                                        >  2:1044842/1‑93 (MQ=255)

                                                     aTTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGc                                      <  1:2148336/93‑1 (MQ=255)

                                                     aTTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGc                                      >  1:4356554/1‑93 (MQ=255)

                                                      ttGCAGACGCCCGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGcc                                     <  1:1627/93‑1 (MQ=255)

                                                        gCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCaa                                   <  2:4636660/93‑1 (MQ=255)

                                                         cAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCaaa                                  >  1:952971/1‑93 (MQ=255)

                                                         cAGACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCaaa                                  <  1:1930138/93‑1 (MQ=255)

                                                           gACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATc                                >  2:4591548/1‑93 (MQ=255)

                                                           gACGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAACCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATc                                >  2:1919032/1‑93 (MQ=255)

                                                             cGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCg                              <  2:2722370/93‑1 (MQ=255)

                                                             cGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCg                              <  1:4092728/93‑1 (MQ=255)

                                                             cGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCg                              <  2:4230818/93‑1 (MQ=255)

                                                             cGCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCg                              <  2:4683609/93‑1 (MQ=255)

                                                              gCACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAGAAAGATAATGCCAAATCCgg                             >  2:3307862/1‑93 (MQ=255)

                                                               cACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGc                            >  2:2043525/1‑93 (MQ=255)

                                                               cACGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGc                            >  1:2820357/1‑93 (MQ=255)

                                                               cACGG‑TGGCCATCCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCa                           <  1:4409835/93‑1 (MQ=255)

                                                                aCGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCa                           >  2:6257454/1‑93 (MQ=255)

                                                                aCGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCa                           <  2:2751806/93‑1 (MQ=255)

                                                                aCGGTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCa                           >  2:3349525/1‑93 (MQ=255)

                                                                   gTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACt                        >  1:3558525/1‑93 (MQ=255)

                                                                   gTTGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACt                        >  1:5398952/1‑93 (MQ=255)

                                                                   gTTGGCCATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACt                        >  2:6064539/1‑93 (MQ=255)

                                                                    ttGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTa                       >  1:1685379/1‑93 (MQ=255)

                                                                    ttGGCCATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTa                       <  2:6160705/93‑1 (MQ=255)

                                                                     tGGCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTaa                      >  1:3982718/1‑93 (MQ=255)

                                                                      ggCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTaaa                     >  1:2410399/1‑93 (MQ=255)

                                                                      ggCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTaaa                     >  2:3020789/1‑93 (MQ=255)

                                                                       gCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCCACTaaaa                    >  1:3362586/1‑93 (MQ=255)

                                                                       gCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTaaaa                    <  2:837327/93‑1 (MQ=255)

                                                                       gCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTaaaa                    <  1:4451686/93‑1 (MQ=255)

                                                                       gCCATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTaaaa                    <  2:2129442/93‑1 (MQ=255)

                                                                        ccATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAAt                   >  2:4053947/1‑93 (MQ=255)

                                                                         cATTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATa                  >  2:2528898/1‑93 (MQ=255)

                                                                         cATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATa                  <  2:4086569/93‑1 (MQ=255)

                                                                          aTTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATaa                 >  2:209102/1‑93 (MQ=255)

                                                                          aTTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATaa                 >  1:6038706/1‑93 (MQ=255)

                                                                          aTTCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATaa                 >  1:645600/1‑93 (MQ=255)

                                                                           ttCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAAc                >  1:5231707/1‑93 (MQ=255)

                                                                           ttCCAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAAc                >  1:3526457/1‑93 (MQ=255)

                                                                           ttCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAAc                >  2:5177676/1‑93 (MQ=255)

                                                                              cAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAg             >  1:6275555/1‑93 (MQ=255)

                                                                              cAGCAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAg             >  1:1154222/1‑93 (MQ=255)

                                                                              cAGCAAATAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAg             <  1:2834119/93‑1 (MQ=255)

                                                                                 cAAATAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCaa          <  2:261939/93‑1 (MQ=255)

                                                                                   aaTAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAgg        >  1:2532430/1‑93 (MQ=255)

                                                                                   aaTAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAgg        >  1:1060227/1‑93 (MQ=255)

                                                                                   aaTAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAgg        >  1:5747304/1‑93 (MQ=255)

                                                                                   aaTAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAgg        <  2:123902/93‑1 (MQ=255)

                                                                                     tAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCa      >  2:1391209/1‑93 (MQ=255)

                                                                                     tAGGCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCa      >  2:1729290/1‑93 (MQ=255)

                                                                                       ggCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAAt    <  1:1263244/93‑1 (MQ=255)

                                                                                       ggCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAAt    <  1:4542643/93‑1 (MQ=255)

                                                                                        gCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATc   <  1:3595327/93‑1 (MQ=255)

                                                                                        gCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATc   <  1:2037614/93‑1 (MQ=255)

                                                                                        gCCGGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATc   <  2:2195810/93‑1 (MQ=255)

GACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCG  >  NC_000913/2613668‑2613849

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