BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,950,306	C→G	5.6%	G73R (GGC→CGC)	aspS ←	aspartate‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,950,306	0	C	G	5.6%	485.5 / 14.7	142	G73R (GGC→CGC)	aspS	aspartate‑‑tRNA ligase
Reads supporting (aligned to +/- strand): ref base C (64/70); new base G (5/3); total (69/73)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.36e-01

Only 100 of 146 total aligned reads displayed.

CAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCACGACGACGGTTGACC  >  NC_000913/1950157‑1950454

  gACGCTATTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTATCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCTCA‑CGTACGGTGCCCg                                                                                                                                                    <  2:1371547/150‑1 (MQ=255)

  gACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCg                                                                                                                                                    >  2:1427884/1‑150 (MQ=255)

  gACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCg                                                                                                                                                    >  2:732199/1‑150 (MQ=255)

  gACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCg                                                                                                                                                    <  2:1065799/150‑1 (MQ=255)

  gACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCg                                                                                                                                                    >  1:1221695/1‑150 (MQ=255)

    cgcgCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTg                                                                                                                                                  <  2:948533/150‑1 (MQ=255)

    cgcgCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTg                                                                                                                                                  >  2:1112190/1‑150 (MQ=255)

    cgcgCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTg                                                                                                                                                  >  2:1187979/1‑150 (MQ=255)

     gcgcTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGa                                                                                                                                                 <  2:1311608/150‑1 (MQ=255)

      cgcTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGAc                                                                                                                                                >  1:265651/1‑150 (MQ=255)

       gcTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAcg                                                                                                                                               <  1:875807/150‑2 (MQ=255)

       gcTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGAcc                                                                                                                                               >  1:622698/1‑150 (MQ=255)

        cttcttCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAcga                                                                                                                                              <  1:548282/150‑3 (MQ=255)

         ttcttcGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTg                                                                                                                                             <  1:253520/150‑1 (MQ=255)

           cttcGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGaa                                                                                                                                           <  1:1493534/150‑1 (MQ=255)

             tcGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATg                                                                                                                                         <  2:1519529/150‑1 (MQ=255)

              cGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGc                                                                                                                                        >  2:1473815/1‑150 (MQ=255)

               ggTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATAGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCa                                                                                                                                       >  2:750566/1‑150 (MQ=255)

                gtgtTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAg                                                                                                                                      <  2:1338025/150‑1 (MQ=255)

                   ttGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                                                                                   <  1:1146502/150‑1 (MQ=255)

                     gACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGAA‑CGTACGGTGCCCGTGACCTGAATGCAGAACta                                                                                                                                 >  1:370036/1‑149 (MQ=255)

                      aCGTGGTTAGAGTCAAGCGGCAGAACATGTTGGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCa                                                                                                                                <  1:719691/150‑1 (MQ=255)

                      aCGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCa                                                                                                                                <  1:159229/150‑1 (MQ=255)

                        gTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCAtt                                                                                                                              <  1:1554099/150‑1 (MQ=255)

                          ggTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTAc                                                                                                                            >  1:941720/1‑150 (MQ=255)

                               gagTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGt                                                                                                                       <  1:881728/150‑1 (MQ=255)

                               gagTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGt                                                                                                                       <  1:260820/150‑1 (MQ=255)

                                 gTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTc                                                                                                                     >  2:339310/1‑150 (MQ=255)

                                   cAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGa                                                                                                                                     >  1:120160/1‑132 (MQ=255)

                                   cAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCag                                                                                                                   >  2:421885/1‑150 (MQ=255)

                                   cAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCag                                                                                                                   >  1:625615/1‑150 (MQ=255)

                                    aaGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCaga                                                                                                                  >  1:361741/1‑150 (MQ=255)

                                    aaGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCaga                                                                                                                  <  1:1514060/150‑1 (MQ=255)

                                     aGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCagag                                                                                                                 >  2:543611/1‑150 (MQ=255)

                                     aGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGGCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCagag                                                                                                                 <  2:376718/150‑1 (MQ=255)

                                        ggCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGcc                                                                                                              >  2:1073583/1‑150 (MQ=255)

                                         gCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCa                                                                                                             >  1:1337126/1‑150 (MQ=255)

                                         gCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCa                                                                                                             >  2:1516372/1‑150 (MQ=255)

                                          cataaCATCTTCGAGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCTACGATCGACATATAGAGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAg                                                                                                            >  2:483887/4‑150 (MQ=255)

                                             aaCATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCtt                                                                                                         <  2:75055/150‑1 (MQ=255)

                                               cATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTa                                                                                                       >  1:45922/1‑150 (MQ=255)

                                                 tCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAc                                                                                                     <  1:974368/150‑1 (MQ=255)

                                                          ttGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCa                                                                                            <  1:1473815/150‑1 (MQ=255)

                                                           tgatgaTAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCAc                                                                                           >  1:1200783/1‑150 (MQ=255)

                                                           tgatgaTAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCAc                                                                                           <  1:732199/150‑1 (MQ=255)

                                                            gatgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCTACACGCGCA‑CGTACGGTGCCCGTGACATGAATGCAGAACTCATTACGCAGTAAAGAGGCCAGCTTAAACGCGTCCGCACg                                                                                          >  1:1069570/1‑150 (MQ=255)

                                                              tgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGAt                                                                                        >  1:406290/1‑150 (MQ=255)

                                                                atAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATcc                                                                                      >  1:1651764/1‑150 (MQ=255)

                                                                   gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGa                                                                                   >  1:17564/1‑150 (MQ=255)

                                                                     cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTa                                                                                                                             >  2:1587969/1‑106 (MQ=255)

                                                                     cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTa                                                                                                                             <  1:1587969/106‑1 (MQ=255)

                                                                     cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCagag                                                                                                                 >  1:488120/1‑118 (MQ=255)

                                                                     cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCagag                                                                                                                 <  2:488120/118‑1 (MQ=255)

                                                                       gCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGa                                                                               <  1:1593435/150‑1 (MQ=255)

                                                                         gAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaa                                                                             <  1:543611/150‑1 (MQ=255)

                                                                         gAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaa                                                                             <  1:339310/150‑1 (MQ=255)

                                                                           ggACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaaaa                                                                           >  1:1037835/1‑150 (MQ=255)

                                                                              cGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc                                                                        >  2:63506/1‑150 (MQ=255)

                                                                              cGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATAc                                                                        >  1:114535/1‑150 (MQ=255)

                                                                                ccAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt                                                                      >  1:1085915/1‑150 (MQ=255)

                                                                                       ttCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAACTACCTGCACGAt                                                               >  2:711537/1‑150 (MQ=255)

                                                                                       ttCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGAt                                                               <  1:1109388/150‑1 (MQ=255)

                                                                                        tCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATa                                                              >  2:1574811/1‑150 (MQ=255)

                                                                                          gATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATAcc                                                            >  2:803380/1‑150 (MQ=255)

                                                                                              tCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCCGCTTTAACGCGACCGCACGATCCGGATCGAAAAATACCAGCACGATACCTTcg                                                        >  1:1017564/1‑150 (MQ=255)

                                                                                                gCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCAGCACGATCCGGATCGAAAAATACCTGAACGATACCTTcgcg                                                      >  2:512588/1‑150 (MQ=255)

                                                                                                 ccGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCgg                                                     <  1:77693/150‑1 (MQ=255)

                                                                                                 ccGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCgg                                                     >  1:205393/1‑150 (MQ=255)

                                                                                                      cGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTcgcg                                                <  1:1436958/150‑1 (MQ=255)

                                                                                                       gCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCgcgc                                               <  1:413893/150‑1 (MQ=255)

                                                                                                       gCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCgcgc                                               >  1:678275/1‑150 (MQ=255)

                                                                                                            atCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCAc                                                                  >  2:104991/1‑126 (MQ=255)

                                                                                                            atCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATc                                          <  1:1574811/150‑1 (MQ=255)

                                                                                                             tCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCg                                         >  2:1584145/1‑150 (MQ=255)

                                                                                                                 ggTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGa                                     <  1:128154/150‑1 (MQ=255)

                                                                                                                  gTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGaa                                    >  1:34766/1‑150 (MQ=255)

                                                                                                                      atatTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATc                                >  1:852913/1‑150 (MQ=255)

                                                                                                                      atatTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATc                                <  2:1207309/150‑1 (MQ=255)

                                                                                                                      atatTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATc                                <  2:34766/150‑1 (MQ=255)

                                                                                                                         tttttttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAgg                             <  1:750566/150‑1 (MQ=255)

                                                                                                                          ttttttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGc                            <  1:1048355/150‑1 (MQ=255)

                                                                                                                           tttttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCt                           <  1:1062179/150‑1 (MQ=255)

                                                                                                                           attcaggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCt                           >  1:1332415/7‑150 (MQ=255)

                                                                                                                              ttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTAcc                        <  2:361741/150‑1 (MQ=255)

                                                                                                                               tCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCa                       >  2:83847/1‑150 (MQ=255)

                                                                                                                                 gTCACGGGCACCGTAC‑GTGCGCGTGTCCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGAACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAg                     >  1:751849/1‑150 (MQ=255)

                                                                                                                                 gTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAg                     >  1:1285436/1‑150 (MQ=255)

                                                                                                                                 gTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGc                            >  1:1219294/1‑143 (MQ=255)

                                                                                                                                 gTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGc                            <  2:1219294/143‑1 (MQ=255)

                                                                                                                                 gTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAg                     <  2:211490/150‑1 (MQ=255)

                                                                                                                                 gTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAg                     >  1:1232460/1‑150 (MQ=255)

                                                                                                                                 gTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAg                     >  1:774662/1‑150 (MQ=255)

                                                                                                                                  tCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGa                    >  1:677115/1‑150 (MQ=255)

                                                                                                                                       cgcA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCacg               <  1:483887/150‑1 (MQ=255)

                                                                                                                                             gTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCacgacgac          <  2:535735/150‑1 (MQ=255)

                                                                                                                                                cGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCACGACGACGGt       >  2:789316/1‑150 (MQ=255)

                                                                                                                                                    gCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCACGACGACGGTTGAc   <  2:622698/150‑1 (MQ=255)

                                                                                                                                                    gCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCACGACGACGGTTGAc   <  2:1539533/150‑1 (MQ=255)

                                                                                                                                                     cccGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCACGACGACGGTTGAcc  <  2:45922/150‑1 (MQ=255)

                                                                                                                                                     cccGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCACGACGACGGTTGAcc  <  2:17564/150‑1 (MQ=255)

CAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCACGACGACGGTTGACC  >  NC_000913/1950157‑1950454

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: `‑`