BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,175,787	G→A	7.3%	noncoding (34/76 nt)	thrT →	tRNA‑Thr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,175,787	0	G	A	7.3%	582.6 / 23.0	191	noncoding (34/76 nt)	thrT	tRNA‑Thr
Reads supporting (aligned to +/- strand): ref base G (80/96); new base A (8/6); total (88/103)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.06e-01

Only 100 of 194 total aligned reads displayed.

GCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGG  >  NC_000913/4175639‑4175935

gCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTgg                                                                                                                                                      <  2:2071363/150‑1 (MQ=255)

gCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTgg                                                                                                                                                      <  2:1430401/150‑1 (MQ=255)

  aGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTa                                                                                                                                                    <  1:2277534/150‑1 (MQ=255)

  aGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTa                                                                                                                                                    <  1:1947393/150‑1 (MQ=255)

   gAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTaa                                                                                                                                                   >  2:185854/1‑150 (MQ=255)

    aaTGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAg                                                                                                                                                  >  1:170700/1‑150 (MQ=255)

    aaTGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTACAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAg                                                                                                                                                  >  2:753763/1‑150 (MQ=255)

      tGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAggg                                                                                                                                                <  1:119211/150‑1 (MQ=255)

      tGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAggg                                                                                                                                                <  2:435965/150‑1 (MQ=255)

      tGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAggg                                                                                                                                                <  2:838675/150‑1 (MQ=255)

       gAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGt                                                                                                                                               <  1:1379683/150‑1 (MQ=255)

       gAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGt                                                                                                                                               <  1:693058/150‑1 (MQ=255)

        aaaaTCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTACCAAGGGTg                                                                                                                                              >  2:1951259/1‑150 (MQ=255)

         aaaTCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGa                                                                                                                                             <  1:23273/150‑1 (MQ=255)

         aaaTCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGa                                                                                                                                             <  1:1902182/150‑1 (MQ=255)

          aaTCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTACCAAGGGTgcg                                                                                                                                            <  2:2202415/150‑3 (MQ=255)

           aTCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAgg                                                                                                                                           <  2:2490660/150‑1 (MQ=255)

            tCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTACCAAGGGTgcgct                                                                                                                                          <  2:57363/150‑5 (MQ=255)

            tCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTACCAAGGGTgcgct                                                                                                                                          <  2:905611/150‑5 (MQ=255)

               ggTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCgg                                                                                                                                       >  2:2168654/1‑150 (MQ=255)

                gTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGc                                                                                                                                      <  1:1795047/150‑1 (MQ=255)

                 tAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCa                                                                                                                                     >  1:212336/1‑150 (MQ=255)

                    ccGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑tt                                                                                                                                 >  2:324878/1‑150 (MQ=255)

                         ttCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAAt                                                                                                                            <  1:164081/150‑1 (MQ=255)

                         ttCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAAt                                                                                                                            >  1:764921/1‑150 (MQ=255)

                          tCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATc                                                                                                                           >  1:514611/1‑150 (MQ=255)

                           ccAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCt                                                                                                                          >  1:614600/1‑150 (MQ=255)

                            cAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTg                                                                                                                         <  1:1624557/150‑1 (MQ=255)

                             aGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGc                                                                                                                        <  1:1544174/150‑1 (MQ=255)

                             aGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGc                                                                                                                        >  2:933596/1‑150 (MQ=255)

                              ggATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGcc                                                                                                                       <  2:1657468/150‑1 (MQ=255)

                              ggATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGcc                                                                                                                       >  2:1271409/1‑150 (MQ=255)

                                 tGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGATTCGAA‑CTGccgac                                                                                                                    <  2:2018079/150‑4 (MQ=255)

                                   cGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCa                                                                                                                  >  1:2478618/1‑150 (MQ=255)

                                     ggCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGc                                                                                                                <  2:2292425/150‑1 (MQ=255)

                                           gTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCAct                                                                                                          >  2:669010/1‑150 (MQ=255)

                                           gTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCAct                                                                                                          >  1:1755064/1‑150 (MQ=255)

                                              taATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCTTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCActtct                                                                                                       >  2:1155943/1‑150 (MQ=255)

                                               aaTGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCActtctt                                                                                                      >  2:97264/1‑150 (MQ=255)

                                                aTGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCttt                                                                                                     <  2:2306026/150‑1 (MQ=255)

                                                 tGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCtttt                                                                                                    >  2:155339/1‑150 (MQ=255)

                                                   gCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTtct                                                                                                  <  1:595525/150‑1 (MQ=255)

                                                     tATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTctcc                                                                                                <  2:2308130/150‑1 (MQ=255)

                                                          ccTCAGACTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCt                                                                                           >  1:2362900/1‑150 (MQ=255)

                                                            tCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGt                                                                                         >  2:1000497/1‑150 (MQ=255)

                                                            tCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGt                                                                                         >  1:1781562/1‑150 (MQ=255)

                                                            tCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGt                                                                                         <  2:1266647/150‑1 (MQ=255)

                                                                ccTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGttttt                                                                                     <  1:1917084/150‑1 (MQ=255)

                                                                      aaGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTc                                                                               <  1:1687549/150‑1 (MQ=255)

                                                                          tgatgaTGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGtt                                                                           <  1:1577876/150‑1 (MQ=255)

                                                                          tgatgaTGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGtt                                                                           <  1:1665366/150‑1 (MQ=255)

                                                                          tgatgaTGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGtt                                                                           >  2:1932190/1‑150 (MQ=255)

                                                                             tgatgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTAt                                                                        >  2:832894/1‑150 (MQ=255)

                                                                              gatgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTAtt                                                                       <  1:1201593/150‑1 (MQ=255)

                                                                               ctgCGGGTTCGAGTCCCGCTGCCCGCTCCACGAGGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTg                                                                      <  1:1086276/149‑1 (MQ=255)

                                                                               atgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTg                                                                      >  1:1841278/1‑150 (MQ=255)

                                                                               atgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTg                                                                      <  2:1679010/150‑1 (MQ=255)

                                                                                tgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGc                                                                     <  1:502325/150‑1 (MQ=255)

                                                                                     gTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTca                                                                <  1:2246729/150‑1 (MQ=255)

                                                                                     gTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTca                                                                >  1:373504/1‑150 (MQ=255)

                                                                                       tCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTcaac                                                              <  1:2173380/150‑1 (MQ=255)

                                                                                        cGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGCGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTcaaca                                                             >  1:2245171/1‑150 (MQ=255)

                                                                                        cGATTCCCGCTGCCCGCTCCAAGATGAGCTTATATAGCTCAGTTGGTAGAGAGCACCCATGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCAGTTCTTTTCTCCTCCCTGTTTTATGCTTCTGTTTATTGCATTcaaca                                                             >  1:2059768/1‑150 (MQ=255)

                                                                                         gATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTcaacaa                                                            >  2:1891505/1‑150 (MQ=255)

                                                                                         gATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTcaacaa                                                            >  2:620973/1‑150 (MQ=255)

                                                                                          aTTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAg                                                           >  1:1956467/1‑150 (MQ=255)

                                                                                            tCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTc                                                         <  1:1648467/150‑1 (MQ=255)

                                                                                             cccGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCg                                                        <  1:933596/150‑1 (MQ=255)

                                                                                              ccGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCgg                                                       >  2:611472/1‑150 (MQ=255)

                                                                                                gCTGCGCGCTGCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTGCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGc                                                     <  2:1087813/150‑1 (MQ=255)

                                                                                                 cTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCTTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCa                                                    <  1:1155943/150‑1 (MQ=255)

                                                                                                   gCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATg                                                  >  1:2342450/1‑150 (MQ=255)

                                                                                                        ctgcAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGTCATGTTGcc                                             <  1:1440973/147‑1 (MQ=255)

                                                                                                        cTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGcc                                             <  1:1019112/150‑1 (MQ=255)

                                                                                                         tCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTTTTTATTGCATTCAACAAGTCGGGCATGTTGCCt                                            <  1:97264/150‑1 (MQ=255)

                                                                                                           cAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGTCAG‑TTCGAATCTGCCTATCAACACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATACTACAAGTCGGGCATGTTGCCTgg                                          >  2:2059154/1‑150 (MQ=255)

                                                                                                           cAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTgg                                          <  2:2342450/150‑1 (MQ=255)

                                                                                                              gATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTg                                       <  1:46613/150‑1 (MQ=255)

                                                                                                                 ttgCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCTGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGAtg                                    <  1:421754/149‑1 (MQ=255)

                                                                                                                     ggATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGGCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATgtggt                                <  1:736027/149‑1 (MQ=255)

                                                                                                                        tataGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGat                             <  1:629725/150‑1 (MQ=255)

                                                                                                                         ataGCTCTGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGata                            <  2:237453/150‑1 (MQ=255)

                                                                                                                           aGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATc                          <  1:287878/150‑1 (MQ=255)

                                                                                                                            gCTCAGTTGGTAGAGCGCACACTTGGTAAGGGTTAGGTAGGCAG‑TTCGAATCTGCATATCAGCACAACTTATTTTATCCTCCCTGTTTTTTCGTTCTGTTTATTGCATTCAACAAGTCGGGCCTGTTGCCTGGTTGATGTGGTGATATCa                         >  2:2384209/1‑150 (MQ=38)

                                                                                                                              tCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCAcc                       >  1:1659694/1‑150 (MQ=255)

                                                                                                                               cAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCg                      <  1:2318091/150‑1 (MQ=255)

                                                                                                                               cAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCg                      <  2:514611/150‑1 (MQ=255)

                                                                                                                                 gTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGAt                    <  1:611472/150‑1 (MQ=255)

                                                                                                                                 gTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGAt                    <  1:2201073/150‑1 (MQ=255)

                                                                                                                                  ttGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTGGAATGTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGAtt                   <  2:176908/150‑1 (MQ=255)

                                                                                                                                  ttGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGAtt                   <  1:883431/150‑1 (MQ=255)

                                                                                                                                      tAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATc               >  1:1016774/1‑150 (MQ=255)

                                                                                                                                       cgacctcACCCTTACCAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATcc              >  1:1480593/7‑150 (MQ=255)

                                                                                                                                        gagCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCg             >  1:1508044/1‑150 (MQ=255)

                                                                                                                                          gcgcACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCgtg           <  1:324878/150‑1 (MQ=255)

                                                                                                                                          cctcACCCTTACCAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCgtg           >  2:78695/4‑150 (MQ=255)

                                                                                                                                          cctcACCCTTACCAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCgtg           >  1:1591604/4‑150 (MQ=255)

                                                                                                                                            gcACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTc         <  1:185854/150‑1 (MQ=255)

                                                                                                                                                   tGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATc                          <  2:472433/126‑1 (MQ=255)

                                                                                                                                                   tGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATc                          >  1:472433/1‑126 (MQ=255)

GCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAG‑TTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGG  >  NC_000913/4175639‑4175935

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`