BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,596,295	A→G	V332A (GTG→GCG)	b3460 ←	LivJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,596,295	0	A	G	99.6%	806.8 / ‑5.7	255	V332A (GTG→GCG)	b3460	LivJ
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (118/136); total (118/137)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 262 total aligned reads displayed.

GTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCACGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGCGCGGCGTAGGTG  >  NC_000913/3596205‑3596384

gtgccgtTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCcgcg                                                                                         >  1:3678895/1‑93 (MQ=255)

  gccgtTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGa                                                                                       <  2:4143375/93‑1 (MQ=255)

     gtTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGtt                                                                                    <  2:4226864/93‑1 (MQ=255)

     gtTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGtt                                                                                    >  2:1165929/1‑93 (MQ=255)

     gtTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGtt                                                                                    >  2:328273/1‑93 (MQ=255)

     gtTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGtt                                                                                    <  1:606252/93‑1 (MQ=255)

     gtTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGtt                                                                                    <  2:1967545/93‑1 (MQ=255)

      ttGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTc                                                                                   >  2:2927656/1‑93 (MQ=255)

       tGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCg                                                                                  >  2:684615/1‑93 (MQ=255)

          cGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCtt                                                                               <  1:3283239/93‑1 (MQ=255)

             gCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCa                                                                            <  2:3586226/93‑1 (MQ=255)

               cAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAgg                                                                          <  1:743832/93‑1 (MQ=255)

               cAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGGATCCGCGGAGTTCGCTTTCAgg                                                                          >  1:149309/1‑93 (MQ=255)

                aGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGt                                                                         >  2:4035652/1‑93 (MQ=255)

                 gTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTa                                                                        <  1:1052232/93‑1 (MQ=255)

                 gTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTa                                                                        <  2:2260562/93‑1 (MQ=255)

                    aaaTACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTAttt                                                                     <  1:1885872/93‑1 (MQ=255)

                    aaaTACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTAttt                                                                     >  2:2519641/1‑93 (MQ=255)

                    aaaTACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTAttt                                                                     <  1:1332838/93‑1 (MQ=255)

                    aaaTACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTAttt                                                                     <  2:1716163/93‑1 (MQ=255)

                       tACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGc                                                                  <  1:154342/93‑1 (MQ=255)

                         cGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGa                                                                >  2:4517915/1‑93 (MQ=255)

                          gCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGAt                                                               <  2:812988/93‑1 (MQ=255)

                             gAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTc                                                            <  2:3471199/93‑1 (MQ=255)

                             gAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTc                                                            <  1:3698747/93‑1 (MQ=255)

                             gAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTc                                                            <  2:3375718/93‑1 (MQ=255)

                              aaCTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCa                                                           <  2:2894285/93‑1 (MQ=255)

                              aaCTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCa                                                           >  1:2468577/1‑93 (MQ=255)

                                cTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGc                                                         <  1:888521/93‑1 (MQ=255)

                                cTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGc                                                         <  2:1418911/93‑1 (MQ=255)

                                   aaaGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCgg                                                      >  1:304496/1‑93 (MQ=255)

                                    aaGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGa                                                     >  1:2986531/1‑93 (MQ=255)

                                    aaGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGa                                                     <  2:3727437/93‑1 (MQ=255)

                                     aGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGAt                                                    <  2:2579302/93‑1 (MQ=255)

                                     aGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGAt                                                    <  2:1352398/93‑1 (MQ=255)

                                       ccTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCg                                                  >  2:907486/1‑93 (MQ=255)

                                        cTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGt                                                 <  1:4369917/93‑1 (MQ=255)

                                        cTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGt                                                 <  2:4303670/93‑1 (MQ=255)

                                         tttCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTc                                                >  1:2813305/1‑93 (MQ=255)

                                          ttCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCa                                               <  2:852224/93‑1 (MQ=255)

                                          ttCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCa                                               <  2:3615772/93‑1 (MQ=255)

                                           tCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAg                                              >  2:4103645/1‑93 (MQ=255)

                                            cAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGa                                             <  2:2211584/93‑1 (MQ=255)

                                            cAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGa                                             <  1:273186/93‑1 (MQ=255)

                                             aGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGAc                                            <  1:2328713/93‑1 (MQ=255)

                                             aGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGAc                                            >  2:2705250/1‑93 (MQ=255)

                                             aGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGAc                                            <  1:1567177/93‑1 (MQ=255)

                                              gATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGAATTTGGCGATTTCAGCCGGATCGTCAGACt                                           >  2:3862788/1‑93 (MQ=255)

                                               aTCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTg                                          >  2:724099/1‑93 (MQ=255)

                                                tCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGa                                         >  2:2011558/1‑93 (MQ=255)

                                                tCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGa                                         >  1:1294815/1‑93 (MQ=255)

                                                tCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGa                                         <  1:3563474/93‑1 (MQ=255)

                                                tCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGCTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGa                                         >  1:1995820/1‑93 (MQ=255)

                                                 cGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGAt                                        <  2:2679493/93‑1 (MQ=255)

                                                  gCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGAtt                                       <  1:3954431/93‑1 (MQ=255)

                                                    cTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCa                                     <  1:4141777/93‑1 (MQ=255)

                                                     tttCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGGCTGATTCAg                                    >  2:1478359/1‑93 (MQ=255)

                                                     tttCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAg                                    <  2:1852284/93‑1 (MQ=255)

                                                      ttCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAgg                                   <  1:4631600/93‑1 (MQ=255)

                                                       tctcATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGc                                  <  1:2093858/93‑1 (MQ=255)

                                                         tcATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGccc                                >  1:2156569/1‑93 (MQ=255)

                                                          cATCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGGCCg                               >  2:894926/1‑93 (MQ=255)

                                                            tCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGcc                             <  1:237153/93‑1 (MQ=255)

                                                            tCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGcc                             >  1:209858/1‑93 (MQ=255)

                                                            tCCCAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTCGCGATTTCAGCCGGATCGTCAGAATGATTCAcgcccgcc                             >  2:994924/1‑93 (MQ=255)

                                                             cccAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCt                            >  2:1070993/1‑93 (MQ=255)

                                                              ccAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTg                           >  1:4025836/1‑93 (MQ=255)

                                                              ccAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTg                           >  2:3007972/1‑93 (MQ=255)

                                                              ccAGGTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTg                           >  1:1675855/1‑93 (MQ=255)

                                                                 ggTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCaa                        >  1:2535006/1‑93 (MQ=255)

                                                                  gTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCaaa                       <  2:12850/93‑1 (MQ=255)

                                                                  gTCAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCaaa                       <  1:4528848/93‑1 (MQ=255)

                                                                    cAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGa                     <  1:1956041/93‑1 (MQ=255)

                                                                    cAGCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGa                     <  1:4608048/93‑1 (MQ=255)

                                                                      gCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGAtt                   >  2:3261770/1‑93 (MQ=255)

                                                                      gCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGAtt                   >  1:2690761/1‑93 (MQ=255)

                                                                      gCGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGAtt                   >  1:2873767/1‑93 (MQ=255)

                                                                       cGGTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTg                  <  2:1289831/93‑1 (MQ=255)

                                                                        ggTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGc                 <  2:2450587/93‑1 (MQ=255)

                                                                        ggTCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGc                 <  1:850243/93‑1 (MQ=255)

                                                                          tCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAg               <  1:1156343/93‑1 (MQ=255)

                                                                          tCCCATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAg               <  1:1795812/93‑1 (MQ=255)

                                                                           cccATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCGGGCCCGCCTGCAAAGATTGCAgc              >  1:2754987/1‑93 (MQ=255)

                                                                           cccATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAgc              <  2:3023950/93‑1 (MQ=255)

                                                                           cccATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAgc              <  2:4246262/93‑1 (MQ=255)

                                                                           cccATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAgc              <  1:66037/93‑1 (MQ=255)

                                                                           cccATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAgc              <  2:1552685/93‑1 (MQ=255)

                                                                             cATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAgcgc            <  1:4429329/93‑1 (MQ=255)

                                                                             cATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAgcgc            >  2:3518471/1‑93 (MQ=255)

                                                                             cATTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAgcgc            <  2:282754/93‑1 (MQ=255)

                                                                              aTTACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGcgcg           >  1:2187013/1‑93 (MQ=255)

                                                                               ttccGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGCgcgg          <  1:545279/90‑1 (MQ=255)

                                                                               ttACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGCgcgg          <  1:165107/93‑1 (MQ=255)

                                                                                tACGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGCgcggc         <  2:4470742/93‑1 (MQ=255)

                                                                                 aCGGTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGCgcggcg        >  2:531307/1‑93 (MQ=255)

                                                                                   ggTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGCGCGGCGTa      <  1:3201119/93‑1 (MQ=255)

                                                                                   ggTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGCGCGGCGTa      >  2:4427297/1‑93 (MQ=255)

                                                                                    gTATCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGCGCGGCGTAg     <  2:2800694/93‑1 (MQ=255)

                                                                                       tCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGCGCGGCGTAggtg  <  2:419891/93‑1 (MQ=255)

                                                                                       tCCGCGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGCGCGGCGTAggtg  >  2:1072281/1‑93 (MQ=255)

GTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCACGGAGTTCGCTTTCAGGTATTTGGCGATTTCAGCCGGATCGTCAGACTGATTCAGGCCCGCCTGCAAAGATTGCAGCGCGGCGTAGGTG  >  NC_000913/3596205‑3596384

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: `‑`