BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,561,348	A→G	C49C (TGT→TGC)	b4563 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,561,348	0	A	G	99.5%	608.5 / ‑3.8	194	coding (147/258 nt)	b4563	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (92/101); minor base . (0/1); total (92/102)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-01

Only 100 of 201 total aligned reads displayed.

ATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGC  >  NC_000913/4561258‑4561437

aTGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCa                                                                                         >  1:3453006/1‑93 (MQ=255)

 tgtgCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAg                                                                                        <  1:6278673/93‑1 (MQ=255)

 tgtgCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAg                                                                                        >  1:2054702/1‑93 (MQ=255)

  gtgCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAgg                                                                                       >  2:111323/1‑93 (MQ=255)

   tgCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGt                                                                                      <  1:4806732/93‑1 (MQ=255)

    gCAATTCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTc                                                                                     <  2:1941685/93‑1 (MQ=255)

     cAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCa                                                                                    <  1:1539698/93‑1 (MQ=255)

     cAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCa                                                                                    <  1:1608516/93‑1 (MQ=255)

     cAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCa                                                                                    <  1:3036156/93‑1 (MQ=255)

       aTGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGc                                                                                  >  1:3279144/1‑93 (MQ=255)

       aTGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGc                                                                                  >  2:5274818/1‑93 (MQ=255)

        tGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGcc                                                                                 >  2:6130317/1‑93 (MQ=255)

        tGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGcc                                                                                 >  2:5457515/1‑93 (MQ=255)

           aGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCaaa                                                                              <  1:3565555/93‑1 (MQ=255)

            gcgcCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAAt                                                                             >  1:5360778/1‑93 (MQ=255)

               ccAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCa                                                                          >  1:4712725/1‑93 (MQ=255)

                  gACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATaa                                                                       <  2:4402347/93‑1 (MQ=255)

                   acacCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAAc                                                                      <  1:4595486/93‑1 (MQ=255)

                    cacCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAAcc                                                                     >  1:6126470/1‑93 (MQ=255)

                    cacCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAAcc                                                                     >  1:3902649/1‑93 (MQ=255)

                     acCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGTATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAAccc                                                                    <  2:4064471/93‑1 (MQ=255)

                     acCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAAccc                                                                    >  2:1424440/1‑93 (MQ=255)

                     acCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAAccc                                                                    <  1:4635230/93‑1 (MQ=255)

                     acCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAAccc                                                                    >  2:253781/1‑93 (MQ=255)

                     acCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAAccc                                                                    <  1:5715719/93‑1 (MQ=255)

                       cATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGc                                                                  <  1:4949973/93‑1 (MQ=255)

                       cATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGc                                                                  <  1:402296/93‑1 (MQ=255)

                        aTCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCa                                                                 <  1:179684/93‑1 (MQ=255)

                        aTCCTGCTGGCAAAAAATTGCCTCAGTCAGCCCAGGGATCCCCCATTGTATTTCCCCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCa                                                                 <  1:284239/93‑1 (MQ=255)

                             gctgGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAgccgc                                                            >  2:2145723/1‑93 (MQ=255)

                              ctgGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAgccgcc                                                           >  1:6205233/1‑93 (MQ=255)

                               tgGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGccgccg                                                          >  1:845663/1‑93 (MQ=255)

                               tgGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGccgccg                                                          <  1:1863360/93‑1 (MQ=255)

                                 gCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCcgcg                                                        >  2:2046011/1‑93 (MQ=255)

                                   aaTAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTa                                                      <  2:5740890/93‑1 (MQ=255)

                                    aTAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTaa                                                     >  1:1898641/1‑93 (MQ=255)

                                    aTAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTaa                                                     >  1:4530558/1‑93 (MQ=255)

                                     tAAATTGACTCAGTCAGCCCAGGGATCTCCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTaaa                                                    <  1:2708762/93‑1 (MQ=255)

                                     tAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTaaa                                                    <  2:5192122/93‑1 (MQ=255)

                                     tAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTaaa                                                    <  2:5970663/93‑1 (MQ=255)

                                        aTTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCa                                                 >  2:3587766/1‑93 (MQ=255)

                                        aTTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCa                                                 >  2:4487365/1‑93 (MQ=255)

                                          tGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGG‑CAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACgg                                              <  2:626806/93‑1 (MQ=255)

                                           gACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACgg                                              <  1:3779866/93‑1 (MQ=255)

                                           gACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACgg                                              >  1:250122/1‑93 (MQ=255)

                                               cagtcagCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGaa                                          <  2:5005137/93‑1 (MQ=255)

                                               cagtcagCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCAAGGGGaa                                          >  2:4301161/1‑93 (MQ=255)

                                               cagtcagCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCACAGCCGCCGCGTAAATCACGGGGaa                                          <  2:3142237/93‑1 (MQ=255)

                                                 gtcagCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAata                                        <  1:713248/93‑1 (MQ=255)

                                                  tcagCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAatat                                       >  2:1712296/1‑93 (MQ=255)

                                                  tcagCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAatat                                       <  2:2068961/93‑1 (MQ=255)

                                                  tcagCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAatat                                       <  1:5598173/93‑1 (MQ=255)

                                                   cagCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATc                                      <  2:785003/93‑1 (MQ=255)

                                                    agCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCg                                     <  1:972203/93‑1 (MQ=255)

                                                    agCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCg                                     >  2:201754/1‑93 (MQ=255)

                                                     gCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGc                                    <  1:2347400/93‑1 (MQ=255)

                                                      cccAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCt                                   <  2:2183198/93‑1 (MQ=255)

                                                      cccAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCt                                   <  1:4839244/93‑1 (MQ=255)

                                                      cccAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCt                                   <  2:6084859/93‑1 (MQ=255)

                                                        cAGGGATCCCCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGa                                 <  1:5272239/93‑1 (MQ=255)

                                                        cAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATACCCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGa                                 <  2:3676008/93‑1 (MQ=255)

                                                        cAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGa                                 <  2:189255/93‑1 (MQ=255)

                                                        cAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGa                                 <  2:4855978/93‑1 (MQ=255)

                                                         aGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGAc                                <  1:4409142/93‑1 (MQ=255)

                                                            gATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGGAGCCGCCGCGTAAATCACGGGGAATATCGCTGAcggg                             >  2:2321573/1‑91 (MQ=255)

                                                              tCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATCACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGaa                           <  2:5358835/93‑1 (MQ=255)

                                                              tCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGaa                           >  1:412831/1‑93 (MQ=255)

                                                               cACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGaaa                          >  1:2979855/1‑93 (MQ=255)

                                                                ccccATTGTATTTCCGCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGaaaa                         <  2:1802415/92‑1 (MQ=255)

                                                                aCCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGaaaa                         <  2:5581263/93‑1 (MQ=255)

                                                                aCCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGaaaa                         >  1:3310366/1‑93 (MQ=255)

                                                                 cccATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAg                        >  1:5226785/1‑93 (MQ=255)

                                                                  ccATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATTACCCGCAGCCGCCGCGTAAATCACGGGGAATATTGATGACGCGAACAgc                       >  2:5932977/1‑92 (MQ=255)

                                                                   cATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTc                      <  2:1358205/93‑1 (MQ=255)

                                                                   cATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTc                      <  1:1196462/93‑1 (MQ=255)

                                                                   cATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTc                      >  1:1003351/1‑93 (MQ=255)

                                                                   cATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTc                      <  1:5438338/93‑1 (MQ=255)

                                                                    aTTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCa                     <  1:381223/93‑1 (MQ=255)

                                                                    aTTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCa                     <  2:3526168/93‑1 (MQ=255)

                                                                    aTTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCa                     <  2:1497502/93‑1 (MQ=255)

                                                                    aTTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCa                     >  1:631719/1‑93 (MQ=255)

                                                                     ttGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTGAt                    >  1:3817352/1‑93 (MQ=255)

                                                                     ttGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCAt                    >  1:3139165/1‑93 (MQ=255)

                                                                     ttGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAAACCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCAt                    >  1:4432670/1‑93 (MQ=255)

                                                                      tGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATg                   <  2:4137405/93‑1 (MQ=255)

                                                                      tGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATg                   <  1:1932033/93‑1 (MQ=255)

                                                                      tGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATg                   <  1:2759690/93‑1 (MQ=255)

                                                                       gcaTTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGc                  >  2:2851756/3‑93 (MQ=255)

                                                                          tttCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCg               >  2:4436638/1‑93 (MQ=255)

                                                                           ttCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCACCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGc              >  1:3383894/1‑93 (MQ=255)

                                                                             ccTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGctg            <  2:46088/93‑1 (MQ=255)

                                                                             ccTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGctg            <  1:5549699/93‑1 (MQ=255)

                                                                              cTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGAAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATCCCCGctgc           >  1:3955617/1‑93 (MQ=255)

                                                                               tCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGctgcc          >  2:3050086/1‑93 (MQ=255)

                                                                                 gCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGctgcctg        >  1:6026411/1‑93 (MQ=255)

                                                                                   cGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCa      >  1:1386180/1‑93 (MQ=255)

                                                                                    ggggggGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCaa     >  2:2943173/1‑93 (MQ=255)

                                                                                     ggggggCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAAc    >  2:3506671/1‑93 (MQ=255)

                                                                                     ggggggCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAAc    <  1:5975140/93‑1 (MQ=255)

                                                                                      ggggGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACg   >  1:2185523/1‑93 (MQ=255)

ATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGC  >  NC_000913/4561258‑4561437

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