BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,210,124	T→C	E31E (GAA→GAG)	yehS ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,210,124	0	T	C	99.2%	767.4 / ‑3.7	247	E31E (GAA→GAG)	yehS	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (117/127); total (118/129)
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.68e-01

Only 100 of 253 total aligned reads displayed.

GAGGAATGACGACAAAACAATGTCCGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCTTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTCAAAATGTAGCGCACGCTGCGCAGAATATCGTTACTTA  >  NC_000913/2210034‑2210212

gAGGAATGACGACAAAACAATGTCCGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTc                                                                                        <  2:6578918/93‑1 (MQ=255)

gAGGAATGACGACAAAACAATGTCCGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTc                                                                                        >  1:7255236/1‑93 (MQ=255)

  ggAATGACGACAAAACAATGTCCGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGa                                                                                      <  1:422499/93‑1 (MQ=255)

  gcaaTGACGACAAAACAATGTCCGGACAACGCTGAACCCTCTCTCCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGa                                                                                      <  2:7296229/91‑1 (MQ=255)

   gAATGACGACAAAACAATGTCCGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGAc                                                                                     >  2:1548487/1‑93 (MQ=255)

    aaTGACGACAAAACAATGTCCGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGAct                                                                                    >  1:3581854/1‑92 (MQ=255)

    aaTGACGACAAAACAATGTCCGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACa                                                                                    <  2:7883263/93‑1 (MQ=255)

    aaTGACGACAAAACAATGTCCGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACa                                                                                    <  1:7497568/93‑1 (MQ=255)

     aTGACGACAAAACAATGTCCGGATAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACAt                                                                                   <  2:4297727/93‑1 (MQ=255)

     aTGACGACAAAACAATGTCCGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACAt                                                                                   >  1:4469712/1‑93 (MQ=255)

       gacgacAAAACAATGTCCGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTa                                                                                 >  2:3481927/1‑93 (MQ=255)

       gacgacAAAACAATGACAGGACAACGCTGAAAACCCTCTTCGTATTCATTACGTAGCCAGACGGCGATCTGTTCCGAGGTGGCATAGACATTa                                                                                 <  2:7029365/93‑1 (MQ=255)

        acgacAAAACAATGTCCGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTAc                                                                                <  1:3736964/93‑1 (MQ=255)

        acgacAAAACAATGTCCGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTAc                                                                                >  1:3195253/1‑93 (MQ=255)

         cgacAAAACAATGTCCGGACAACGCTGAAATCCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTAcc                                                                               <  1:2249735/93‑1 (MQ=255)

          gacAAAACAATGTCCGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTAccc                                                                              <  2:140102/93‑1 (MQ=255)

          gacAAAACAATGTCCGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTAccc                                                                              >  1:3540469/1‑93 (MQ=255)

          gacAAAACAATGTCCGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTAccc                                                                              <  1:2056431/93‑1 (MQ=255)

           acAAAACAATGTCCGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCa                                                                             <  1:2104004/93‑1 (MQ=255)

            cAAAACAATGTCCGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAg                                                                            <  2:3027395/93‑1 (MQ=255)

            cAAAACAAGGTCCGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCGGTTCCGCGGTGGCCTCGACATTACCCAg                                                                            <  2:4785825/93‑1 (MQ=255)

              aaaCAATGTCCGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAgcg                                                                          >  2:6851065/1‑93 (MQ=255)

              aaaCAATGTCCGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAgcg                                                                          <  2:7213623/93‑1 (MQ=255)

              aaaCAATGTCCGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAgcg                                                                          >  1:5933249/1‑93 (MQ=255)

              aaaCAATGTCCGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAgcc                                                                          >  1:1207108/1‑92 (MQ=255)

              aaaCAATGTCCGGACAACGCTGAAAACCCCCTTCGTCTTCTTTACGTAGCCAGACGGCGTTCTGTTCCGCGTTGGCCTCGACATTACCCAgcg                                                                          >  1:1538258/1‑93 (MQ=255)

                 cAATGTCCGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCa                                                                       <  2:4743696/93‑1 (MQ=255)

                  aaTGTCCGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCACCGGCAg                                                                      >  1:265873/1‑93 (MQ=255)

                   aTGTCCGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGa                                                                     <  2:3658811/93‑1 (MQ=255)

                     gTCCGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAAt                                                                   >  1:6386480/1‑93 (MQ=255)

                      tCCGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATa                                                                  <  1:7845674/93‑1 (MQ=255)

                      tCCGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATa                                                                  <  2:6312480/93‑1 (MQ=255)

                      gccGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATa                                                                  <  2:2299902/92‑1 (MQ=255)

                       ccGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATAc                                                                 <  2:5785946/93‑1 (MQ=255)

                        cGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGGCTCGACATTACGCAGCGCCAGAATACg                                                                >  2:4467170/1‑93 (MQ=255)

                        cGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACg                                                                <  2:7766593/93‑1 (MQ=255)

                        cGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACg                                                                >  2:4467172/1‑93 (MQ=255)

                         ggACAACGCTGAAACCCCTCTTCGTCTTCTGTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCTTCGACATTACCCAGCGCCAGAATACGc                                                               <  2:3181998/93‑1 (MQ=255)

                         ggACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGc                                                               <  1:4259424/93‑1 (MQ=255)

                          gACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCa                                                              >  2:4225157/1‑93 (MQ=255)

                          gACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCa                                                              <  1:175676/93‑1 (MQ=255)

                          gACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCa                                                              >  1:3078999/1‑93 (MQ=255)

                            cAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCAcc                                                            <  1:5088948/93‑1 (MQ=255)

                            cAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCAcc                                                            <  1:4374011/93‑1 (MQ=255)

                             aaCGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCa                                                           <  2:6769688/93‑1 (MQ=255)

                             aaCGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCa                                                           >  2:6006855/1‑93 (MQ=255)

                              aCGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAg                                                          <  2:4135870/93‑1 (MQ=255)

                               cGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAgg                                                         >  1:5304343/1‑93 (MQ=255)

                               cGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAgg                                                         <  1:7782590/93‑1 (MQ=255)

                                gCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGt                                                        >  2:3811686/1‑93 (MQ=255)

                                 cTGAAACCCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTc                                                       <  2:6523513/93‑1 (MQ=255)

                                 cTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTc                                                       <  1:4627821/93‑1 (MQ=255)

                                 cTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTc                                                       <  1:4981331/93‑1 (MQ=255)

                                 cTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTc                                                       <  1:2431064/93‑1 (MQ=255)

                                     aaaCCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCatta                                                   >  1:954171/1‑93 (MQ=255)

                                      aaCCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCattat                                                  >  2:1209869/1‑93 (MQ=255)

                                       aCCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCattatt                                                 >  1:2763722/1‑93 (MQ=255)

                                        cccTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCattatta                                                <  1:6775733/93‑1 (MQ=255)

                                        cccTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCattatta                                                <  2:536002/93‑1 (MQ=255)

                                          ctctTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCattattatt                                              <  1:6502877/93‑1 (MQ=255)

                                          ctctTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCattattatt                                              <  1:4816387/93‑1 (MQ=255)

                                           tctTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTg                                             >  2:5298112/1‑93 (MQ=255)

                                           tctTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTg                                             >  2:1710125/1‑93 (MQ=255)

                                            tttCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTgg                                            <  1:7067512/92‑1 (MQ=255)

                                             ttCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGc                                           >  2:7549145/1‑93 (MQ=255)

                                             ttCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGc                                           >  2:7074873/1‑93 (MQ=255)

                                              tCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCt                                          >  2:5956818/1‑93 (MQ=255)

                                              tCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCt                                          >  1:7577183/1‑93 (MQ=255)

                                              tCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCt                                          >  2:1302459/1‑93 (MQ=255)

                                               cGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCtt                                         <  2:552115/93‑1 (MQ=255)

                                                gTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTAATTATTATTGGCttt                                        >  1:621803/1‑93 (MQ=255)

                                                 tcttctTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTc                                       <  2:2755788/93‑1 (MQ=255)

                                                  cttcttTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTCa                                      >  2:2115875/1‑93 (MQ=255)

                                                   ttcttTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTCaa                                     >  1:7132310/1‑93 (MQ=255)

                                                    tcttTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTCaaa                                    <  2:7097691/93‑1 (MQ=255)

                                                     cttTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTCaaaa                                   <  2:6204819/93‑1 (MQ=255)

                                                     cttTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTCaaaa                                   >  2:5055335/1‑93 (MQ=255)

                                                          cGTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTCAAAATGTAg                              >  1:4881628/1‑93 (MQ=255)

                                                           gTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTCAAAATGTAgc                             >  2:4317298/1‑93 (MQ=255)

                                                           gTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTCAAAATGTAgc                             <  1:5433797/93‑1 (MQ=255)

                                                           gTAGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTCAAAATGTAgc                             >  1:2891244/1‑93 (MQ=255)

                                                             aGCCAGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTCAAAATGTAgcgc                           >  1:1017863/1‑93 (MQ=255)

                                                                 aGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTCAAAATGTAGCGCACGc                       <  2:1956571/93‑1 (MQ=255)

                                                                 aGACGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTCAAAATGTAGCGCACGc                       <  1:5672159/93‑1 (MQ=255)

                                                                   aCGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTCAAAATGTAGCGCACGCTg                     >  2:2421624/1‑93 (MQ=255)

                                                                   aCGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTCAAAATGTAGCGCACGCTg                     >  2:5163780/1‑93 (MQ=255)

                                                                    cGGCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTCAAAATGTAGCGCACGCTgc                    >  1:5536625/1‑93 (MQ=255)

                                                                     ggCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTCAAAATGTAGCGCACGCTgcg                   >  2:5622595/1‑93 (MQ=255)

                                                                      gCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTCAAAATGTAGCGCACGCTgcgc                  <  2:1572130/93‑1 (MQ=255)

                                                                      gCGATCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTCAAAATGTAGCGCACGCTgcgc                  <  1:6309562/93‑1 (MQ=255)

                                                                         aTCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTCAAAATGTAGCGCACGCTGCGCAGa               >  2:1848689/1‑93 (MQ=255)

                                                                         aTCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTCAAAATGTAGCGCACGCTGCGCAGa               >  2:7815060/1‑93 (MQ=255)

                                                                          tCTGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTCAAAATGTAGCGCACGCTGCGCAGaa              >  2:3172274/1‑93 (MQ=255)

                                                                            tGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTCAAAATGTAGCGCACGCTGCGCAGAata            <  1:1281889/93‑1 (MQ=255)

                                                                            tGTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTCAAAATGTAGCGCACGCTGCGCAGAata            >  1:5716322/1‑93 (MQ=255)

                                                                             gTTCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTCAAAATGTAGCGCACGCTGCGCAGAatat           <  1:3358463/93‑1 (MQ=255)

                                                                              ttCCGCGGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTCAAAATGTAGCGCACGCTGCGCAGAATATc          >  1:9944/1‑93 (MQ=255)

                                                                                   cgGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTCAAAATGTAGCGCACGCTGCGCAGAATATCGTTAc     <  2:5845427/93‑1 (MQ=255)

                                                                                   cgGTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTCAAAATGTAGCGCACGCTGCGCAGAATATCGTTAc     <  2:7201955/93‑1 (MQ=255)

                                                                                     gTGGCCTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGTCTTTCAAAATGTAGCGCACGCTGCGCAGAATATCGTTACtt   <  1:1664461/93‑1 (MQ=255)

GAGGAATGACGACAAAACAATGTCCGGACAACGCTGAAAACCCTCTTCGTCTTCTTTACGTAGCCAGACGGCGATCTGTTCCGCGGTGGCTTCGACATTACCCAGCGCCAGAATACGCACCAGGTCATTATTATTGGCTTTCAAAATGTAGCGCACGCTGCGCAGAATATCGTTACTTA  >  NC_000913/2210034‑2210212

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