BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,735,666	T→C	S181P (TCC→CCC)	malS →	alpha‑amylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,735,666	0	T	C	99.1%	720.6 / NA	214	S181P (TCC→CCC)	malS	alpha‑amylase
Reads supporting (aligned to +/- strand): ref base T (2/0); new base C (112/100); total (114/100)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 231 total aligned reads displayed.

ATTGCGATAGTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACATCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGACCAGAGTTAC  >  NC_000913/3735576‑3735755

aTTGCGATAGTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACAccc                                                                                         <  1:3856131/93‑1 (MQ=255)

aTTGCGATAGTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACAccc                                                                                         <  1:5603937/93‑1 (MQ=255)

 ttGCGATAGTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCg                                                                                        >  1:1995840/1‑93 (MQ=255)

  tGCGATAGTGAAGAACGGTCAAATAATGTTACACCCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGc                                                                                       <  2:1566674/93‑1 (MQ=255)

  tGCAATGGTAAAGAACGGTCAAATAATGTTACACCCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGc                                                                                       <  1:4833029/93‑1 (MQ=255)

    cGATAGTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGccc                                                                                     <  2:2290700/93‑1 (MQ=255)

    cGATAATGAAGAGCGGTCAAATAATGTTACACCCCGCTGCCACCAGCACCGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGccc                                                                                     <  2:3225545/93‑1 (MQ=255)

     gATAGTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGcccc                                                                                    <  2:5296580/93‑1 (MQ=255)

      aTAGTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCt                                                                                   <  1:3990640/93‑1 (MQ=255)

       aaGTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCtt                                                                                  <  2:4114953/92‑1 (MQ=255)

        aGTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCttt                                                                                 <  2:413623/93‑1 (MQ=255)

         gtaaagaaCGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTc                                                                                <  2:4775067/90‑1 (MQ=255)

         gTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTGCTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTc                                                                                <  1:3091242/93‑1 (MQ=255)

         gTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACATCCGCTCCTTTc                                                                                >  1:2498361/1‑93 (MQ=255)

         gTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTc                                                                                >  1:3132924/1‑93 (MQ=255)

         gTGAAGAACGGTCAAATAATGTTACAAACCGCTGCCACCAGCAACGGATTACTCCTGCTGGAACGGGTAGTAAATGAGACACCCACACCTTTc                                                                                >  1:1511470/1‑93 (MQ=255)

           gaagaaCGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGa                                                                              <  2:3791424/93‑1 (MQ=255)

           gaagaaCGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGa                                                                              >  1:3569057/1‑93 (MQ=255)

            aagaaCGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGAc                                                                             >  1:6020226/1‑93 (MQ=255)

            aagaaCGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGAc                                                                             <  1:1576647/93‑1 (MQ=255)

             agaaCGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACt                                                                            >  2:6025941/1‑93 (MQ=255)

               aaCGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTgg                                                                          <  2:2857153/93‑1 (MQ=255)

                aCGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGc                                                                         <  1:3340558/93‑1 (MQ=255)

                  ggTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCAt                                                                       >  1:1692901/1‑93 (MQ=255)

                  ggTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCAt                                                                       <  2:5621802/93‑1 (MQ=255)

                    tCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATaa                                                                     >  2:3519669/1‑93 (MQ=255)

                     cAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAAc                                                                    >  1:2891787/1‑93 (MQ=255)

                     cAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAAc                                                                    >  2:5577423/1‑93 (MQ=255)

                     cAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAAc                                                                    >  2:5390670/1‑93 (MQ=255)

                      aaaTAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACg                                                                   >  2:1537715/1‑93 (MQ=255)

                      aaaTAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACg                                                                   >  2:200445/1‑93 (MQ=255)

                       aataatGTTACAACCCGCTGCCACCAGCAACGGTTCACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGc                                                                  <  1:5186479/93‑1 (MQ=255)

                       aataatGTAACAACCCGCTGCCGCCAACAACGGTTTGCTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGc                                                                  <  1:4369261/93‑1 (MQ=255)

                        ataatGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACATCCGCCCCTTTCGACTGGCATAACGcc                                                                 >  2:6067012/1‑93 (MQ=255)

                        ataatGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGcc                                                                 >  2:188467/1‑93 (MQ=255)

                         taatGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCa                                                                >  2:3851420/1‑93 (MQ=255)

                         taatGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCa                                                                >  1:183230/1‑93 (MQ=255)

                         taatGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCa                                                                <  2:5701832/93‑1 (MQ=255)

                          aatGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCAc                                                               <  2:2978038/93‑1 (MQ=255)

                          aatGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCAc                                                               <  1:2545111/93‑1 (MQ=255)

                          aatGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCAc                                                               >  1:3492486/1‑93 (MQ=255)

                          aatATTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCAc                                                               >  1:5333829/1‑93 (MQ=255)

                           atGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACg                                                              >  1:2401277/1‑93 (MQ=255)

                            tGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATACCGCCACgg                                                             <  1:6311975/93‑1 (MQ=255)

                              ttACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGCCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGtt                                                           <  2:2063063/93‑1 (MQ=255)

                               tccAACCGGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAAAAACTGACACCCCCGCCCCTTTCGACTGGCATAACGCCACGGttt                                                          <  2:148641/91‑1 (MQ=255)

                               tACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGttt                                                          >  2:825137/1‑93 (MQ=255)

                                 caccccGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTAc                                                        <  2:5710587/90‑1 (MQ=255)

                                 caccccGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTAc                                                        <  2:6368367/90‑1 (MQ=255)

                                  aaCCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACt                                                       >  2:6472923/1‑93 (MQ=255)

                                   aCCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACCCCACGGTTTACtt                                                      >  1:3412346/1‑93 (MQ=255)

                                    cccGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACttt                                                     <  1:1846355/93‑1 (MQ=255)

                                     ccGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTtg                                                    >  1:1717327/1‑93 (MQ=255)

                                      cGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTtgt                                                   <  1:4886003/93‑1 (MQ=255)

                                      cGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTtgt                                                   <  2:1387589/93‑1 (MQ=255)

                                      cGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTtgt                                                   <  1:278315/93‑1 (MQ=255)

                                         tGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCt                                                <  2:1772430/93‑1 (MQ=255)

                                          gCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTg                                               >  1:3667703/1‑93 (MQ=255)

                                             accaGCAACGGTTTACTCCTGCTGGACCGGGGAGAAACTGACACACCCGCTCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACa                                            <  2:4520004/93‑1 (MQ=255)

                                             accaGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACa                                            >  2:2521064/1‑93 (MQ=255)

                                              ccaGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAg                                           >  1:3141059/1‑93 (MQ=255)

                                              ccaGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAg                                           >  2:4672986/1‑93 (MQ=255)

                                                aGCAACTGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGAt                                         <  2:3093756/93‑1 (MQ=255)

                                                aGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGAt                                         <  2:6231333/93‑1 (MQ=255)

                                                aGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGAt                                         >  2:1191163/1‑93 (MQ=255)

                                                aGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGAt                                         <  1:1192549/93‑1 (MQ=255)

                                                 gCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATc                                        >  2:2952172/1‑93 (MQ=255)

                                                 gCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATc                                        >  2:3045913/1‑93 (MQ=255)

                                                 gCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATc                                        <  2:5450198/93‑1 (MQ=255)

                                                     cGGTTTCCTCCTGCTGGAACGGGCAGAAACTGACACCCCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGttt                                    <  1:648904/93‑1 (MQ=255)

                                                      ggTTTACTCCTGCTGGAACGGGCAGAAACTGACACCCCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTc                                   <  2:4318021/93‑1 (MQ=255)

                                                       gTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCg                                  >  1:6085692/1‑93 (MQ=255)

                                                        tttACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGa                                 >  2:5508549/1‑93 (MQ=255)

                                                        tttACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGa                                 >  1:4705914/1‑93 (MQ=255)

                                                            cTCCTGCTGGAACGGGCAGAAACTGACGCACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAAc                             <  2:2496991/93‑1 (MQ=255)

                                                            cTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAAc                             <  2:2326531/93‑1 (MQ=255)

                                                             tCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACg                            >  1:1810034/1‑93 (MQ=255)

                                                             tCCTGCGGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACg                            <  2:1375924/93‑1 (MQ=255)

                                                               ctgctgGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGc                          >  2:3717603/1‑93 (MQ=255)

                                                               ctgctgGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGc                          >  2:1709833/1‑93 (MQ=255)

                                                               ctgctgGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGc                          >  1:4204607/1‑93 (MQ=255)

                                                                tgctgGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCg                         <  2:5239025/93‑1 (MQ=255)

                                                                tgctgGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCg                         <  1:274607/93‑1 (MQ=255)

                                                                 gctgGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGa                        >  1:2273537/1‑93 (MQ=255)

                                                                  ctgGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGAt                       <  1:5234211/93‑1 (MQ=255)

                                                                   tgGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATc                      >  2:1031689/1‑93 (MQ=255)

                                                                     gAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATccc                    >  1:6090451/1‑93 (MQ=255)

                                                                       aCGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAg                  <  1:5727010/93‑1 (MQ=255)

                                                                       aCGGACAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAg                  <  1:2678077/93‑1 (MQ=255)

                                                                         gggCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTa                <  1:881440/93‑1 (MQ=255)

                                                                          ggCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTaa               >  1:1157801/1‑93 (MQ=255)

                                                                            cAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATg             >  2:1318533/1‑93 (MQ=255)

                                                                             aGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGa            >  2:338252/1‑93 (MQ=255)

                                                                             aGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGa            <  1:131135/93‑1 (MQ=255)

                                                                                  cTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGACCaga       <  1:2205269/93‑1 (MQ=255)

                                                                                    gACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGACGAGAGt     >  1:3367653/1‑93 (MQ=255)

                                                                                    gACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGACCAGAGt     <  1:4640682/93‑1 (MQ=255)

                                                                                    gACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGACCAGAGt     >  2:2126218/1‑93 (MQ=255)

                                                                                      cacaCCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGACCAGAGTTa   >  1:4510477/1‑93 (MQ=255)

                                                                                      cacaCCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGACCAGAGTTa   <  2:5836702/93‑1 (MQ=255)

ATTGCGATAGTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACATCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGACCAGAGTTAC  >  NC_000913/3735576‑3735755

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