BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,935,895	A→G	T15A (ACC→GCC)	rbsR →	RbsR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,935,895	0	A	G	85.0%	701.7 / 103.6	300	T15A (ACC→GCC)	rbsR	RbsR
Reads supporting (aligned to +/- strand): ref base A (19/26); new base G (138/117); total (157/143)
Fisher's exact test for biased strand distribution p-value = 1.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 301 total aligned reads displayed.

TACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTA  >  NC_000913/3935803‑3935986

  ccGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGcc                                                                                           >  1:4207189/1‑93 (MQ=255)

  ccGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGcc                                                                                           >  1:4555692/1‑93 (MQ=255)

  ccGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGcc                                                                                           >  2:2246748/1‑93 (MQ=255)

  ccGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGcc                                                                                           <  2:17674/93‑1 (MQ=255)

      ggCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCaa                                                                                       <  1:909040/93‑1 (MQ=255)

       gCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAAc                                                                                      <  2:2290418/93‑1 (MQ=255)

        cGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACa                                                                                     <  1:2563589/93‑1 (MQ=255)

        cGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACa                                                                                     >  1:3204459/1‑93 (MQ=255)

          tGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGt                                                                                   >  2:7045066/1‑93 (MQ=255)

           gAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGtt                                                                                  >  2:7398976/1‑93 (MQ=255)

           gAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGtt                                                                                  >  1:4864763/1‑93 (MQ=255)

           gAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGtt                                                                                  <  1:6414732/93‑1 (MQ=255)

           gAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGtt                                                                                  <  2:612190/93‑1 (MQ=255)

            aaGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCCGGCGGGCGTTTCTACCTCAACAGttt                                                                                 >  1:1278973/1‑93 (MQ=255)

             agagATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTtc                                                                                <  1:551693/93‑1 (MQ=255)

             agagATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTtc                                                                                <  1:7011731/93‑1 (MQ=255)

                gaTCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCa                                                                             >  2:6628285/1‑93 (MQ=255)

                    gACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAGCAGTTTCTCACGtt                                                                         >  2:4381492/1‑93 (MQ=255)

                     aCGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTa                                                                        >  1:3641101/1‑93 (MQ=255)

                       gCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATc                                                                      >  2:1204773/1‑93 (MQ=255)

                       gCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATc                                                                      <  1:6026639/93‑1 (MQ=255)

                         aTTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTATCTCACGTTATCaa                                                                    <  2:6453124/93‑1 (MQ=255)

                          tttttAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAAt                                                                   >  1:1247533/1‑93 (MQ=255)

                          tttttAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAAt                                                                   >  1:5133934/1‑93 (MQ=255)

                          tttttAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCCGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAAt                                                                   <  1:7189926/93‑1 (MQ=255)

                           ttttAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCCGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATa                                                                  <  2:3210639/93‑1 (MQ=255)

                            tttAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATaa                                                                 >  2:5462954/1‑93 (MQ=255)

                            tttAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATaa                                                                 <  2:2465571/93‑1 (MQ=255)

                             ttAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATaaa                                                                >  1:2878142/1‑93 (MQ=255)

                             ttAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCCGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATaaa                                                                <  2:349723/93‑1 (MQ=255)

                                gACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGAt                                                             <  1:6261867/93‑1 (MQ=255)

                                gACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCCGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGAt                                                             >  1:4081498/1‑93 (MQ=255)

                                  cAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCg                                                           <  1:2236533/93‑1 (MQ=255)

                                  cAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCg                                                           <  2:5244436/93‑1 (MQ=255)

                                   aGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGc                                                          >  1:478201/1‑93 (MQ=255)

                                   aGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGc                                                          >  2:3643855/1‑93 (MQ=255)

                                     gCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCCGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCtt                                                        <  2:5563353/93‑1 (MQ=255)

                                       agagGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCg                                                      >  2:3070810/1‑93 (MQ=255)

                                       agagGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCg                                                      >  1:4612449/1‑93 (MQ=255)

                                       agagGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCg                                                      >  2:4086034/1‑93 (MQ=255)

                                          ggTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCa                                                   <  1:4015177/93‑1 (MQ=255)

                                          ggTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCa                                                   >  2:6260005/1‑93 (MQ=255)

                                          ggTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCa                                                   <  1:2951585/93‑1 (MQ=255)

                                          ggTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCa                                                   <  2:1577520/93‑1 (MQ=255)

                                              aCGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGa                                               >  1:5852158/1‑93 (MQ=255)

                                              aCGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGa                                               <  2:6573189/93‑1 (MQ=255)

                                              aCGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGa                                               >  1:2736896/1‑93 (MQ=255)

                                              aCGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGa                                               <  2:1463218/93‑1 (MQ=255)

                                               cGCTTGGCTACAATGAAAGATGTTGCCCGCCCGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGaa                                              <  2:5010717/93‑1 (MQ=255)

                                               cGCTTGGCTACAATGAAAGATGTTGCCCGCCCGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGaa                                              <  1:1033333/93‑1 (MQ=255)

                                                gCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGAGAAg                                             <  2:768577/93‑1 (MQ=255)

                                                 cTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGc                                            <  1:5348866/93‑1 (MQ=255)

                                                  ttGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCg                                           >  2:6222829/1‑93 (MQ=255)

                                                  ttGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCg                                           >  1:3409456/1‑93 (MQ=255)

                                                  ttGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGGTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCg                                           >  1:4458416/1‑93 (MQ=255)

                                                    ggCTACAATGAAAGATGTTGCCCGCCCGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGAt                                         <  1:6344564/93‑1 (MQ=255)

                                                     gCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGAtt                                        >  2:1609057/1‑93 (MQ=255)

                                                     gCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGAtt                                        <  2:3389285/93‑1 (MQ=255)

                                                       tACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTAc                                      <  2:4544384/93‑1 (MQ=255)

                                                        aCAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTAcc                                     >  1:438471/1‑93 (MQ=255)

                                                        aCAATGAAAGATGTTGCCCGCCCGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTAcc                                     >  2:5967163/1‑93 (MQ=255)

                                                          aaTGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGc                                   >  1:1729593/1‑93 (MQ=255)

                                                            tGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCa                                 >  2:6885814/1‑93 (MQ=255)

                                                             gAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCaa                                <  1:3052461/93‑1 (MQ=255)

                                                             gAAAGATGTTGCCCGCCCGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCaa                                >  1:3073579/1‑93 (MQ=255)

                                                              aaaGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGTTTCGTCAGTGAAGCGATTACCGCCaaa                               <  2:5171568/93‑1 (MQ=255)

                                                              aaaGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCaaa                               <  2:5079120/93‑1 (MQ=255)

                                                              aaaGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCaaa                               <  2:3186695/93‑1 (MQ=255)

                                                              aaaGATGTAGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCaca                               >  1:7408050/1‑91 (MQ=255)

                                                               aaGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATTAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAg                              >  2:4809396/1‑93 (MQ=255)

                                                                aGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGt                             <  2:729742/93‑1 (MQ=255)

                                                                 gATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGtt                            >  1:569991/1‑93 (MQ=255)

                                                                  aTGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTg                           <  2:3547033/93‑1 (MQ=255)

                                                                  aTGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTg                           <  2:2088438/93‑1 (MQ=255)

                                                                   tGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGa                          >  2:5673991/1‑93 (MQ=255)

                                                                   tGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGa                          >  2:2943222/1‑93 (MQ=255)

                                                                   tGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGa                          >  2:3774018/1‑93 (MQ=255)

                                                                    gTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGaa                         <  2:1528658/93‑1 (MQ=255)

                                                                     ttGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAg                        >  2:3906419/1‑93 (MQ=255)

                                                                      tGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAgc                       <  2:6980609/93‑1 (MQ=255)

                                                                      tGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAgc                       >  2:1815890/1‑93 (MQ=255)

                                                                       gCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAgcg                      <  1:6772977/93‑1 (MQ=255)

                                                                        cccGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAgcgg                     >  1:1163734/1‑93 (MQ=255)

                                                                         ccGCCCGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAgcggc                    <  2:4919810/93‑1 (MQ=255)

                                                                           gCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGa                  >  2:1040517/1‑93 (MQ=255)

                                                                           gCCTGGCGGGCGTTTCTGCCCCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGa                  <  1:7506089/93‑1 (MQ=255)

                                                                             cTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGAtt                >  1:5576147/1‑93 (MQ=255)

                                                                             cTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGAtt                >  1:5954598/1‑93 (MQ=255)

                                                                             cTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGAtt                >  1:6469369/1‑93 (MQ=255)

                                                                              tGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTa               >  1:5076988/1‑93 (MQ=255)

                                                                              tGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTa               >  2:6723955/1‑93 (MQ=255)

                                                                               ggcgggcgTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTaa              <  2:4222906/93‑1 (MQ=255)

                                                                                  gggcgTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGa           <  2:1984546/93‑1 (MQ=255)

                                                                                   ggcgTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGaa          >  1:3776489/1‑93 (MQ=255)

                                                                                   ggcgTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGaa          >  2:2712206/1‑93 (MQ=255)

                                                                                   ggcgTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGaa          >  1:3929993/1‑93 (MQ=255)

                                                                                    gcgTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAAc         <  1:7229022/93‑1 (MQ=255)

                                                                                      gTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTc       <  1:6947876/93‑1 (MQ=255)

                                                                                      gTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTc       <  1:1988296/93‑1 (MQ=255)

                                                                                         tCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAAt    >  1:2503544/1‑93 (MQ=255)

TACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTA  >  NC_000913/3935803‑3935986

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