BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,515,312	T→C	M181V (ATG→GTG)	b2395 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,515,312	0	T	C	99.3%	471.8 / ‑4.6	147	M181V (ATG→GTG)	b2395	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (70/76); total (70/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.85e-01

Only 100 of 156 total aligned reads displayed.

CTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCATGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTA  >  NC_000913/2515222‑2515401

ctTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCAcgt                                                                                         >  1:1719801/1‑93 (MQ=255)

ctTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCAcgt                                                                                         <  1:4214575/93‑1 (MQ=255)

 tttGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCAcgtt                                                                                        >  1:4737236/1‑93 (MQ=255)

  ttGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCAcgttg                                                                                       <  2:733110/93‑1 (MQ=255)

      ccAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTaaa                                                                                   <  2:1222070/93‑1 (MQ=255)

       cAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAAt                                                                                  >  1:2817110/1‑93 (MQ=255)

         acgacgGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAAtaa                                                                                >  1:226677/1‑91 (MQ=255)

         acgacgGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCa                                                                                <  2:3472445/93‑1 (MQ=255)

          cgacgGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAg                                                                               >  2:2057515/1‑93 (MQ=255)

           gacgGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGc                                                                              >  2:4329475/1‑93 (MQ=255)

             cgGAGCGATAACCCTGCACCCCAATAAATGCGCACAGTGGTGAATTACAATCATGAGAGTGAGATAGTAGAAACGAACGTTGTAAATCAGcac                                                                            <  1:1002159/93‑1 (MQ=255)

              ggAGCGATATCCCTGCGCCACAATATCTGCGCCAAGTGGGCACTTAGAATCATGCTAGTGCGATAGTCGAAAAGCACGATGGAAATCAGcaca                                                                           >  2:440459/1‑93 (MQ=255)

              ggAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCCAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGcaca                                                                           >  1:3877832/1‑93 (MQ=255)

               gAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAg                                                                          >  2:2285599/1‑93 (MQ=255)

               gAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAg                                                                          >  2:4180279/1‑93 (MQ=255)

                aGCTATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGc                                                                         <  1:3505624/93‑1 (MQ=255)

                aGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGc                                                                         <  2:2502684/93‑1 (MQ=255)

                aGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGc                                                                         <  2:3069574/93‑1 (MQ=255)

                aGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGc                                                                         <  2:2004098/93‑1 (MQ=255)

                aGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGc                                                                         <  2:3697034/93‑1 (MQ=255)

                 gCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCg                                                                        <  2:2715266/93‑1 (MQ=255)

                 gCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCg                                                                        <  2:4049842/93‑1 (MQ=255)

                 gCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCg                                                                        <  2:2970199/93‑1 (MQ=255)

                   gATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGt                                                                      >  2:1478335/1‑93 (MQ=255)

                    atatCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTg                                                                     <  2:780690/93‑1 (MQ=255)

                     tatCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGa                                                                    >  2:4339564/1‑93 (MQ=255)

                     tatCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTCTAAATCAGCACAGCGGTGa                                                                    >  2:933703/1‑93 (MQ=255)

                      atCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGaa                                                                   <  2:886210/93‑1 (MQ=255)

                      atCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGaa                                                                   >  1:3497786/1‑93 (MQ=255)

                       tCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGaaa                                                                  <  1:3868247/93‑1 (MQ=255)

                        cccTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGaaaa                                                                 <  2:1762842/93‑1 (MQ=255)

                        cccTGCGCCACAATAGCTGCGCAAAGGTGTGACTTACAAGCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGaaaa                                                                 <  1:3082460/93‑1 (MQ=255)

                          cTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAgg                                                               >  1:1771579/1‑93 (MQ=255)

                           cgcgcCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGc                                                              <  1:3554082/92‑1 (MQ=255)

                             cgcCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCtt                                                            <  2:3191345/93‑1 (MQ=255)

                              gcCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTa                                                           >  1:718831/1‑93 (MQ=255)

                               ccACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAg                                                          <  2:330974/93‑1 (MQ=255)

                               ccACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAg                                                          >  1:3031302/1‑93 (MQ=255)

                                cacaATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGAACGTTGTAAATCAGCAAAGCGGTGAAAAGGCTTAga                                                         >  1:3551164/1‑92 (MQ=255)

                                 acaATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCa                                                        <  1:1526983/93‑1 (MQ=255)

                                 acaATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCa                                                        >  2:1556317/1‑93 (MQ=255)

                                  caATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCaa                                                       <  2:2641797/93‑1 (MQ=255)

                                  caATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCaa                                                       >  1:1964431/1‑93 (MQ=255)

                                     tatCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATc                                                    <  1:2204636/93‑1 (MQ=255)

                                     tatCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATc                                                    <  1:2452531/93‑1 (MQ=255)

                                         tGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACg                                                >  2:4137746/1‑93 (MQ=255)

                                          gcgcAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGa                                               <  1:2427364/93‑1 (MQ=255)

                                           cgcAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGAc                                              <  1:1033852/93‑1 (MQ=255)

                                           cgcAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGAc                                              <  2:3784096/93‑1 (MQ=255)

                                            gcAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGAcc                                             <  1:3841650/93‑1 (MQ=255)

                                             cAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCg                                            >  1:3282501/1‑93 (MQ=255)

                                                aGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGa                                         >  2:2223400/1‑93 (MQ=255)

                                                 gTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGaa                                        >  1:644586/1‑93 (MQ=255)

                                                   ggggACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGaaaa                                      <  1:2343962/93‑1 (MQ=255)

                                                   ggggACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGaaaa                                      >  1:3115914/1‑93 (MQ=255)

                                                    gggACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAAt                                     <  2:377724/93‑1 (MQ=255)

                                                    gggACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAAt                                     <  2:4690514/93‑1 (MQ=255)

                                                     ggACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAATTg                                    >  1:2358544/1‑93 (MQ=255)

                                                      gACTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATgg                                   >  2:3940835/1‑93 (MQ=255)

                                                       aCTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGc                                  >  2:3014021/1‑93 (MQ=255)

                                                        cTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCa                                 >  2:2401852/1‑93 (MQ=255)

                                                        cTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCa                                 <  1:3429954/93‑1 (MQ=255)

                                                        cTTACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCa                                 >  1:4250383/1‑93 (MQ=255)

                                                          tACAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATc                               >  2:1273828/1‑93 (MQ=255)

                                                           aCAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATcc                              >  1:2569693/1‑93 (MQ=255)

                                                           aCAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATcc                              <  1:3718821/93‑1 (MQ=255)

                                                           aCAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATcc                              <  2:516733/93‑1 (MQ=255)

                                                           aCAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATcc                              <  2:3144110/93‑1 (MQ=255)

                                                           aCAATCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATcc                              <  1:4017911/93‑1 (MQ=255)

                                                               tCATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCAt                          >  1:1395060/1‑93 (MQ=255)

                                                                cATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATc                         <  2:29361/93‑1 (MQ=255)

                                                                cATGCGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATc                         <  1:1613323/93‑1 (MQ=255)

                                                                    cGAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCg                     <  2:1804032/93‑1 (MQ=255)

                                                                     gAGTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGa                    >  2:1810964/1‑93 (MQ=255)

                                                                      aGTGAGATAGTAGAAAAGCACGTTGTAAATAAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGaa                   >  2:2215588/1‑93 (MQ=255)

                                                                       gTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGaaa                  >  1:866921/1‑93 (MQ=255)

                                                                       gTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGaaa                  >  2:3853284/1‑93 (MQ=255)

                                                                       gTGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGaaa                  >  1:2913352/1‑93 (MQ=255)

                                                                        tGAGATAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGaaaa                 <  1:1163099/93‑1 (MQ=255)

                                                                         gagaTAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGaaaaa                >  2:947079/1‑93 (MQ=255)

                                                                         gagaTAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGaaaaa                <  2:639018/93‑1 (MQ=255)

                                                                         gagaTAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGaaaaa                <  1:4506877/93‑1 (MQ=255)

                                                                         gagaTAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGaaaaa                >  2:2188679/1‑93 (MQ=255)

                                                                          agaTAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAg               >  2:369730/1‑93 (MQ=255)

                                                                           gaTAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGGTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAgg              >  1:1290866/1‑93 (MQ=255)

                                                                           gaTAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAgg              >  2:3672259/1‑93 (MQ=255)

                                                                            aTAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAggg             >  2:2477800/1‑92 (MQ=255)

                                                                            aTAGTAGAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGt             <  1:4054944/93‑1 (MQ=255)

                                                                            aTAGTAGAAAAGCAAGTTGTAAATCAGCACAGCGGTGAAAAGGCTTATCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGt             <  2:1682871/93‑1 (MQ=255)

                                                                             aagtaaAAAAGAACGTTGTAAATCAGCACAGCGGTAAAAAGGCTTAGCAAATCAACGCCCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTg            <  2:3285611/92‑1 (MQ=255)

                                                                              agtagAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTgg           >  1:1141591/1‑93 (MQ=255)

                                                                               gtagAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATTGCCGAAAAAGGTGGa          >  2:1033223/1‑93 (MQ=255)

                                                                               gtagAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGa          >  2:3103046/1‑93 (MQ=255)

                                                                                tagAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGat         <  1:4627499/93‑1 (MQ=255)

                                                                                tagAAAAGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGat         <  2:4514875/93‑1 (MQ=255)

                                                                                    aaaGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCt     >  1:1468137/1‑93 (MQ=255)

                                                                                     aaGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCtt    >  2:1310761/1‑93 (MQ=255)

                                                                                      aGCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCttt   <  2:3050968/93‑1 (MQ=255)

                                                                                       gCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTa  <  1:3739296/93‑1 (MQ=255)

                                                                                       gCACGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTa  <  1:2366505/93‑1 (MQ=255)

CTTTGCCCAACGACGGAGCGATATCCCTGCGCCACAATATCTGCGCAAAGTGGGGACTTACAATCATGCGAGTGAGATAGTAGAAAAGCATGTTGTAAATCAGCACAGCGGTGAAAAGGCTTAGCAAATCAACGACCGTGAAAATGGCATCCGCATCGCCGAAAAAGGTGGATATCTTTA  >  NC_000913/2515222‑2515401

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