BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,014,125	G→A	A36A (GCG→GCA)	udp →	uridine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,014,125	0	G	A	99.5%	768.6 / ‑5.9	216	A36A (GCG→GCA)	udp	uridine phosphorylase
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (99/116); total (100/116)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 222 total aligned reads displayed.

TTTTCATCTCGGCCTCACTAAAAACGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCGCTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTGGATGGTAAACCTGTTATCGTCTGCTC  >  NC_000913/4014035‑4014214

 tttCATCTCGGCCTCACTAAAAACGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTg                                                                                        >  2:3358449/1‑93 (MQ=255)

 tttCATCTCGGCCTCACTAAAAACGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTg                                                                                        >  2:4635637/1‑93 (MQ=255)

  ttCATCTCGGCCTCACTAAAAACGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGa                                                                                       <  1:3609158/93‑1 (MQ=255)

  ttCATCTCGGCCTCACTAAAAACGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGa                                                                                       <  2:1613258/93‑1 (MQ=255)

  ttCATCTCGGCCTCACTAAAAACGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGa                                                                                       <  2:1741650/93‑1 (MQ=255)

   tCATCTCGGCCTCACTAAAAACGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGAt                                                                                      <  2:754999/93‑1 (MQ=255)

   tCATCGCGGCCTCACTAAAAACGATTTACAAGGGGCTACGCTCGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGAt                                                                                      <  2:740817/93‑1 (MQ=255)

     aTCTCGGCCTCACTAAAAACGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATgg                                                                                    <  1:2908826/93‑1 (MQ=255)

     aTCTCGGCCTCACTAAAAACGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATgg                                                                                    >  2:1784734/1‑93 (MQ=255)

      tctcGGCCTCACTAAAAACGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGa                                                                                   >  1:2760895/1‑93 (MQ=255)

      tctcGGCCTCACTAAAAACGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGa                                                                                   <  1:2079065/93‑1 (MQ=255)

       ccccgCCCCACCAAAAACGATTTACCAGGGGCCACCCCTGCCAACGGCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGAt                                                                                  <  1:1178556/89‑1 (MQ=255)

        tcGGCCTCACTAAAAACGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATa                                                                                 >  1:3348103/1‑93 (MQ=255)

         cGGCCTCACTAAAAACGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATaa                                                                                <  1:3207402/93‑1 (MQ=255)

           gCCTCACTAAAAACGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGc                                                                              >  1:3166707/1‑93 (MQ=255)

            ccTCACTAAAAACGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGcc                                                                             <  1:4323025/93‑1 (MQ=255)

              tCACTAAAAACGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCgg                                                                           <  2:2538919/93‑1 (MQ=255)

               cACTAAAAACGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATGGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGt                                                                          <  2:3732665/93‑1 (MQ=255)

               cACTAAAAACGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGt                                                                          >  1:4151698/1‑93 (MQ=255)

               cACTAAAAACGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGt                                                                          <  1:1557045/93‑1 (MQ=255)

                aCTAAAAACGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGtt                                                                         <  1:3252147/93‑1 (MQ=255)

                aCTAAAAACGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGtt                                                                         <  2:1486457/93‑1 (MQ=255)

                aCTAAAAACGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGtt                                                                         >  1:3092694/1‑93 (MQ=255)

                  tAAAAACGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTaa                                                                       <  2:3025016/93‑1 (MQ=255)

                  tAAAAACGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTaa                                                                       <  1:953791/93‑1 (MQ=255)

                    aaaaCGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGc                                                                     <  2:2119498/93‑1 (MQ=255)

                    aaaaCGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGc                                                                     <  1:2545556/93‑1 (MQ=255)

                    aaaaCGATTTACAAGGGGCTACGCTGGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGc                                                                     <  2:1552838/93‑1 (MQ=255)

                     aaaCGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCt                                                                    <  1:1101756/93‑1 (MQ=255)

                       aCGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTgg                                                                  >  1:1546125/1‑93 (MQ=255)

                       aCGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTgg                                                                  <  2:4420475/93‑1 (MQ=255)

                         gATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGct                                                                >  2:3482905/1‑92 (MQ=255)

                         gATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCa                                                                >  1:2565013/1‑93 (MQ=255)

                           tttACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCAtc                                                              <  1:132590/93‑1 (MQ=255)

                           tttACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCAtc                                                              >  1:3916670/1‑93 (MQ=255)

                             tACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCAtctc                                                            >  1:1722853/1‑93 (MQ=255)

                              aCAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCa                                                           <  2:3707/93‑1 (MQ=255)

                              aCAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCa                                                           <  2:1861448/93‑1 (MQ=255)

                                aaGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCAcc                                                         <  2:2439666/93‑1 (MQ=255)

                                 aGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACcg                                                        >  2:260922/1‑93 (MQ=255)

                                 aGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGACCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAACTGGCATCTCACcg                                                        >  1:2493143/1‑93 (MQ=255)

                                   gggCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGGTAAGCTGGCATCTCACcgcg                                                      >  2:3639545/1‑93 (MQ=255)

                                       tACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAAtt                                                  <  1:4084334/93‑1 (MQ=255)

                                       tACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAAtt                                                  <  2:4016184/93‑1 (MQ=255)

                                       tACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAAtt                                                  >  2:4246015/1‑93 (MQ=255)

                                        aCGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTc                                                 >  2:3872697/1‑93 (MQ=255)

                                         cGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCa                                                >  1:1795551/1‑93 (MQ=255)

                                          gCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCAc                                               >  1:4439096/1‑93 (MQ=255)

                                            ttGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTa                                             >  2:4669719/1‑93 (MQ=255)

                                            ttGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACAGCGAATTCACTa                                             >  1:4210053/1‑93 (MQ=255)

                                             tGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTAc                                            <  2:555979/93‑1 (MQ=255)

                                             tGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTAc                                            <  2:2225289/93‑1 (MQ=255)

                                               ccATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCt                                          <  2:1044908/93‑1 (MQ=255)

                                                cATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTg                                         >  2:188976/1‑93 (MQ=255)

                                                 aTCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTgg                                        <  1:511623/93‑1 (MQ=255)

                                                 aTCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTgg                                        >  2:2164871/1‑93 (MQ=255)

                                                 aTCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTgg                                        >  2:2527233/1‑93 (MQ=255)

                                                  tCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGc                                       <  2:3099530/93‑1 (MQ=255)

                                                  tCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGc                                       <  2:1731907/93‑1 (MQ=255)

                                                    gTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGt                                     <  2:4381955/93‑1 (MQ=255)

                                                      cccTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGc                                   >  2:2946093/1‑93 (MQ=255)

                                                       ccTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCa                                  <  1:2687997/93‑1 (MQ=255)

                                                       ccTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCa                                  <  1:4544093/93‑1 (MQ=255)

                                                        cTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCag                                 >  2:229865/1‑93 (MQ=255)

                                                        cTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCag                                 >  1:2626778/1‑93 (MQ=255)

                                                         tGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCaga                                <  2:4197107/93‑1 (MQ=255)

                                                         tGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCaga                                <  1:1078560/93‑1 (MQ=255)

                                                         tGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCaga                                <  1:4641982/93‑1 (MQ=255)

                                                         tGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCaga                                >  1:694240/1‑93 (MQ=255)

                                                         tGGCGACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAAATCACTACCTGGCGTGCaga                                >  2:3212916/1‑93 (MQ=255)

                                                             gACCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTg                            <  2:492693/93‑1 (MQ=255)

                                                              aCCCGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTgg                           >  1:3473062/1‑93 (MQ=255)

                                                               cccGGATCGTGTGTGAAAGATCGCCGGACTGATGGATAAGCCTGGTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTGGa                          <  1:851584/93‑1 (MQ=255)

                                                               cccGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTGGa                          <  2:409653/93‑1 (MQ=255)

                                                                ccGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTGGAt                         >  1:1045371/1‑93 (MQ=255)

                                                                ccGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTGGAt                         <  2:4132019/93‑1 (MQ=255)

                                                                 cGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGTGCTGGATg                        <  2:4574260/93‑1 (MQ=255)

                                                                 cGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTGGATg                        >  2:265228/1‑93 (MQ=255)

                                                                 cGGATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTGGATg                        >  1:369463/1‑93 (MQ=255)

                                                                  ggATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTGGATgg                       >  2:1469260/1‑93 (MQ=255)

                                                                  ggATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTGGATgg                       >  2:2466045/1‑93 (MQ=255)

                                                                  ggATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTGGATgg                       >  2:2876290/1‑93 (MQ=255)

                                                                   gATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTGGATGGt                      >  2:4118811/1‑93 (MQ=255)

                                                                   gATCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGATGGATGGt                      >  2:857823/1‑93 (MQ=255)

                                                                     tCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTGGATGGTaa                    <  1:1950589/93‑1 (MQ=255)

                                                                     tCGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTGGATGGTaa                    <  2:2795460/93‑1 (MQ=255)

                                                                      cGTGTGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTGGATGGTaaa                   <  1:2652535/93‑1 (MQ=255)

                                                                       gtgtGGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTGGATGGTAAAc                  >  1:1694552/1‑93 (MQ=255)

                                                                        tgtgGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTGGATGGTAAAcc                 <  2:1987585/93‑1 (MQ=255)

                                                                        tgtgGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTGGATGGTAAAcc                 >  1:2024607/1‑93 (MQ=255)

                                                                        tgtgGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTGGATGGTAAAcc                 >  2:2703583/1‑93 (MQ=255)

                                                                          tgGAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTGGATGGTAAACCTg               <  1:3329271/93‑1 (MQ=255)

                                                                           ggAAAAGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTGGATGGTAAACCTGt              <  2:1317062/93‑1 (MQ=255)

                                                                            gAAAAGATCGCAGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTGGATGGTAAACCTGtt             <  2:3557307/93‑1 (MQ=255)

                                                                             aaaaGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTGGATGGTAAACCTGTTa            <  1:2276531/93‑1 (MQ=255)

                                                                             aaaaGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTGGATGGTAAACCTGTTa            >  2:2342576/1‑93 (MQ=255)

                                                                              aaaGATCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTGGATGGTAAACCTGTTAt           <  2:395900/93‑1 (MQ=255)

                                                                                  aTCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTGGATGGTAAACCTGTTATCGTc       <  2:614907/93‑1 (MQ=255)

                                                                                   tCGCCGCACTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTGGATGGTAAACCTGTTATCGTCt      <  1:3461349/93‑1 (MQ=255)

                                                                                       cgcaCTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTGGATGGTAAACCTGTTATCGTCTGctc  <  1:3947077/93‑1 (MQ=255)

TTTTCATCTCGGCCTCACTAAAAACGATTTACAAGGGGCTACGCTTGCCATCGTCCCTGGCGACCCGGATCGTGTGGAAAAGATCGCCGCGCTGATGGATAAGCCGGTTAAGCTGGCATCTCACCGCGAATTCACTACCTGGCGTGCAGAGCTGGATGGTAAACCTGTTATCGTCTGCTC  >  NC_000913/4014035‑4014214

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