BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,220,652	A→G	Y199H (TAT→CAT)	iclR ←	IclR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,220,652	0	A	G	87.5%	810.6 / 97.1	321	Y199H (TAT→CAT)	iclR	IclR
Reads supporting (aligned to +/- strand): ref base A (18/22); new base G (145/136); total (163/158)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.84e-01

Only 100 of 330 total aligned reads displayed.

TGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCCGCACCGGAAGCGTGCA  >  NC_000913/4220561‑4220741

 gTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATggg                                                                                         >  1:6961573/1‑93 (MQ=255)

   tCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCa                                                                                       >  1:1837723/1‑93 (MQ=255)

    ccTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCAt                                                                                      >  1:85797/1‑93 (MQ=255)

     cTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATg                                                                                     >  2:5740327/1‑93 (MQ=255)

      tCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGt                                                                                    <  1:2855505/93‑1 (MQ=255)

       cATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTa                                                                                   <  2:3301733/93‑1 (MQ=255)

        aTCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTaa                                                                                  >  1:5154967/1‑93 (MQ=255)

         tCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAAc                                                                                 >  1:3999285/1‑93 (MQ=255)

         tCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAAc                                                                                 >  1:2425977/1‑93 (MQ=255)

          cGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAAcc                                                                                <  1:4980194/93‑1 (MQ=255)

           gTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAAccc                                                                               <  1:4536750/93‑1 (MQ=255)

           gTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAAccc                                                                               <  1:5353047/93‑1 (MQ=255)

           gTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGAATGGGCA‑GGAACCCt                                                                              >  1:6269317/1‑93 (MQ=255)

            tCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGCAACCCt                                                                              <  2:7094141/93‑1 (MQ=255)

              aaaTGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCttt                                                                            <  1:6305648/93‑1 (MQ=255)

               aaTGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTg                                                                           >  1:2727217/1‑93 (MQ=255)

               aaTGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTg                                                                           >  2:735233/1‑93 (MQ=255)

               aaTGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTCTg                                                                           <  2:3087314/93‑1 (MQ=255)

                aTGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGc                                                                          >  1:1091503/1‑93 (MQ=255)

                 tGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCg                                                                         >  2:4381191/1‑93 (MQ=255)

                  gAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCgg                                                                        >  2:2027752/1‑93 (MQ=255)

                   aaTAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGt                                                                       <  1:2632706/93‑1 (MQ=255)

                   aaTAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGt                                                                       >  2:7173061/1‑93 (MQ=255)

                     tAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTgc                                                                     <  1:366752/93‑1 (MQ=255)

                      aaCCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTgca                                                                    >  2:699327/1‑93 (MQ=255)

                       aCCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTgcag                                                                   >  1:4172714/1‑93 (MQ=255)

                        ccccGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTCTGGGCATGTAACCCTTTGCGGTgcagc                                                                  <  2:5164010/93‑1 (MQ=255)

                         cccGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTgcagca                                                                 <  1:5240902/93‑1 (MQ=255)

                         cccGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTgcagca                                                                 >  2:1135220/1‑93 (MQ=255)

                          ccGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGcagcag                                                                >  2:4105000/1‑93 (MQ=255)

                          ccGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGcagcag                                                                >  1:1329961/1‑93 (MQ=255)

                           cGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCagcagc                                                               <  2:6184940/93‑1 (MQ=255)

                            gTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCt                                                              >  1:2326375/1‑93 (MQ=255)

                            gTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCt                                                              <  1:5863777/93‑1 (MQ=255)

                             tttGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCtt                                                             >  2:4228279/1‑93 (MQ=255)

                             tttGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCtt                                                             <  1:7306850/93‑1 (MQ=255)

                             tttGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCtt                                                             >  1:4576018/1‑93 (MQ=255)

                              ttGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTc                                                            <  1:5647172/93‑1 (MQ=255)

                              ttGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTc                                                            >  2:2831690/1‑93 (MQ=255)

                              ttGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTc                                                            <  2:3857020/93‑1 (MQ=255)

                               tGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCg                                                           <  2:1433057/93‑1 (MQ=255)

                                    gTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCAcc                                                      >  1:3757181/1‑93 (MQ=255)

                                    gTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCAcc                                                      >  2:5590503/1‑93 (MQ=255)

                                        gggCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGtt                                                  <  1:3916589/93‑1 (MQ=255)

                                         ggCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGttc                                                 >  1:1879322/1‑93 (MQ=255)

                                         ggCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGttc                                                 >  2:3381705/1‑93 (MQ=255)

                                          gCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGttct                                                <  1:5207897/93‑1 (MQ=255)

                                          gCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGttct                                                <  1:3219440/93‑1 (MQ=255)

                                           cGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGttctt                                               <  2:3909192/93‑1 (MQ=255)

                                            gagaTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGttcttc                                              >  2:4916287/1‑93 (MQ=255)

                                            gagaTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGttcttc                                              >  1:7094407/1‑93 (MQ=255)

                                            gagaTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGttcttc                                              >  1:2522373/1‑93 (MQ=255)

                                            gagaTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGttcttc                                              >  2:3672343/1‑93 (MQ=255)

                                            gagaTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGttcttc                                              <  2:6391483/93‑1 (MQ=255)

                                              gaTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGc                                            >  1:6180916/1‑93 (MQ=255)

                                               aTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCt                                           >  1:3042025/1‑93 (MQ=255)

                                               aTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCt                                           >  2:6226131/1‑93 (MQ=255)

                                                tcttctTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTc                                          >  2:5623294/1‑93 (MQ=255)

                                                 cttcttTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCa                                         >  1:7456860/1‑93 (MQ=255)

                                                 cttcttTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCa                                         >  1:2031753/1‑93 (MQ=255)

                                                   tcttTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGt                                       >  2:4654708/1‑93 (MQ=255)

                                                   tcttTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGt                                       >  2:4554817/1‑93 (MQ=255)

                                                   tcttTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGt                                       <  1:5441452/93‑1 (MQ=255)

                                                    cttTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGtt                                      >  2:5002917/1‑93 (MQ=255)

                                                    cttTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGtt                                      <  2:2373769/93‑1 (MQ=255)

                                                    cttTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGtt                                      <  1:2159515/93‑1 (MQ=255)

                                                     ttttAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTg                                     <  2:2277900/93‑1 (MQ=255)

                                                      tttAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTgg                                    >  1:6046384/1‑93 (MQ=255)

                                                      tttAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTgg                                    <  1:1514218/93‑1 (MQ=255)

                                                         aaaTGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCt                                 >  1:3553115/1‑93 (MQ=255)

                                                            tGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTaaa                              <  2:1655103/93‑1 (MQ=255)

                                                             gCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTaaaa                             <  2:5869967/93‑1 (MQ=255)

                                                             gCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTaaaa                             <  2:5660051/93‑1 (MQ=255)

                                                              cacaGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTaaaaa                            <  2:4197224/93‑1 (MQ=255)

                                                               acaGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAg                           >  1:1036886/1‑93 (MQ=255)

                                                                caGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAgg                          >  2:7267195/1‑93 (MQ=255)

                                                                 aGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGc                         <  1:4956570/93‑1 (MQ=255)

                                                                  ggAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCt                        <  2:6482126/93‑1 (MQ=255)

                                                                   gagaCACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCtt                       <  1:4312000/93‑1 (MQ=255)

                                                                   gagaCACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCtt                       >  1:4735327/1‑93 (MQ=255)

                                                                   gagaCACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCtt                       >  2:3545339/1‑93 (MQ=255)

                                                                    agaCACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCttt                      <  1:5974052/93‑1 (MQ=255)

                                                                    agaCACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCttt                      <  1:6744240/93‑1 (MQ=255)

                                                                    agaCACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCttt                      <  1:1443806/93‑1 (MQ=255)

                                                                    agaCACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCttt                      <  1:3463406/93‑1 (MQ=255)

                                                                     gaCACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTa                     <  2:1276202/93‑1 (MQ=255)

                                                                     gaCACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTa                     <  2:3215780/93‑1 (MQ=255)

                                                                         ccAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCCg                 >  1:5525390/1‑93 (MQ=255)

                                                                         ccAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCCg                 >  2:3112255/1‑93 (MQ=255)

                                                                           aGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCCGCa               >  2:6721313/1‑93 (MQ=255)

                                                                            gCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCCGCAc              <  2:7516622/93‑1 (MQ=255)

                                                                              gTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCCGCACCg            >  2:6549756/1‑93 (MQ=255)

                                                                              gTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCCGCACCg            >  2:5337959/1‑93 (MQ=255)

                                                                               ttGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCCGCACCgg           >  1:5562015/1‑93 (MQ=255)

                                                                                  cGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCCGCACCGGAAg        <  1:5146123/93‑1 (MQ=255)

                                                                                   gTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCCGCACCGGAAGc       >  2:6800571/1‑93 (MQ=255)

                                                                                       gTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCCGCACCGGAAGCGTGc   >  1:2084016/1‑93 (MQ=255)

                                                                                       gTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCCGCACCGGAAGCGTGc   >  2:496674/1‑93 (MQ=255)

                                                                                        tatgGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCCGCACCGGAAGCGTGCa  >  2:2988617/1‑93 (MQ=255)

                                                                                        tataGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCCGCACCGGAAGCGTGCa  >  1:4240250/1‑93 (MQ=255)

TGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCCGCACCGGAAGCGTGCA  >  NC_000913/4220561‑4220741

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