BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	318,936	C→T	G106D (GGC→GAC)	ykgC ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	318,936	0	C	T	99.1%	815.5 / ‑4.6	216	G106D (GGC→GAC)	ykgC	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base C (1/1); new base T (100/114); total (101/115)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 223 total aligned reads displayed.

TATTAATAAAAATTTTCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGCCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCATTTTTACGCTGTATG  >  NC_000913/318847‑319025

tATTAATAAAAATTTTCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGt                                                                                        <  1:6382492/93‑1 (MQ=255)

tATTAATAAAAATTTTCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGt                                                                                        <  1:699867/93‑1 (MQ=255)

 aTTAATAAAAATTTTCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGgtcgtc                                                                                       <  1:3196351/93‑1 (MQ=255)

  ttAATAAAAATTTTCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGtcgtcg                                                                                      <  1:733451/93‑1 (MQ=255)

    aataaAAATTTTCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGAt                                                                                    <  2:5932497/93‑1 (MQ=255)

    aataaAAATTTTCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGAt                                                                                    <  2:1807146/93‑1 (MQ=255)

     ataaAAATTTTCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATc                                                                                   <  1:1893859/93‑1 (MQ=255)

      taaAAATTTTCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCa                                                                                  >  2:1192420/1‑93 (MQ=255)

       aaaaaTTTTCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCAc                                                                                 <  2:6204537/93‑1 (MQ=255)

       aaaaaTTTTCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCAc                                                                                 <  1:896128/93‑1 (MQ=255)

       aaaaaTTTTCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCAc                                                                                 >  2:2447573/1‑93 (MQ=255)

         aaaTTTTCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGt                                                                               >  1:1244050/1‑93 (MQ=255)

            ttttCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGa                                                                            <  2:3041814/93‑1 (MQ=255)

            ttttCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGa                                                                            <  1:4059917/93‑1 (MQ=255)

                ctcGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATAtt                                                                        >  1:2036871/1‑93 (MQ=255)

                ctcGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATAtt                                                                        <  1:5784756/93‑1 (MQ=255)

                 tcGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTg                                                                       >  2:1916003/1‑93 (MQ=255)

                 tcGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTg                                                                       <  1:4902615/93‑1 (MQ=255)

                 tcGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTg                                                                       >  1:3839378/1‑93 (MQ=255)

                  cGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATATACTCCGCCTGGGCGTCGATTACGTCGATATTgg                                                                      >  1:535197/1‑93 (MQ=255)

                     cATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCa                                                                   <  1:2596957/93‑1 (MQ=255)

                       tGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCata                                                                 >  2:792541/1‑93 (MQ=255)

                       tGAATTTCCAGATTCCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCata                                                                 <  2:5136202/93‑1 (MQ=255)

                        gAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCatat                                                                <  2:2931493/93‑1 (MQ=255)

                         aaTTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATc                                                               <  1:3678033/93‑1 (MQ=255)

                         aaTTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATc                                                               <  1:5914077/93‑1 (MQ=255)

                          aTTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATcc                                                              <  2:4409470/93‑1 (MQ=255)

                          aTTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATcc                                                              >  1:3955805/1‑93 (MQ=255)

                           tttCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCg                                                             >  1:2979638/1‑93 (MQ=255)

                             tCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCa                                                           >  2:702267/1‑93 (MQ=255)

                             tCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCa                                                           >  1:6022494/1‑93 (MQ=255)

                              ccAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCaa                                                          <  1:4789683/93‑1 (MQ=255)

                              ccAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCaa                                                          <  2:505086/93‑1 (MQ=255)

                              ccAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCaa                                                          >  2:1096804/1‑93 (MQ=255)

                               cAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAg                                                         >  1:446178/1‑93 (MQ=255)

                                aGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGa                                                        <  2:1691045/93‑1 (MQ=255)

                                aGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGa                                                        <  1:6197797/93‑1 (MQ=255)

                                 gATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGAt                                                       >  1:6433110/1‑93 (MQ=255)

                                 gATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGAt                                                       <  1:1154425/93‑1 (MQ=255)

                                 gATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGAt                                                       <  1:2969657/93‑1 (MQ=255)

                                 gATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGAt                                                       <  2:4855841/93‑1 (MQ=255)

                                 gATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGAt                                                       >  2:540280/1‑93 (MQ=255)

                                  aTTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGAtt                                                      <  2:3446018/93‑1 (MQ=255)

                                  aTTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGAtt                                                      <  1:374681/93‑1 (MQ=255)

                                  aTTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGAtt                                                      <  2:3839378/93‑1 (MQ=255)

                                  aTTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGAtt                                                      >  1:5439239/1‑93 (MQ=255)

                                  aTTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGAtt                                                      >  1:4174406/1‑93 (MQ=255)

                                   tttCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTa                                                     <  1:2525598/93‑1 (MQ=255)

                                    ttCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTAt                                                    <  1:1337545/93‑1 (MQ=255)

                                     tCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATg                                                   >  2:1021030/1‑93 (MQ=255)

                                      cccTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGa                                                  >  2:4696822/1‑93 (MQ=255)

                                       ccTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGaa                                                 >  1:4707083/1‑93 (MQ=255)

                                       ccTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGaa                                                 >  2:3553926/1‑93 (MQ=255)

                                          cAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAAt                                              <  2:5332316/93‑1 (MQ=255)

                                           aGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGaaaagt                                             >  2:4336478/1‑91 (MQ=255)

                                           aGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAAtt                                             >  2:2399052/1‑93 (MQ=255)

                                             gCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTTGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCt                                           >  2:4169555/1‑93 (MQ=255)

                                             gCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGATAATTCt                                           <  1:5359476/93‑1 (MQ=255)

                                             gCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCt                                           >  2:2126596/1‑93 (MQ=255)

                                                gATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCttat                                        >  1:1023831/1‑93 (MQ=255)

                                                gATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCttat                                        >  1:6132152/1‑93 (MQ=255)

                                                caTGACCACGCAGGCTATGATTATTGATAACCTCCGCCTGGTAGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCttat                                        <  1:5801551/92‑1 (MQ=255)

                                                  tGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCttatta                                      >  2:4962367/1‑93 (MQ=255)

                                                   gAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTAc                                     <  2:4108913/93‑1 (MQ=255)

                                                     acacGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGt                                   <  1:2364650/93‑1 (MQ=255)

                                                       acGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTaa                                 <  1:2999643/93‑1 (MQ=255)

                                                       acGCAGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTaa                                 >  2:3697866/1‑93 (MQ=255)

                                                           aGGCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTaaaaaa                             <  1:3341903/93‑1 (MQ=255)

                                                            ggCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAAt                            >  1:5231293/1‑93 (MQ=255)

                                                             gCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAAtt                           <  1:2173775/93‑1 (MQ=255)

                                                             gCTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAAtt                           >  2:4495395/1‑93 (MQ=255)

                                                              cTATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTa                          >  2:2464017/1‑93 (MQ=255)

                                                               tATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAACATTaa                         >  2:4399025/1‑93 (MQ=255)

                                                               tATGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTaa                         <  2:13689/93‑1 (MQ=255)

                                                               tATGATTATTGATAAACTCCGCCTGGCCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTaa                         <  2:2609137/93‑1 (MQ=255)

                                                                aTGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAACAAATTAAc                        <  1:6321069/93‑1 (MQ=255)

                                                                aTGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAAc                        >  2:6217141/1‑93 (MQ=255)

                                                                 tGATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAAcc                       <  1:3461871/93‑1 (MQ=255)

                                                                  gATTATTGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCa                      >  2:4563746/1‑93 (MQ=255)

                                                                   attattGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAATTTCTTATTACGTAAAAAATTAACCAc                     <  1:4509928/93‑1 (MQ=255)

                                                                   attattGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCAc                     <  1:1812833/93‑1 (MQ=255)

                                                                   attattGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCAc                     >  2:225908/1‑93 (MQ=255)

                                                                    ttattGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACt                    <  1:5586883/93‑1 (MQ=255)

                                                                    ttattGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACt                    <  2:5621401/93‑1 (MQ=255)

                                                                     tattGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACtt                   <  2:3238923/93‑1 (MQ=255)

                                                                       ttGATAACCTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCa                 <  1:2298119/93‑1 (MQ=255)

                                                                       ttGATAACCTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCa                 <  2:2343748/93‑1 (MQ=255)

                                                                       ttGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCa                 <  1:5060340/93‑1 (MQ=255)

                                                                        tGATAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCAt                <  2:3374311/93‑1 (MQ=255)

                                                                          aTAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCAttt              <  1:1825149/93‑1 (MQ=255)

                                                                           tAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCAtttt             >  2:1565299/1‑93 (MQ=255)

                                                                           tAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCAtttt             >  2:3926272/1‑93 (MQ=255)

                                                                           tAAACTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCAtttt             <  2:5791358/93‑1 (MQ=255)

                                                                             aaCTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCATTTTTa           >  2:5268938/1‑93 (MQ=255)

                                                                              aCTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCATTTTTAc          >  1:3075289/1‑93 (MQ=255)

                                                                              aCTCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCATTTTTAc          >  2:473134/1‑93 (MQ=255)

                                                                                tCCGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCATTTTTACGc        <  1:2026333/93‑1 (MQ=255)

                                                                                  cGCCTGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCATTTTTACGCTg      >  1:4633896/1‑93 (MQ=255)

                                                                                      tGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCATTTTTACGCTGTATg  <  1:5698434/93‑1 (MQ=255)

                                                                                      tGGTCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCATTTTTACGCTGTATg  >  2:634922/1‑93 (MQ=255)

TATTAATAAAAATTTTCTCGCCATGAATTTCCAGATTTCCCTCAGGCCGATGAACACGCAGGCTATGATTATTGATAAACTCCGCCTGGCCGTCGATCACGTCGATATTGGGCATATCCGCAAGATTATGAAAATTCTTATTACGTAAAAAATTAACCACTTCATTTTTACGCTGTATG  >  NC_000913/318847‑319025

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